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1.0  EXECUTIVE  SUMMARY 


Four  wave  mixing  (FWM)  performance  in  Na  vapor  was  evaluated  in  the 
strong  pumping  regime  (I»Iggt)  required  for  typical  applications  related 
to  point-to-point  optical  links.  Gain  (■  reflectivity),  imaging  fidelity, 
and  noise  crosstalk  were  measured  and  analyzed  under  several  conditions  to 
obtain  an  engineering  data  base  for  potential  tactical  applications  such  as 
communications,  surveillance,  and  guidance  systems. 

Results  include  >230%  reflectivity  and  15  line  pair/mm  imaging 
resolution,  both  the  highest  reported  to  date  with  cw  pumping  in  low 
density  Na  vapor.  The  observed  resolution  was  limited  by  the  pixel 
spacing  of  our  video  camera.  Ten  Ip/mm  resolution  was  obtained  after 
passing  the  input  image  through  a  65  A  aberrator.  Noise  crosstalk  was 
observable  only  when  the  combined  signal/noise  powers  were  large  enough  to 
saturate  the  phase  conjugate  return  signal  via  pump  depletion.  Under 
these  conditions,  a  degenerate,  coherent  noise  beam  3x  stronger  than  the 
signal  beam  produced  only  a  21%  reduction  in  the  intensity  of  the  phase 
conjugate  signal,  and  this  degree  of  crosstalk  decreased  rapidly  as  the 
pump/noise  angle  increased. 

The  experiments  were  performed  with  pump  and  probe  beams  that  were 
derived  from  a  single,  800  mW,  CW  ring  dye  laser  (Coherent  699-21)  oper¬ 
ating  narrowband  (—1  MHz)  or  broadband  (2GHz)  near  the  589  nm  D2  resonance 
line  in  Na.  Pump  intensities  ranged  from  60-130  w/cm^;  focused  beam 
diameters  of  0.2-1  mm  l/e^  produced  intensities  of  10^-10^  which 

power  broadened  the  Na  linewidths  to  30-100x  of  their  collision-free  width 
of  62.5  MHz.  An  on-line  Lamb  Dip  diagnostic  monitored  laser  frequency 
relative  to  the  resonance  line  with  50  MHz  accuracy.  Vapor  densities  in 
the  10^^-10^^/cm^  range  were  tested  in  quartz  and  pyrex  cells  with  path- 
lengths  varying  from  1  mm  to  10  cm.  Highest  reflectivity  was  obtained 
with  no  buffer  gas;  He  buffer  pressures  up  to  50  torr  were  used. 

The  reflectivity  of  DFWM  in  Na  vapor  was  measured  versus  frequency 
offset  from  resonance,  pump  intensity,  pump /probe  angle,  and  pump  band¬ 
width.  With  a  300  mW,  1  MHz  pump  two  reflectivity  peaks  were  observed  on 
the  Na  D2  line,  one  on  the  low  frequency  side  of  the  2-3  hyperfine 
transition  and  one  on  the  high  frequency  side  of  1-0  transition.  Their 
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FWHMs  were  -2.2  GHz  and  1.4  GHz  respectively;  the  2-3  reflectivity  peak 
was  about  20%  of  the  1-0  maximum.  Little  or  no  reflectivity  was  observed 
on  the  other  four  hyperfine  transitions  of  the  D2  line  because  optical 
pumping  depletes  the  ground  state  populations  of  these  lines.  The  strong 
linear  absorption  of  the  medium  (OqL  »  200)  eliminated  a  conjugate  returh 
in  a  2.5  GHz  frequency  range  that  overlaps  the  line  centers  of  the  2-3  and 
1-0  lines  which  are  separated  by  1.7  GHz.  An  increase  of  pump  power  from 
300  to  450  mW  increased  the  2-3  and  1-0  reflectivity  maxima  by  20%  and  12% 
respectively,  the  corresponding  bandwidths  went  up  by  18%  and  50%  respec¬ 
tively,  and  the  gap  between  the  two  peaks  narrowed  to  1.2  GHz.  The  latter 
result  is  due  to  increased  saturation  of  the  linear  absorption  in  that 
region.  A  maximum  reflectivity  of  230%  was  obtained  on  the  1-0  hyperfine 
transition  with  the  peak  occurring  0.9  GHz  on  the  high  frequency  side  of 
the  transition  line  center. 

Self  action  effects  were  observed  as  increases  in  the  divergences 
of  both  the  transmitted  pump  beams  and  the  phase  conjugate  return.  Self 
focusing,  which  occurs  on  the  high  frequency  side  of  a  transition, 
increases  the  effective  pump  intensity  on  the  1-0  line  and  is  the  reason 
the  reflectivity  there  is  higher  than  on  the  2-3  line.  The  low  frequency 
side  of  the  2-3  line  is  subject  to  self  defocusing,  which  reduces  the 
effective  pump  intensity  there.  As  expected,  the  2-3  line,  which  has  a 
larger  oscillator  strength  and  degeneracy  factor,  exhibits  the  higher 
reflectivity  at  lower  pump  Intensities  where  self  action  effects  are 
negligible. 

The  reflectivity  FWHM  versus  pump/probe  angle  was  measured  to  be  -2.2 
mrad  with  1/sin^  8  dependence  for  angles  greater  than  13  mrad.  Tliis 
dependence,  which  is  predicted  for  a  homogeneously  broadened  system, 
is  also  shown  to  be  expected  for  angles  <  Awp^/Awp,  the  ratio  of  the  power - 
broadened  homogeneous  width  to  the  inhomogeneous  width  (=*  0.2  at  our  test 
conditions).  With  2  GHz  bandwidth  pumps,  reflectivities  as  high  as  20% 
were  measured  and  the  FWHM  vs  pump/probe  angle  doubled  to  -  4.5  mrad. 
Simultaneous  degenerate  and  nondegenerate  FWM  interactions  are  believed  to 
contribute  to  the  broadband  response. 

An  active  tracking  demonstration  showed  automatic  tracking  over  a  50 
mrad  field  of  view  at  angular  scan  rates  in  excess  of  100  rad/s. 
tracking  error  was  less  than  our  detection  limit  of  1  /irad. 


Residual 


Noise  mechanisms  intrinsic  to  the  FWM  process  were  studied  by 
introducing  a  second  probe  (noise)  beam  into  the  phase  conjugation  region. 
Intensity  crosstalk,  which  occurs  when  the  presence  of  a  noise  beam 
changes  the  phase  conjugate  reflectivity  of  the  signal  beam,  was  observed 
and  characterized.  Saturation  of  the  nonlinearity  and/or  pump  depletioi 
produce  intensity  crosstalk  when  high  probe  intensities  are  used  (Ip^obe  - 
0.1  IpuQp) .  A  maximum  21%  reduction  in  the  phase  conjugate  signal  was 
observed  at  a  nolse/slgnal  intensity  ratio  of  3:1  and  noise/signal  angular 
separation  of  13  mrad.  Variation  of  signal,  noise,  and  pump  polarizations 
had  no  affect  on  the  degree  of  crosstalk.  Reflected  phase  conjugate  power 
saturates  as  the  probe/pump  intensity  ratio  is  increased  from  0.02  to  0.3; 
this  saturation  curve  can  be  used  to  predict  intensity  crosstalk  levels. 
Our  analysis  indicates  there  is  no  mechanism  for  phase  crosstalk  between  a 
signal  and  noise  beam,  and  none  was  oberved. 

Spatial  crosstalk  is  defined  as  the  transfer  of  transverse  intensity 
patterns  from  a  noise  beam  to  the  phase  conjugate  signal  wave.  This  was 
studied  by  imprinting  a  2.8  line  pair/mm,  three  bar  image  from  an  Air 
Force  resolution  target  onto  a  6  mW  noise  beam,  and  looking  for  that 
image  pattern  in  the  phase  conjugate  return  of  a  1  mU  signal;  none  was 
observed.  Intensity  crosstalk  reduced  the  contrast  in  the  phase  conjugate 
signal  image,  a  three  bar  image  orthogonal  to  that  in  the  noise  beam,  by 
46%,  but  there  was  no  reduction  in  the  sharpness  of  the  edges  and  there 
was  no  evidence  of  an  orthogonal  pattern  overlaid  onto  the  original 
signal.  The  noise  and  signal  inputs  were  focused  into  the  Na  cell  so  the 
intensity  patterns  in  the  conjugation  region  were  the  Fourier  transforms 
of  the  orthogonal  bar  Images.  Only  the  central  lobes  of  the  signal 
and  noise  inputs,  which  contain  the  low  spatial  frequency  components  of 
the  original  Images,  were  coincident.  Thus  the  crosstalk  led  to  contrast 
reduction  with  no  loss  of  edge  definition.  Our  image  fidelity  for  these 
experiments  was  limited  by  the  pixel  resolution  of  our  video  camera  to 
about  15  Ip/mm;  pump  beam  1/e^  diameters  were  about  1100  urn.  These 
measurements  were  made  at  pump  intensities  where  self  action  effects  were 
negligible  and  phase  conjugate  reflectivities  were  10-25%. 

The  operating  region  of  interest  for  field  applications  involves 
strong  pumping  to  obtain  high  reflectivity  of  weak  signals.  In  general 
the  theory  in  this  regime  is  not  amenable  to  available  analytical  or 


numerical  treatments.  Therefore  an  analytical  framework  for  understanding 
FVm  In  a  strongly  pumped,  Doppler-broadened,  saturable  absorber  was 
established.  This  was  supplemented  with  development  of  simple  models  for 
homogeneously  broadened,  nonresonant  and  resonant  systems  that  can  later 
be  upgraded  to  Incorporate  Inhomogeneous  broadening  and  wave  optics 
considerations.  For  high  reflectivity  cases  where  there  Is  strong  linear 
absorption  and  strong  optical  saturation  we  show  that  the  laser  frequency 
for  maximum  nonlinearity  Is  shifted  off  resonance  by  an  amount  dependent 
on  optical  Intensity  and  on  the  ratio  of  signal  and  piomp  beam  Intensities. 

The  frequency  offset  for  maximum  FWM  reflectivity  In  a  strongly 
pumped  saturable  absorber  Is  explained  as  follows.  The  Intensity  contrast 
In  the  alternating  spatial  regions  where  the  Input  pump  and  signal  waves 
Interfere  produces  spatial  variations  In  the  nonlinear  susceptibility  x  of 
the  medlvim.  Fringe  contrast  Is  high  when  beams  of  equal  Intensity 
Interfere,  but  the  dark  regions  remain  strongly  Illuminated  when  beams  of 
unequal  Intensity  Interfere.  In  a  strongly  saturated  system  with  equal 
Intensity  beams,  high  fringe  contrast  produces  large  spatial  variations  In 
X  that  are  maximized  near  resonance.  When  weak  signal  beams  are  combined 
with  a  strong  pump  Intensity  that  Is  >  Igat*  fJ^inge  contrast  Is  low  and 
both  the  light  and  dark  regions  remain  strongly  saturated  near  resonance 
with  a  result  that  x  variations  are  also  small.  Then  maximum  nonlinearity 
Is  obtained  by  moving  off  resonance  to  a  frequency  where  the  effective 
I/Isat  “  because  this  Is  the  regime  where  x  Is  most  sensitive  to  small 
changes  In  Illumination  Intensity. 

Doppler  broadened  media  can  be  understood  by  estimating  the  effects 
of  experimental  parameters  on  the  population  of  the  velocity  group  In 
which  the  Interference  grating  Is  written.  Key  parameters  are  frequency 
offset,  relative  pump  and  signal  angles,  and  optical  Intensities.  A 
procedure  for  Identifying  the  velocity  group  that  Is  simultaneously 
resonant  with  the  Input  optical  waves  has  been  outlined. 

The  optical  field  Interacts  resonantly  only  with  absorbers  that  have 
velocity  components  that  shift  their  (stationary)  resonance  frequency  to 
that  of  the  optical  field.  In  the  standard,  degenerate  FWM  geometry  with 
weak,  counterpropagatlng  pumps;  maximum  reflectivity  occurs  on  resonance 
and  only  the  v-0  velocity  group  can  contribute  to  the  nonlinearity.  When 
strong  pumping  Is  used,  the  fringe  contrast  considerations  mentioned  above 
move  the  maximum  nonlinearity  to  a  nonresonant  frequency.  Now  the  various 
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input  waves,  which  typically  enter  at  different  angles,  optimally  interact 
with  velocity  groups  moving  in  different  directions.  One  result  is  that 
maximum  reflectivity  can  occur  with  nondegenerate  rather  than  degenerate 
FWM  in  this  case  if  the  frequency  of  one  of  the  counterpropagating  pumps 
is  shifted  by  an  amount  equal  to  when  the  frequencies  of  the  other 
pump  and  the  signal  are  +bu  from  resonance. 

In  general,  FWM  performance  tends  to  optimize  when  frequency  offsets 
and  angles  are  set  so  that  the  velocity  group  involved  is  near  the  thermal 
velocity.  The  situation  becomes  very  complex,  however,  when  the  details 
of  Na  hyperfine  spectroscopy,  hyperfine  optical  pumping,  power  broadening, 
and  fringe  contrast  considerations  are  included. 

The  Na  hyperfine  transitions  are  closely  spaced  relative  to  the 
Doppler  line  width,  so  each  transition  can  be  interacting  resonantly 
through  separate  velocity  groups.  Under  low  intensity  illumination  the 
hyperfine  separations  are  greater  than  the  homogeneous  linewidths,  but  in 
our  strong  pumping  regime  the  homogeneous  widths  become  greater  than  the 
hyperfine  spaclngs.  The  increased  homogeneous  bandwidth  increases  the 
rate  at  which  thermal  washout  of  the  Interference  grating  occurs,  because 
it  permits  Interaction  with  more  velocity  groups  that  are  moving  perpendic¬ 
ular  to  the  grating  lines.  In  the  pure  Doppler  limit  only  velocities 
parallel  to  the  Interference  grating  lines  are  permitted.  Finally  thermal 
diffusion  conteracts  the  optical  pumping  that  occurs  in  4  out  of  6  of  the 
Na  D2  hyperfine  transitions.  The  preceding  phenomena  have  variable, 
opposing  Impacts  on  the  net  nonlinearity  obseirved.  Our  physical  picture, 
which  provides  an  outline  for  qualitatively  understanding  these  effects , 
can  serve  as  the  basis  for  further  development  of  a  comprehensive  model. 

The  contract  SOW  called  for  experimental  measurements  that  address 
technical  issues  relating  to  phase  conjugated  optical  links.  No  analysis 
beyond  that  required  to  empirically  relate  results  to  possible  application 
scenarios  was  specified.  The  required  reflectivity,  crosstalk,  and 
fidelity  experiments  that  were  performed  are  presented  in  this  report. 
During  the  program  it  was  determined  that  analysis  of  FWM  physics  ’^eyond 
that  required  by  the  contract  was  deslreable.  In  particular,  noise 
crosstalk  mechanisms  were  investigated  on  TRW  IR&D,  and  this  theoretical 
description  is  included  in  Section  5  because  it  substantiates  our  conclu¬ 
sion  that  crosstalk  will  not  seriously  impact  FWM  applications. 
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2.0  INTRODUCTION 


This  final  report  describes  results  obtained  during  the  first  phase 
of  the  Nonlinear  Optics  Technology  program,  a  one  year  effort  to  invest¬ 
igate  phase  conjugation  by  four  wave  mixing  (FWM)  for  enhancing  perform¬ 
ance  of  electro -optical  systems  in  tactical  applications.  There  were  two 
primary  tasks: 

1.  An  experimental  investigation  of  the  characteristics  of  FWM  in 
the  high  reflectivity  regime  required  by  many  field  applica¬ 
tions.  Sodium  vapor  was  chosen  as  the  FWM  medivim.  Studied  were 
the  spectral  response  of  reflectivity  for  strong  CW  pumping 
conditions,  field  of  view,  noise  crosstalk,  image  fidelity,  and 
the  impact  of  finite  pump  laser  bandwidth. 

2.  An  analytical  modeling  task  to  develop  the  understanding 
necessary  for  projecting  potential  mission  performance  based  on 
laboratory  measurements.  Numerical  one  dimensional  and  two 
dimensional  models,  as  well  as  a  qualitative  physical  model, 
were  developed  to  explain  the  observed  behavior. 

In  addition  to  these  tasks,  some  effort  was  directed  to  identifying 
potential  tactical  missions  that  could  benefit  from  the  application  of 
phase -conjugation  FWM  technology.  Mission  factors  that  could  be 
substantially  impacted  by  FWM  were  identified.  These  include  optics 
jitter,  image  quality,  optical  link  tracking,  laser  power  requirements, 
and  others.  Several  potential  mission  applications  were  identified,  and 
one  of  the  most  promising,  tactical  optical  communications  links,  was 
briefly  evaluated. 

This  section  of  the  report  summarizes  the  overall  program  objectives 
and  guiding  considerations  and  introduces  the  subject  of  nonlinear  phase 
conjugation  via  FWM  interactions.  Section  3  presents  a  general  discussion 
of  FWM  in  sodium  vapor,  and  the  remaining  sections  describe  specific 
measurements  and  analyses  that  were  performed. 

2.1  PROGRAM  APPROACH/GOALS 

FWM  in  sodium  vapor  has  been  the  subject  of  many  investigations  over 
the  last  ten  years.  While  quite  elegant,  most  of  this  work  has 
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concentrated  on  the  basic  spectroscopy  of  the  FUM  process  and  its 
development  as  a  nonlinear  tool  under  experimental  conditions  tractable 
for  theoretical  description.  The  focus  of  the  work  reported  here  has  been 
to  Investigate  FVM  parameters  that  are  critical  for  effective  mission 
performance,  under  experimental  conditions  that  are  optimized  for  best 
mlsslon-orlented  performance  rather  than  analytical  convenience. 

Parameters  that  were  studied  Included  spectral  response,  field  of  view, 
crosstalk  between  different  Input  beams  to  a  FVM,  and  conjugation 
fidelity.  An  effort  was  made  to  acquire  data  on  these  parameters  In  the 
experimental  regime  of  high  FUM  phase  conjugate  reflectivity,  this  regime 
being  a  particularly  Important  one  for  practical  Implementation  of  this 
technology.  It  was  a  program  goal  to  understand  better  the  dependence  of 
these  mlsslon-orlented  performance  characteristics  on  experimental 
conditions  In  the  high  reflectivity  regime. 

2.2  PRACTICAL  CONSIDERATIONS/REQUIREMENTS 

High  reflectivity  is  usually  achieved  using  Intense  peak  pump  powers 
available  from  pulsed  laser  sources.  However,  use  of  very  short  pulse 
pvimps  Is  not  a  pragmatic  approach  when  viewed  in  light  of  many  mission 
applications  involving  detection  and  conjugation  of  a  return  reflection 
from  a  remote  target.  Use  of  short  pulse  pumps  requires  a  veiry  accurate 
knowledge  of  range  to  target  if  the  pump  pulses  are  to  coincide  in  time 
with  the  return  pulse  in  the  FUM.  For  this  reason,  practical  application 
of  FUM  technology  for  remote  sounding  applications  requires  a  FUM  medium 
that  can  provide  high  conjugate  reflectivity  at  the  low  to  moderate  pump 
powers  available  from  long-pulse  or  CW  pump  sources .  Such  materials  are 
also  necessary  to  reduce  the  size  and  weight  of  phase  conjugating 
electro-optical  systems .  Sodium  vapor  was  employed  as  the  FVJM  medium  in 
our  experiments  because  it  has  a  strong  nonlinearity  that  can  lead  to  high 
reflectivity  for  modest  pump  Intensities  in  the  visible.  This  made  it 
possible  to  perform  the  experiments  with  relatively  low  Intensities 
available  from  CW  narrowband  dye  laser  sources . 

2.3  OPTICAL  PHASE  CONJUGATION  BY  FWM 

A  number  of  different  nonlinear  optical  techniques  have  been 
developed  for  optical  phase  conjugation,  the  process  in  which  an  optical 


wave  has  Its  propagation  direction  and  relative  phase  reversed  so  that  it 
exactly  retraces  its  path.  Four  wave  mixing  is  one  approach  particularly 
well-suited  for  missions  where  the  return  or  input  signal  is  weak  because 
the  process  has  no  threshold.  Reference  2-1  provides  an  extensive  review 
of  phase  conjugation,  including  discxissions  of  the  FWM  process. 


2.3.1  General  Principles  of  FWM 

Phase  conjugation  by  FWM  is  typically  performed  in  the  backward,  or 
counterpropagating  pump  geometry  shown  in  Figure  2-1.  Three  input  waves 
incident  on  a  material  with  a  nonlinear  response  are  mixed  to  form  a 
fourth  wave  that  is  radiated  from  the  medium  with  a  frequency  and 
direction  determined  by  the  frequencies  and  propagation  directions  of  the 
three  input  waves.  The  nonlinear  polarization  that  gives  rise  to  this 
fourth  wave  is 


p(3)  -  v<3)e3  2-1 

NL  ^  T 

where  E.J.  is  the  total  electric  field  present  in  the  nonlinear  medium  and 
is  the  third- order  nonlinear  susceptibility  of  the  medium.  Given 
input  waves  of  the  form 

-  1/2  exp[i(wjt  -  kj^z)]  +  complex  conjugate  2-2 

there  are  many  terms  in  the  third-order  nonlinear  polarization,  several  of 
which  have  the  form 
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If  Ej^  and  E2  are  phase  conjugates  of  one  another,  as  they  are  for  the 
two  counterpropagating  plane  wave  pumps  in  Figure  2-1,  then  this 
polarization  gives  rise  to  a  field  propagating  in  the  -k^  direction  that 
is  the  phase  conjugate  of  the  input  field  defined  by  Equation  2-2.  The 
interaction  can  be  degenerate,  i.e.  all  fieldi.  have  the  same  frequency 
w,  or  nondegenerate,  with  the  frequency  of  the  fourth  field  being  the 
sum  of  the  two  pump  fields  minus  the  input  signal  field  as  shown  in  the 
first  term  of  the  exponential  in  Equation  2-3.  The  dominant  four  wave 
mixing  process  out  of  all  of  the  terms  in  the  third- order  nonlinear 
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Figure  2-1.  Counterpropagating  pump  (EfEb)  geometry  for 
phase  conjugating  an  input  probe  beam  (Ep) 
using  FUM. 


polarization  is  determined  by  the  wavevector  of  the  generated  wave  defined 
by  the  second  factor  in  the  exponential.  This  wavevector  of  the  generated 
fourth  wave  needs  to  meet  the  phase -matching  condition  that  k^-nw^/c 
for  efficient  transfer  of  energy  from  the  pump  fields  to  the  input  signal 
wave  and  its  conjugate. 


2.3.2 


The  technical  data  base  on  the  performance  of  FWM  phase  conjugators 
under  conditions  pertinent  to  field  applications  was  inadequate  for  even 
simple  systems  engineering  projections.  High  power,  high  reflectivity 
data  exists  for  nanosecond  pulse  lengths,  and  low  power  spectroscopic  data 
is  available,  but  little  is  known  about  long  pulse  or  CW  behavior  of  media 
with  large  nonlinear  susceptibilities  in  the  high  reflectivity  regime. 

This  program  collected  data  needed  for  further  engineering  of  FWM 
phase -conjugated  electro-optical  systems. 
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3.0  FWM  IN  SODIUM  VAPOR 


To  establish  an  analytical  framework  for  Interpreting  the  experi¬ 
mental  results  of  Sections  4,  5,  and  6,  this  section  presents  a  discussion 
of  the  physics  of  four-wave  Interactions  In  sodium  vapor.  Instead  of  a 
purely  theoretical  treatment  with  limited  applicability  to  the  strong 
saturation  regime  of  our  experiments ,  Section  3 . 1  concentrates  on 
developing  a  physical  model  of  the  Important  processes  occur Ing  In 
Doppler -broadened  media  when  exposed  to  a  combination  of  weak  and  strong 
optical  fields.  Section  3.2  extends  this  discussion  to  the  specific  case 
of  the  sodium  D2  line.  Section  3.3  outlines  the  approach  taken  for 
analyzing  experimental  results  In  light  of  the  saturation  regime  In  which 
the  data  was  obtained. 

3.1  GENERAL  DESCRIPTION 


Nonlinear  optical  Interactions  In  Inhomogeneously-broadened  saturable 
absorbers  are  complex  experiments  Involving  only  selected  subgroups  of  the 
total  available  absorber  number  density  and  a  variety  of  frequency 
dependent  saturation  effects.  Current  analytical  solutions  are  valid  only 
In  the  limit  of  low  optical  Intensity.  In  this  section,  we  develop  a 
model  for  understanding  FWM  Interactions  In  a  strongly  pumped,  Doppler- 
broadened  saturable  absorber.  For  the  case  of  high  FWM  reflectivity, 
where  there  Is  strong  linear  absorption  and  strong  optical  saturation,  we 
show  that  the  expected  results  are  substantially  different  than  those 
predicted  by  analytical  models  In  either  the  homogeneous  or  low-power 
Doppler  limits.  In  this  general  discussion.  It  Is  shown  that: 

1.  the  maximum  nonlinearity  Is  shifted  away  from  the  laser  frequency 
with  a  shift  that  depends  on  the  optical  Intensity  and  the  ratio  of 
power  In  the  signal  and  pump  beams. 

2.  much  of  FWM  performance  In  Doppler-broadened  media  can  be 
understood  by  estimating  the  effects  of  experimental  parameters  on 
the  population  of  the  velocity  group  In  which  the  Interference 
grating  Is  written. 

3.  deleterious  thermal  motion  effects,  such  as  washout  of  the 
Interference  grating,  can  be  aggravated  as  the  homogeneous  width 


Increases.  Such  Increases  are  expected  due  to  power  broadening 
accompanying  the  strong  saturation. 

3.1.1  FUM  in  Saturable  Absorbers 


An  important  class  of  nonlinear  optical  materials  are  saturable 
absorbers.  A  strong  optical  nonlinear  response  is  possible  at  low  optical 
irradiance  because  the  refractive  index  can  be  changed  significantly  near 
an  optical  transition  due  to  saturation  of  the  absorption. 

A  physical  basis  for  understanding  enhancement  near  optical 
transitions  can  be  obtained  by  considering  the  behavior  of  the 
susceptibility,  x>  near  an  optical  transition.  It  is  exactly  this 
susceptibility  that  determines  the  polarization  of  the  medium, 


P  -  XE  -  x^^^E  +  x^^^E^  +  . . .  (3-1) 

where  the  nonlinear  Interaction  is  due  to  higher  orders  of  x  and  the  E 
field,  and  the  refractive  index  of  the  medium  is  proportional  to  x-  In 
general,  the  susceptibility  has  a  real  part,  x' .  giving  rise  to  a 
polarization  component  in  phase  with  E,  and  an  imaginary  component,  x" < 
which  is  out  of  phase.  It  is  this  out  of  phase  component  that  is 
responsible  for  energy  exchange  between  the  fields  and  the  medium,  i.e., 
is  the  absorptive  contribution  to  the  total  susceptibility  and  refractive 
index.  The  phase  shift  arises  from  the  finite  response  time  associated 
with  the  transition. 

Both  components  of  x  saturate  in  homogeneously  broadened 
transitions.  In  the  approximation  of  stationary  atoms,  the  dispersive  and 
absorptive  components  of  x  saturate  according  to 
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where  Aw  is  the  laser  detuning  from  line  center,  r  is  the  homogeneous 
line  halfwidth,  I®  transition  moment,  AN  is  the  population 


Is  the 


difference  between  the  two  levels  involved,  and 
characteristic  saturation  intensity.  Figure  3-1  shows  these  components 
near  an  optical  resonance  for  three  optical  intensities  (0,  1,  and  100 
times  the  saturation  intensity,  Ig^t^ •  illustrating  that  a  substantial 
intensity  dependence  of  the  refractive  index  can  be  achieved  for 
frequencies  close  to  the  optical  resonance.  Note  that  the  vertical  scale 
for  I/Ig^^-100  has  been  multiplied  by  100,  and  that  the  absorptive 
component  is  essentially  completely  saturated  for  many  homogeneous 
halfwidths,  F.  An  important  feature  in  understanding  FWM  in  saturable 
absorbers  is  that  the  absorptive  component,  x".  saturates  more  rapidly 
than  the  dispersive  component.  In  fact,  for  media  in  the 
Doppler -broadened  limit,  where  the  power-broadened  homogeneous  width, 

AWp^,  is  much  less  than  the  Doppler  width,  Aw^,  only  the 
absorptive  x*  saturates  appreciably.  The  dispersive  component  (x' )  is 
insensitive  to  field  intensity  because  it  consists  of  many  overlapping 
tails  of  dispersion  curves  from  the  many  homogeneous  subcomponents  of  the 
line . 

It  is  the  field-dependent  saturation  characteristic  of  the 
susceptibility  that  gives  rise  to  the  optical  nonlinearity.  Equation  3-2 
illustrates  additional  features  of  the  nonlinear  process  when  viewed  in 
light  of  particular  nonlinear  experiments.  In  the  case  of  four  wave 
mixing,  where  the  interference  of  the  input  waves  generates  a  spatial 
variation  of  optical  intensity,  the  nonlinearity  depends  on  spatial 
differences  of  the  susceptibility.  From  Equation  3-2  and  Figure  3-1,  it 
can  be  seen  that  the  maximum  difference  in  x  between  two  spatial  regions 
will  occur  if  one  region  is  completely  dark  and  the  other  is  illuminated, 
a  situation  that  will  occur  only  if  the  interfering  beams  have  the  same 
intensity.  In  this  case  of  equal  pump  and  signal  intensities,  the  maximum 
"contrast"  or  spatial  variation  in  x  occurs  close  to  the  resonance 
transition  frequency,  changing  from  A«  -  0  at  low  optical  intensity  to 
Aw  -  r  when  the  intensity  is  high  enough  to  strongly  saturate  x” • 

This  is  because  the  susceptibility  is  saturated  and  close  to  zero  for  all 
frequencies  in  the  light  regions,  but  is  nonzero  close  to  the  transition 
frequency  in  the  dark  regions.  Thus  for  input  beams  with  similar 
intensities,  the  maximum  nonlinearity  is  within  one  homogeneous  width  of 
I  the  transition  frequency.  This  will  also  be  true  for  input  Intensities 

I 

I 

I 

> 

I 


that  are  small  or  comparable  to  irrespective  of  the  relative  beam 

intensities.  It  should  be  noted  that  this  discussion  ignores  the  effects 
of  residual  absorption,  so  that  the  maximum  nonlinear  response  may  occur 
at  a  different  frequency  than  where  the  contrast  in  x  is  largest. 

However,  if  strong,  unequal  input  beams  are  involved,  the  situation 
is  quite  different.  In  such  a  case,  both  "dark"  and  "light"  regions  are 
strongly  saturated.  Near  the  transition  frequency,  the  two  regions 
saturate  identically  and  there  is  no  nonlinearity  "contrast" ;  only  for 
frequencies  detuned  on  the  order  of  (1/  homogeneous 

linewidths  from  resonance  will  the  two  regions  begin  to  saturate 
differently.  Thus  the  best  nonlinearity  is  achieved  for  substantial 
detunings  from  resonace  under  conditions  of  high  optical  intensity  with 
low  contrast  between  light  and  dark  regions.  It  is  this  regime  that  is  of 
particular  interest  to  this  program  because  we  are  interested  in  weak 
signal  beams  mixing  with  very  strong  ptunp  beams.  Consideration  of  this 
aspe' t  of  strong  saturation  in  both  light  and  dark  regions  is  particularly 
important  in  understanding  FWM  in  Doppler>broadened  saturable  absorbers 
such  as  sodium. 

3.1.2  FWM  in  Doppler -Broadened  Saturable  Absorbers 

In  a  homogeneously  broadened  transition,  the  resonance  condition  is 
given  by 


^  2r  (3-3) 

i.e.,  the  optical  field  frequency  must  be  within  a  homogeneous 
linewidth  of  the  isolated,  stationary  absorber  transition  frequency 
wq.  Only  part  of  the  absorbers  in  a  Doppler-broadened  transition  can 
interact  with  the  optical  fields  because  of  the  Doppler  shifts  induced  by 
motion  of  the  absorbers.  In  this  case,  the  resonance  condition  now 
becomes 


(3-4) 


where  kj^  is  the  optical  wavevector  and  v  is  the  absorber  velocity.  The 
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vector  dot  product  kL*^  velocity  component  along  the  optical 

field  propagation  direction  which  gives  rise  to  the  Doppler  shift  of  the 
resonance  frequency  wq.  Thus,  the  optical  field  interacts  resonantly 
only  with  absorbers  that  have  velocity  components  that  shift  their 
(stationary)  resonance  frequency  to  that  of  the  optical  field. 

The  consequences  of  this  selective  velocity  interaction  in  FWM 
experiments  are  several.  Consider  the  standard  counterpropagating  pump 
beam  FWM  geometry  illustrated  in  Figure  3-2.  The  forward  pump,  Ej, 
interacts  resonantly  with  absorbers  with  velocity  v-(w£-WQ)/k£. 

The  signal  beam,  E^,  which  intersects  E^  at  a  small  angle 

must  Interact  with  the  same  group  of  absorbers  if  an  interference 

grating  is  to  be  produced;  this  requires  that 

|«g-a)f|  <  2r.  (3-5) 

Since  the  maximum  contrast  in  x  is  close  to  the  Doppler- shifted 
resonance  frequency  of  the  absorbers  simultaneously  interacting  with  E^ 
and  Eg,  the  efficient  readout  of  the  spatial  grating  requires  that  Ej^ 
be  at  this  frequency  where  the  x  contrast  is  largest.  For  low 
intensities,  this  is  not  at  W£  because  the  back  pump  is 
counterpropagating  to  Ej.  Hence  the  absorbers  interacting  with  E£  and 
Eg  are  only  resonant  with  E|j  if 

|w£+(k£/kjj)c^-(l+(k£/kjj))wQ|  <  2r  (3-6) 

In  the  case  of  degenerate  four  wave  mixing  (DFWM) ,  the  above  resonance 
condition  can  only  be  satisfied  for  -  wq  (or  v-0)  .  Figure  3-3 
illustrates  this  for  DFWM,  showing  the  velocity  groups  that  the  different 
fields  interact  with.  Figure  3-3  also  shows  x'  frequency  space 
showing  the  offset  of  the  back  pump  from  the  frequency  where  x'  is 
largest.  This  illustrates  that  DFWM  can  be  performed  with  the  laser 
detuned  from  wq,  but  that  x  is  sampled  off- resonance  by  2k*v, 
resulting  in  a  weak  nonlinearity. 

In  fact  the  very  low  power,  low  Na  density  experiments  reported  in 
Reference  3-1  show  DFWM  reflectivity  exactly  on  the  stationary  resonance 
(v-0)  with  a  linewidth  substantially  narrower  than  the  Doppler  width.  This 
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is  as  expected  since  under  these  conditions  the  deleterious  effects  of  the 
small  absorption  at  line  center  (ug,  v~0)  are  unimportant  and  the 
nonlinearity  Is  so  small  that  It  must  be  sampled  on- resonance  to  achieve 
any  DFWM  reflectivity.  Humphrey,  et.  al.  (Reference  3-2)  obtain  similar 
results  with  llnewldths  approaching  the  homogeneous  width  of  Na  In  the 
collisionless  regime. 

The  situation  under  high  reflectivity  conditions  of  strong  linear 
absorption  and  high  pump  intensity  is  quite  different.  When  there  is 
strong  linear  absorption  near  tog,  higher  reflectivity  is  actually 
achieved  for  off  resonance  operation.  It  Is  better  to  reduce  the 
absorptive  loss  by  detuning  the  laser  frequency  and  "read"  the 
nonlinearity  nonresonantly .  This  is  discussed  in  detail  in  Section  4.3.1. 

Another  consideration  in  Doppler-broadened  saturable  absorbers  is  the 
fraction  of  absorbers  participating  in  the  experiment.  Figure  3-4 
illustrates  that  absorbers  with  velocities  within  a  range  F/k^  of 
V“Aw£/k£,  where  Aw£  is  (<)£-cog,  interact  with  E£.  These 

absorbers  define  a  disc  in  3 -dimensional  velocity  space  with  a  thickness 
along  the  direction  of  optical  propagation  of  r/k£,  and  a  radial 
distribution  given  by  the  Maxwellian  velocity  distribution  in  the 
directions  orthogonal  to  the  propagation  axis.  A  similar  disc  defines 
absorbers  interacting  with  Eg.  It  is  absorbers  in  the  velocity  volume 
defined  by  the  intersection  of  these  two  discs  that  participate  resonantly 
in  establishing  the  spatial  interference  grating  between  E£  and  Eg . 

The  nvunber  of  absorbers  in  the  intersecting  volume  is  determined  by 
their  absolute  velocity,  i.e.  the  length  of  the  vectors  in  Figure  3-4  that 
originate  at  v-0  and  terminate  within  the  selected  volume.  The  magnitude 
of  these  velocities  are,  in  turn,  determined  by  the  laser  detuning  from 
resonance.  Aw,  and  the  angle  between  the  two  fields  in  question.  The 
fraction  of  absorbers  with  appropriate  velocities,  can  then 

be  determined  from  the  Maxwell  velocity  distribution 

Nex/Ntot  -  (M/2xRT)3/2^expI-«('^x^+^y^+^z^)/2RT]dv^dVydv^  (3 

where  the  integration  is  over  the  intersecting  volume.  This  fraction  is  a 
key  parameter  in  determining  the  efficiency  of  the  nonlinear  process,  and 
so  considerable  insight  into  the  dependence  of  the  FWM  signal  on  various 
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parameters  can  be  gained  by  considering  the  impact  of  experimental 
conditions  both  on  the  size  of  the  intersection  volume  and  on  the 
magnitude  of  velocities  contained  within  this  volume. 

As  an  illustration,  consider  Figure  3-5  which  shows  in  two  dimensions 
the  intersecting  volumes  for  degenerate  optical  fields  crossing  at 
different  angles.  The  shaded  diamond- shaped  areas  indicate  the 
populations  which  have  the  proper  velocities  to  simultaneously  interact 
with  the  signal  and  forward  or  backward  pump  fields.  Except  in  the  case 
of  exactly  counterpropagating  signal  and  pump  fields,  there  is  always  a 
velocity  group  that  is  resonant  with  both  fields  regardless  of  Au> 
(counterpropagating  fields  require  A(<>-0)  .  However,  as  the  crossing 
angle  increases,  the  required  velocity  increases  from  a  minimum  of  Aw/k 
to  ®.  This  results  in  a  rapid  dropoff  in  absorber  population 
participating  resonantly  in  the  experiment,  particularly  if  Aw  is 
comparable  to  the  Doppler  width,  Awq.  This  is  because  for  this 
detuning,  the  minimum  velocity  (for  colinear  beams)  already  is 
approximately  the  average  thermal  velocity,  v,j.,  a  velocity  where  the 
Maxwellian  distribution  changes  rapidly.  Such  detunings  were  found  to  be 
necessary  experimentally  for  achieving  high  reflectivity.  The  thermal 
velocity  is  shown  as  the  circle  in  Figure  3-5. 

As  an  example  of  the  power  of  this  physical  model  for  achieving  a 
basic  understanding  of  FWM  physics,  one  can  conclude  from  examining  Figure 
3-5  that  the  interference  grating  formed  by  the  back  pump,  E^,  and  Eg 
("the  back  grating")  will  be  increasingly  unimportant  in  DFWM  experiments 
in  Doppler -broadened  systems  when  Aw  is  changed  from  zero.  The  very 
large  crossing  angle  requires  that  the  absorbers  which  can  record  the 
spatial  interference  pattern  of  the  back  grating  have  a  very  large 
velocity.  Therefore  only  a  small  fraction  of  the  total  population  can 
contribute  to  the  FWM  signal  via  the  back  grating  under  these  conditions. 

3.1.3  Other  Effects  due  to  Thermal  Motion 

In  addition  to  the  Doppler  effects  discussed  in  the  last  section, 
thermal  motion  of  the  absorbers  also  impacts  a  FWM  experiment  by  "washing 
out"  the  interference  gratings  formed  by  Ef+Eg  and  E^+Eg.  These 
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Figure  3-5.  Velocity  groups  in  which  the  forward  (E,/E  )  and  back  (E^/E  )  gratings 
are  written  for  DFWM.  The  circle  corresponds  to  the  thermal  velocity. 
The  back  grating  is  usually  unimportant  because  of  the  large  velocity 
required  to  be  simultaneously  resonant  with  E^  and  E^. 


gratings  contain  the  unique  spatial  variation  in  refractive  index  that 
leads  to  scattering  of  pump  energy  into  the  signal  and  conjugate 
directions  with  exactly  the  correct  phase  relationship.  The  grating 
spacing  for  two  plane  waves  Ef  and  crossing  at  angle  is 
given  by 


Af3  -  (A/2)/sin(tffg/2).  (3-8) 

For  arbitrary  waves,  the  spacing  between  grating  lines  can  vary  somewhat 
because  of  spatial  phase  variations.  The  grating  is  oriented  so  that 
grating  lines  are  parallel  to  the  bisector  of  i.e.  with  a  grating 

vector  defined  by  k£-kg.  It  is  the  velocity  component  along  the 
grating  vector  perpendicular  to  the  grating  lines  that  leads  to  thermal 
washout.  Interestingly,  the  Doppler  resonance  condition  selects  absorbers 
with  velocities  primarily  along  the  grating  lines,  i.e.,  with  velocities 
that  do  not  carry  absorbers  from  light  to  dark  regions.  It  is  only 
because  of  finite  homogeneous  width  that  there  is  any  cross  gracing  line 
velocity,  as  illustrated  in  Figure  3-6. 

The  maximum  cross  gracing  velocity  in  the  resonant  velocity  group  is 
just  v-(r/k)/sin(tf/2) ,  however  there  are  only  two  points  in  velocity 
space,  at  two  of  the  apices  of  the  population  diamond,  that  correspond  to 
resonant  absorbers  with  that  cross  grating  velocity.  The  average  grating 
destroying  velocity  is  the  weighted  average  over  the  upper  (lower)  half  of 
the  population  diamond 
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The  Integration  is  over  only  one  half  of  the  diamond  since  velocities  in 
either  direction  perpendicular  to  the  grating  are  equally  effective  in 
washing  out  the  grating  (<Vj.g>  is  zero  integrated  over  the  entire 
resonant  velocity  diamond,  but  the  gratings  do  indeed  wash  out!).  Note 
that  in  the  Doppler  limit,  T  ^  AWpjj  «  Awq,  the  velocity  weighting 
in  the  direction  of  the  grating  lines  (  v^ggonant  Figure  3-6)  is 
uniform  across  the  population  diamond  and  the  average  grating  destroying 
velocity  is  given  by 


<v 


((2RT/(irM))^/2 


(l.expl  *MrV{2RTk2sin2(«/2) )  ]  ^ 

erf (2RT)^/^sin(«/2) ) ] 


(3-10) 


Note  that  for  the  case  where  the  two  fields  are  colinear,  <Vg^>— (2RT/»rM)^'^^ , 
which  is  (2/ir)^/^  times  the  one -dimensional  rms  velocity. 

The  thermal  washout  time,  then,  is  the  time  required  for  absorbers  to 
travel  one-half  the  fringe  spacing,  i.e.. 


^washout  -  V[4sin(tf/2)<Vg^>]  (3-11) 

Several  aspects  of  grating  washout  due  to  thermal  motion  can  be 
appreciated  from  Equations  3-10  and  3-11.  First,  grating  washout  requires 
a  finite  homogeneous  width  r  for  because  it  is  this  width 

that  gives  rise  to  the  population  diamond  and  resonant  absorbers  with 
velocities  different  than  Aw/kcos(^/2) .  Secondly,  the  average  grating 
destroying  velocity  increases  with  F  if  one  is  still  in  the  limit  where 
the  mean  free  path  of  the  molecules  is  large  compared  to  the  physical 
dimensions  of  the  interaction  volume.  This  would  be  the  case,  for 
instance,  if  power-broadening  was  the  dominant  mechanism  for  homogeneous 
broadening,  in  which  case  F  should  be  replaced  by  tM^-^/2  in  the 
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above  equations . 


3.2  SODIUM  HYPERFINE  SPECTROSCOPY 

Many  of  the  experimental  results  discussed  in  Section  4  depend 
strongly  on  the  characteristics  of  the  particular  sodium  transition  used 
in  these  studies.  All  of  the  work  reported  here  was  performed  at  or  near 
the  resonance  of  the  ^^1/2*  ^3/2  transition,  commonly  called  the 
D2  line.  Because  Na  has  a  nuclear  spin  of  1-3/2,  the  purely  electronic 
transition  is  broken  into  six  closely  spaced  hyperfine  components. 
Therefore,  calculations  of  saturation  intensity,  transition  moment,  linear 
absorption,  etc.  need  to  account  for  the  hyperfine  components  that  are 
involved.  This  is  particularly  the  case  when  using  laser  sources  with 
linewidths  small  compared  to  the  hyperfine  splittings. 

Figure  3-7  presents  the  hyperfine  level  structure  of  the  D2  line 
and  a  stick  spectrum  showing  the  frequencies  and  relative  line  strengths 
of  the  six  dipole -allowed  hyperfine  transitions.  The  ground  state 
consists  of  two  hyperfine  levels,  F-1  and  F-2,  separated  by  1.77  GHz. 

This  splitting  in  the  ground  state  accounts  for  the  separation  of  the 
hyperfine  spectrum  into  two  separate  groups  of  three  closely  spaced 
transitions,  the  F-1  to  F'-0,l,2  group  and  another  corresponding  to  the 
F-2  to  F'-l,2,3  lines.  This  second  group  is  stronger  in  absorption,  both 
because  the  transition  moments  are  larger  and  because  the  F— 2  state  has 
5/8  of  the  ground  state  population  due  to  the  relative  degeneracies  of  the 
two  ground  state  hyperfine  levels.  Relative  line  intensities,  effective 
Einstein  A  coefficents,  and  dipole  transition  moments  of  the  six 
transitions  are  listed  in  Figure  3-8.  Because  of  the  magnetic  degeneracy 
of  the  F  and  F'  hyperfine  levels,  one  must  actually  consider  the  dipole 
moments  of  transitions  between  magnetic  sublevels  within  a  particular 
hyperfine  component  to  obtain  the  correct  transition  moment.  This  results 
from  the  dependence  of  magnetic  sublevel  selection  rules  on  the 
polarization  of  the  field,  i.e.  AM-0  for  linear  polarization,  +1  for 
right  circular  polarization,  and  -1  for  left  circular  polarization.  The 
transition  moments  for  Individual  magnetic  sublevel  transitions  in  the  F-2 
to  F'-3  and  F-1  to  F'-O  transitions  are  listed  in  Figure  3-8  as  examples 
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ATOMIC  SPECIES  -  Na  -  5890  Angstrom  Line 

UPPER  STATE  ELECTRONIC  ANG.  MOMENTUM  -  3/2 
LOWER  STATE  ELECTRONIC  ANG.  MOMENTUM  -  1/2 
NUCLEAR  MOMENTUM  -  3/2  -  , 

ELECTRONIC  TRANSITION  A  COEFFICIENT  -  6.22  x  10  s” 

TOTAL  ELECTRONIC  TRANSITION  DIPOLE  MOMENT  -  1.27  xl0“  esu-cm 
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Figure  3-8.  Spectroscopic  parameters  of  the  Na  line. 
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since  these  two  transitions  are  the  most  important  in  the  FWM  work 
reported  here. 

The  frequency  separations  of  the  h3rperfine  transitions  are  large 
compared  to  the  bandwidth  of  the  CW  ring  dye  laser  used  in  our  experiments 
(<  1  MHz),  however  the  effects  of  Doppler  broadening  make  it  impossible  to 
isolate  individual  hyperfine  components.  Figure  3-9  shows  the  overlapping 
Doppler-broadened  lineshapes  of  the  six  transitions  (T  -  500  K) ,  along 
with  their  sum,  which  is  the  absorption  lineshape  expected  for  the  D2 
line.  It  is  readily  apparent  that  a  laser  at  frequency  can  be 
simultaneously  resonant  with  several  different  hyperfine  transitions, 
albeit  with  different  velocity  groups  because  Aw  will  be  different  for 
the  different  transitions.  Each  transition  will  contribute  to  the  FWM 
response  at  w^  according  to  its  dipole  moment  and  Aw,  which 
determines  the  population  fraction  and  the  amount  off- resonance  the 
nonlinearity  is  "read"  in  a  DFWM  experiment.  This  concept  is  illustrated 
in  Figure  3-lOa,  which  shows  the  two-dimensional  population  diamonds 
resonant  with  the  different  hyperfine  transitions  when  the  laser  is  tuned 
to  the  1-0  frequency,  Clearly,  there  is  substantial  population 

with  velocities  such  that  the  laser  is  also  resonant  with  the  1-1  and  1-2 
transitions  because  their  (angular)  frequency  offsets  from  W]^_q  (100 
and  314  MHz  respectively)  are  small  compared  to  the  Doppler  width  of  10.7 
GHz  at  500  K.  The  fact  that  the  F— 2  to  F'— 1,2,3  transitions  are  about 
10.7  GHz  away  suggests  that  they  also  contribute  very  weakly  due  to 
velocity  components  in  the  wings  of  those  lines. 

The  contribution  of  the  1-1  and  1-2  transitions  to  the  FWM  response 

at  Wl“‘^1-0  limited  however  by  the  effects  of  optical  pumping. 

Optical  pumping  arises  when  population  pumped  from  the  lower  hyperfine 
state  (F-1)  into  a  particular  upper  hyperfine  state  can  radiate  to  other 

lower  state  hyperfine  levels  (in  this  case  F-2) .  Since  the  F-2  state  is 

not  connected  with  any  of  the  upper  hyperfine  levels  by  the  optical  field, 
the  effect  of  driving  the  1-1  and  1-2  transitions  is  to  pump  population 
from  one  of  the  groundstate  hyperfine  components  (F-1)  to  the  other 
(F"2) .  The  F— 1  population  is  removed  in  those  velocity  groups  where  the 
1-1  and  1-2  transitions  are  resonant  with  the  laser,  and  hence  the 
contribution  of  these  transitions  to  the  FWM  response  is  reduced  to  zero. 

Actually,  the  F-1  population  is  not  completely  removed  due  to  the 


mitigating  impact  of  diffusion.  Na  atoms  with  the  appropriate  velocities 
diffuse  from  unradiated  regions  into  the  laser  beam;  the  steady  state  F-1 
population  is  determined  by  the  optical  pumping  rate  compared  to  the 
diffusional  rate.  As  a  result,  the  1-1  and  1-2  transitions  will  have  a 
nonzero  contribution  to  the  FWM  response.  For  low  optical  power  DFWM 
experiments ,  the  contribution  will  be  very  small  because  the  nonlinearity 
is  being  "read"  off  resonance  as  discussed  above  in  Section  3.1.  Another 
mechanism  for  mitigating  the  effects  of  optical  pumping  is  to  consider  the 
overlap  of  the  wings  of  the  F-2  to  F'-l,2  transitions  with  the  laser 
frequency.  This  provides  a  mechanism  for  optical  pumping  from  F-2  to  F-1, 
counteracting  the  optical  pumping  in  the  other  direction.  However, this 
mechanism  can  only  occur  when  the  homogeneous  width  is  substantial  because 
the  two  velocity  groups  involved  in  the  F-1  transitions  and  the  F-2 
transitions  are  very  different.  Velocity  changing  collisions  in  fact 
would  lead  to  the  large  homogeneous  width  required,  however  such 
conditions  were  not  realized  in  the  experiments  reported  here. 

When  strong  optical  fields  are  employed,  the  situation  is  somewhat 
more  complex  than  illustrated  in  Figure  3-lOa,  Under  these  conditions, 
the  homogeneous  width  increases  due  to  power  broadening  according  to 

A«pb  -  2r(l+I/l3^^)V2  (3.;L2) 

where  F  is  the  homogeneous  half  width  at  low  power  and 

^sat  ”  (3-13) 

is  the  local  saturation  Intensity  for  velocity  group  v  and  laser  detuning 
Aw.  Figure  3-lOb  illustrates  the  case  when  the  power -broadened 
homogeneous  width  exceeds  the  frequency  separation  of  the  1-0, 1-1,  and  1-2 
transitions  when  Wj^  is  tuned  to  While  the  population 

diamond  for  the  1-0  transition  is  substantially  larger  because  of  the 
increased  homogeneous  width,  much  of  the  population  is  subject  to  optical 
pumping  because  of  the  overlapping  diamonds  of  the  1-1  and  1-2 
transitions.  As  discussed  above,  the  F-1  population  for  velocities  in  the 
overlapping  regions  is  reduced. 

In  addition  to  the  overlapping  populations  in  velocity  space,  one 


must  consider  the  effect  of  strong  optical  fields  on  the  nonlinearity. 

For  strongly  saturating  fields,  the  maximum  nonlinear  response  In 
frequency  space  can  occur  many  homogeneous  llnewldths  off  resonance  (see 
discussion  In  Section  3.1).  As  a  result,  the  on-resonance  population 
diamonds  Illustrated  In  Figure  3 -10a  and  10b  may  be  contributing  only 
weakly  to  the  FWM  response.  Off -resonance  population  diamonds,  where  the 
optical  field  can  sample  x  the  frequency  where  the  contrast  Is  larg¬ 
est,  can  contribute  more  If  they  have  a  sufficient  population  fraction. 

The  three  hyperflne  components  still  remain  coupled  In  the  same  way  as 
Illustrated  In  Figure  3-lOb.  but  It  Is  the  overlap  of  off -resonance 
diamonds  that  Is  pertinent.  Because  the  power  broadening  Is  dependent  on 
Aw-k*v,  the  populations  of  diamonds  decrease  with  Increasing  detun¬ 
ing  from  the  laser;  this  Is  Illustrated  In  Figure  3-11.  Actually,  It  Is 
only  the  area  (volume)  of  the  population  diamonds  that  decreases;  the 
velocity  distribution  must  also  be  considered  to  estimate  the  actual 
change  In  participating  population. 

The  coupling  of  the  1-0,  1-1,  and  1-2  transitions  Is  just  the  overlap 
of  three  such  grids,  each  centered  on  Its  own  resonant  velocity  group.  It 
should  be  noted  that  the  spaclngs  In  the  three  grids  are  not  Identical 
because  Ig^^(Aw-0 ,v-0)  Is  different  for  each  of  the  transitions. 
Furthermore ,  It  should  be  appreciated  that  the  population  diamonds  are  not 
hard  "boxes";  rather  the  grid  lines  In  Figure  3-11  represent  the 
Doppler-shifted  velocities  corresponding  to  the  half-height  of  Lorentzlan 
distributions  In  frequency  space.  Thus  there  Is  some  overlap  of  tails  in 
the  population  diamond  grid  that  is  not  adequately  represented  in  diagrams 
of  population  diamonds . 

3.3  ANALYTICAL  APPROACH 


A  major  objective  of  this  program  was  to  determine  FWM  character- 
isitics  and  capabilities  in  sodium  vapor  under  the  high  CW  reflectivity 
conditions  that  many  missions  are  likely  to  require.  As  a  result,  the 
experiments  were  neither  in  the  homogeneous  or  Doppler  limits  even  though 
the  laser  bandwidth  was  generally  very  narrow  in  comparison  with  F,  the 
intrinsic  homogeneous  width  of  the  system.  This  is  because  of  substantial 
power-broadening  of  the  lines  (c.f.  Equation  3-12)  Induced  by  the  strong 
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optical  pumping  fields.  For  the  50-100  W/cm  pump  fields  employed  In 
our  experiments,  I/Ig^t  10^-10^  giving  rise  to  power -broadened 
llnewldths,  30-100  times  the  collision- free  Intrinsic 

homogeneous  width  of  r-62.5  MHz.  Thus  the  homogeneous  width  was  20-60% 
of  the  Doppler  width,  too  large  to  satisfy  the  Doppler  limit  condition 
that  the  homogeneous  llnewldth  (Including  power  broadening)  be  very  small 
compared  to  the  Doppler  width.  The  Doppler  limit  Is  tantamount  to  saying 
the  field  Interacts  only  with  absorbers  within  a  very  narrow  velocity 
spread,  l.e.,  the  effects  of  thermal  motion  In  the  ensemble  can  be 
Ignored.  The  Doppler -broadened  regime  Is  Inherently  an  "on-resonance" 
small-signal  regime  so  the  nonlinearity  Is  primarily  an  absorptive  ix") 
one.  For  the  experimental  conditions  In  our  experiments,  so  few 
homogeneous  subcomponents  span  the  line  -  .2Awq)  the 

assumption  that  x'  doesn't  saturate  Is  no  longer  valid.  In  fact, 

Woerdman  and  Schuurmans  (Reference  3-3)  present  results  for  experiments  In 
the  same  experimental  regime  Indicating  that  to  first  order,  the  nonlinear 
polarization  obtained  In  the  absence  of  strong  linear  absorption  or  ptunp 
depletion  can  be  modeled  reasonably  well  as  a  purely  dispersive  suscept¬ 
ibility  with  saturation  effects. 

Our  experimental  conditions  for  high  reflectivity  did  not  satisfy  the 
homogeneous  limit  either.  Because  ^p^j  <  total  ensemble 

of  sodium  can  behave  homogeneously  only  when  the  laser  detuning  Aw  is 
much  larger  than  Awjj,  a  situation  where  the  complicating  effects  of 
thermal  motion  can  be  Ignored  (except  for  grating  washout) .  All 
velocities  respond  to  the  field  essentially  identically  because  Aw  » 
k»Vj^  for  all  1.  In  this  case,  the  nonlinearity  is  strictly 
dispersive  because  the  Interaction  is  far  off- resonance .  The 

Intermediate  case,  where  F  <  ^pjj  <  Awjj  is  substantially  more 
difficult  to  treat  because  many  velocity  groups  contribute  to  varying 
degrees,  just  as  one  would  expect  from  the  physical  picture  described 
above.  TRW  started  an  effort  to  describe  this  experimental  region 
theoretically;  the  effort  is  described  in  Section  7.  However,  to  date 
there  is  no  tractable  analytic  solution  for  this  regime.  For  this  reason, 
the  experimental  results  will  be  discussed  primarily  in  terms  of  the 
physical  picture.  The  Importance  of  FWM  in  homogeneously  broadened 
systems  was  not  discounted,  however,  and  several  models  were  developed 


optical  pumping  fields.  For  the  50-100  V/cm  pump  fields  employed  In 
our  experiments,  I/lg^t  is  10^-10^  giving  rise  to  power -broadened 
llnewldths,  dWpjj,  30-100  times  the  collision-free  Intrinsic 
homogeneous  width  of  r-62.5  MHz.  Thus  the  homogeneous  width  was  20-60% 
of  the  Doppler  width,  too  large  to  satisfy  the  Doppler  limit  condition 
that  the  homogeneous  llnewldth  (Including  power  broadening)  be  very  small 
compared  to  the  Doppler  width.  The  Doppler  limit  Is  tantamount  to  saying 
the  field  Interacts  only  with  absorbers  within  a  very  narrow  velocity 
spread,  l.e.,  the  effects  of  thermal  motion  In  the  ensemble  can  be 
Ignored.  The  Doppler -broadened  regime  is  inherently  an  "on-resonance" 
small-signal  regime  so  the  nonlinearity  is  primarily  an  absorptive  (x") 
one.  For  the  experimental  conditions  in  our  experiments,  so  few 
homogeneous  subcomponents  span  the  line  --  .  2^0)^)  the 

assumption  that  x'  doesn't  saturate  Is  no  longer  valid.  In  fact, 

Woerdman  and  Schuurmans  (Reference  3-3)  present  results  for  experiments  in 
the  same  experimental  regime  Indicating  that  to  first  order,  the  nonlinear 
polarization  obtained  In  the  absence  of  strong  linear  absorption  or  pump 
depletion  can  be  modeled  reasonably  well  as  a  purely  dispersive  suscept¬ 
ibility  with  saturation  effects. 

Our  experimental  conditions  for  high  reflectivity  did  not  satisfy  the 
homogeneous  limit  either.  Because  AUp^  <  Awq,  the  total  ensemble 
of  sodium  can  behave  homogeneously  only  when  the  laser  detuning  Is 
much  larger  than  a  situation  where  the  complicating  effects  of 

thermal  motion  can  be  Ignored  (except  for  grating  washout)  .  A.11 

velocities  respond  to  the  field  essentially  identically  because  Aw  » 
k*Vj^  for  all  1.  In  this  case,  the  nonlinearity  is  strictly 
dispersive  because  the  Interaction  Is  far  off- resonance .  The 
intermediate  case,  where  F  <  AWpj^  <  Awp  Is  substantially  more 
difficult  to  treat  because  many  velocity  groups  contribute  to  varying 
degrees,  just  as  one  would  expect  from  the  physical  picture  described 
above.  TRW  started  an  effort  to  describe  this  experimental  region 
theoretically;  the  effort  is  described  in  Section  7.  However,  to  date 
there  is  no  tractable  analytic  solution  for  this  regime.  For  this  reason, 
the  experimental  results  will  be  discussed  primarily  in  terms  of  the 
physical  picture.  The  importance  of  FWM  in  homogeneously  broadened 
systems  was  not  discounted,  however,  and  several  models  were  developed 
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during  this  program.  They  are  also  discussed  In  Section  7. 
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4.0  REFLECTIVITY  MEASUREMENTS 


This  section  reports  experimental  data  and  analysis  of  the 
observed  functional  dependence  of  FWM  reflectivity  on  parameters  of 
practical  concern  such  as  pump  intensity,  frequency  offset,  and  signal 
viewing  angle  (field  of  view).  Data  was  obtained  for  both  narrowband 
(Ai/  ~  1  MHz)  and  broadband  (Av  *•  2  GHz)  laser  operation,  and 
interpreted  in  terms  of  the  model  presented  in  Section  3. 

4.1  REFLECTIVITY  MEASUREMENT  OBJECTIVES 

Many  proposed  mission  applications  of  four  wave  mixing  technology 
require  phase  conjugation  of  optical  signals  returning  from  a  remote 
target.  Good  FWM  reflectivity  is  generally  necessary,  but  from  a 
practical  point-of-vlew  must  be  obtained  using  CW  or  long-pulsed  pump 
sources  because  of  target  range  uncertainty  that  translates  into  timing 
uncertainties  for  pump  pulses  in  the  FWM  medium.  For  this  reason,  the 
major  emphasis  of  our  reflectivity  experiments  was  to  determine  character¬ 
istics  of  sodium  FWM  using  modest  CW  pump  powers  in  a  configuration  that 
resulted  in  high  (>  100%)  reflectivity.  Frequency  offset"  of  peak 
reflectivity,  field-of -view,  and  intensity  dependence  were  measured 
because  all  have  significant  Impact  on  the  practical  implementation  of 
this  technology.  This  same  motivation  of  determining  FWM  reflectivity 
under  conditions  that  can  be  easily  implemented  also  gave  rise  to  our 
objective  to  understand  reflectivity  under  broadband  laser  excitation. 

Such  information  allows  system  designers  to  trade  off  FWM  performance  with 
complexity  and  fleldability  of  the  pump  lasers. 

4.2  EXPERIMENTAL  DESIGN 


All  degenerate  four  wave  mixing  experiments  shared  a  common 
skeletal  arrangement  which  is  diagrammed  in  Figure  4-1.  A  Coherent  20  W 
Ar'*'  pump  laser  and  699-21  ring  dye  laser  typically  generated  700-800  mW 
single  frequency  (measured  to  be  <1.5  MHz  FWHM)  or  1.3  W  broadband  (-2  GHz 
FWHM)  output  at  the  Na  D2  wavelength.  The  dye  laser  output  was 
s-polarized  to  all  reflecting  optics  in  both  the  pump  and  probe/noise  beam 


37 


paths  to  Insure  a  high  degree  of  copolarization  at  the  Na  cell .  Turning 
mirrors  (M)  were  visible  dielectric  or  MAXBRIte  coated  Zerodur 
substrate  optics  with  twentieth  wave  or  better  surface  figures.  A  50-50 
beamsplitter  (BSl)  served  to  generate  the  two  pvimp  beams  and  was  angle 
tuned  to  power  match  both  pump  legs  to  within  5%.  Lens  LI,  a  90.0  cm 
focal  length  spherical  achromat,  imparted  equal  convergence  to  each  pump 
leg  resulting  in  equal  intensity  counterpropagating  pump  fields  at  the  Na 
cell.  BS2  was  a  4%  first  surface.  AR  coated  second  surface  beamsplitter 
which  provided  a  sampling  of  the  main  beam  for  frequency  and/or  bandwidth 
analysis.  The  probe  and  noise  beams  originated  from  a  high  quality 
variable  beamsplitter  (VBS)  inserted  into  the  main  beam  path  which, 
depending  on  the  particular  experiment  at  hand,  would  allocate  the 
necessary  amount  of  power  into  the  probe/noise  leg. 

Laser  bandwidth  was  measured  with  either  a  150  MHz  free  spectral 
range  (FSR)  confocal  Fabry-Perot  etalon  with  a  finesse  of  100  for  single 
frequency  operation,  or  a  2  GHz  FSR  Fabry-Perot  of  finesse  50  for 
broadband  output.  Rough  laser  frequency  measurements  to  ~  1  GHz  accuracy 
were  made  by  a  Burleigh  wavemeter;  for  higher  accuracy  frequency 
measurements,  an  on  line  Lamb  dip  frequency  diagnostic  was  assembled. 
Following  the  experimental  approach  of  Hansch  ^  al.  (Reference  4-1),  a 
reference  beam  and  probe  beam  of  equal  intensities  were  produced  from  the 
two  surfaces  of  an  uncoated  flat  (see  Figure  4-1) .  The  two  beams  were 
directed  through  a  10  cm  Na  vapor  cell  which  typically  operated  at  400-410 
K.  A  third  saturating  beam  was  sent  through  the  Na  cell  in  a 
counterpropagating  direction  to  the  probe  and  reference  beams  and 
carefully  overlapped  the  probe  beam  path  in  the  cell.  The  output  of  the 
probe  and  reference  beams  was  detected  by  matched  photodiodes  and  sent  to 
a  high  sensitivity  differential  amplifier  whose  output  was  displayed  on  an 
oscilloscope.  With  proper  cell  temperature  and  beam  overlap,  two  lines 
were  clearly  detected  (lower  trace.  Figure  4-2)  representing  transitions 

O 

originating  from  the  hyperflne  split  3s  ground  state  of  the  Na 

D2  line.  Since  the  positions  of  these  transitions  are  accurately  known, 
frequency  accuracy  to  50  MHz  was  possible. 
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4.3  NARROWBAND  LASER  PUMPING  RESULTS  (A»/  <  1  MHz) 


4.3.1  DFWM  Reflectivity  Frequency  Spectrtim 

The  observed  frequency  dependence  of  the  DFWM  reflectivity  in  the 
narrowband  dye  laser  experiments  will  be  examined  in  this  section.  For 
small  pump  intensities,  it  has  been  shown  that  optical  pumping  effects 
play  an  important  role  in  the  FWM  of  the  sodium  D2  resonance  line 
(Reference  4-2).  When  resonantly  driving  the  D2  line,  only  the  2-3  and 
1-0  hyperfine  transitions  are  free  from  optical  pumping  and  are, 
therefore,  the  only  ones  which  contribute  significantly  to  the 
steady-state  DFWM  reflectivity  at  low  pump  powers.  This  was  discussed  in 
some  detail  in  Section  3.2.  At  high  pump  powers,  which  are  required  when 
large  DFWM  signals  are  desired,  all  six  D2  hyperfine  transitions  can 
contribute  to  the  signal,  but  in  a  complicated  way  that  depends  on 
optical  pumping,  diffusion  in  and  out  of  the  beam,  the  explicit  values  for 
the  detuning  from  each  of  the  hyperfine  transition  frequencies  (Aw^  - 
WL'«i)  I  *>^d  the  homogeneous  power-broadened  widths, 

Figure  4-3  shows  the  frequency  dependence  of  FWM  reflectivity  for 
experimental -conditions  giving  very  high  reflectivites  (R-  160-230%).  In 
this  case,  as  well  as  the  conditions  in  Figure  4-4  for  conditions  of  lower 
R,  the  maximum  reflectivity  of  the  2-3  transition  occurs  on  the  low 
frequency  side  of  its  Doppler-broadened  absorption  line  whereas  the 
maximum  reflectivity  for  the  1-0  transition  occurs  on  the  high  frequency 
side  of  its  Doppler -broadened  absorption. 

This  spectral  dependence  can  be  understood  as  follows.  In  all  of 
these  cases,  the  on-resonance  absorptive  component  is  well  saturated 
because  I/^sat  Because  the  dispersive  component  of  the 

nonlinearity  is  zero  at  resonance  (c.f.  Figure  3-1),  one  expects  a  dip  in 
FWM  reflectivity  at  resonance;  such  dips  have  been  reported  previously 
(References  4-3,  4-4).  In  our  experiments,  the  reflectivity  does  indeed 
drop  off  near  the  2-3  and  1-0  resonances  for  this  reason,  but  does  not 
recover  for  frequencies  between  those  two  resonances  because  of  the  very 
strong  linear  absorption  still  present  in  the  center  portion  of  the 
Doppler-broadened  D2  line.  The  absorption  is  strongly  saturated  in  this 
region,  being  10-100  times  smaller  than  the  small-signal  attenuation, 
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DFWM  REFLECTIVITY 


Figure  4-3.  Typical  frequency  spectrum  of  the  DFWM  reflectivity  under  high 
(R  >  100%)  conditions.  The  spectrum  was  obtained  by  scanning 
a  narrowband  dye  laser  through  the  Na  resonance.  Shown  for 
reference  is  the  relative  Doppler  absorption  profile  of  the  D„ 
line. 
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DFWM  Reflectivity  R  -  77.9% 


Figure  4-A.  Frequency  spectra  of  DFWM  reflectivity  under  lower 
R  conditions  a).  Under  low  R  conditions >  the  peak 
corresponding  to  the  2-3  transition  is  stronger,  as 
expected  from  state  degeneracy  and  oscillator  strength 
arguments  b).  At  higher  R,  the  peak  corresponding  to 
the  1-0  transition  becomes  larger.  On-resonance 
dips  in  reflectivity  (2-3  in  Figure  4-4a  and  1-0 
in  4-4b)  are  observed,  indicating  that  the  absorptive 
component  of  the  susceptibility,  x"»  is  still  contrib¬ 
uting  to  FWM  R  under  these  conditions. 
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Oq,  but  is  still  substantial  since  was  about  200  cm’^  in 
these  experiments.  At  lower  sodium  densities  (and  hence  lower  R)  than 
investigated  here,  the  linear  absorption  is  lower  and  FWM  reflectivity  can 
be  observed  at  these  intermediate  frequencies  (see  Reference  4-3).  Thus 
it  is  the  combination  of  higher  nonlinear  response  off- resonance  and  the 
strong  linear  absorption  necessarily  present  for  achieving  high 
reflectivity  that  accounts  for  the  observed  spectral  dependence. 

A  more  detailed  examination  of  Figure  4-3  shows  that  while  the 
reflectivity  goes  to  zero  at  the  2-3  resonance,  there  is  substantial 
reflectivity  at  the  1-0  resonance  frequency.  Two  factors  account  for 
this.  First,  the  linear  absorption  is  highest  at  the  2-3  frequency, 
remaining  higher  than  that  at  the  1-0  frequency  throughout  the  wing  of  the 
2-3  transition.  A  second  factor  is  the  relative  contributions  from  the 
other  hyperfine  transitions.  Consider  the  FWM  response  at 
absorptive  component  of  the  1-0  transition  is  well  saturated  and  its 
dispersive  component  is  zero.  However,  at  this  frequency,  the  dispersive 
components  of  the  1-1  and  1-2  transitions  are  not  zero.  The  transition 
moments  of  the  1-1  and  1-2  transitions  are  both  stronger  than  the  1-0  and, 
as  discussed  in  Section  3.2,  they  can  contribute  signal  if  diffusion  into 
the  optical  field  is  sufficiently  fast  to  offset  the  effects  of  optical 
pumping.  A  similar  argument  holds  for  response  at  W23  the  2-2  and 
2-1  transitions  are  weak  in  comparison  to  the  2-3  transition  and  their 
contribution  is  further  degraded  by  the  stronger  linear  absorption. 

The  role  of  linear  absorption  in  determining  the  spectral  response  is 
again  evident  in  Figure  4-5,  which  presents  the  observed  intensity 
dependence  of  the  FWM  reflectivity  spectrum.  As  the  pump  power  increases, 
the  hole  between  the  two  reflectivity  peaks  begins  to  narrow  because  of 
increasing  saturation  of  the  linear  absorption.  The  reflectivity  peaks 
broaden  mostly  towards  the  center  of  the  £>2  line  because  that  is  where 
the  absorption  has  the  strongest  effect.  Another  factor  leading  to  the 
increased  FWM  bandwidth  is  that  the  homogeneous  width  increases  with 
intensity.  This  enlarges  the  velocity  population  diamond,  increasing  the 
effective  number  den.sity  interacting  in  the  experiment  for  a  given  laser 
frequency.  This  mechanism  also  predicts  that  the  broadening  is  primarily 
towards  linecenter  of  the  D2  line.  Because  of  the  Maxwellian  velocity 
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Figure  4-5.  Observed  dependence  of  DFWM  frequency  spectrum  as  a  function 
of  pump  power.  The  spectra  were  obtained  by  scanning  a 
narrowband  dye  laser  through  the  Na  line  resonance. 
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distribution,  the  number  density  increases  more  rapidly  on  the  side  of  the 
population  diamond  that  moves  towards  lower  velocity  (and  hence  line 
center)  than  on  the  side  moving  away  towards  higher  Doppler  velocity. 

Another  interesting  feature  of  the  intensity  depenaence  of  the  FWM 
reflectivity  spectrum  can  be  seen  by  comparing  the  spectra  in  Figures  A- 3 
and  4-4.  There  is  dramatic  change  in  relative  amplitude  of  the  two 
reflectivity  peaks  as  the  pvunp  Intensity  is  increased.  Despite  the 
complicated  dependence  of  the  DFWM  signal  on  the  parameters  mentioned  at 
the  beginning  of  this  section,  transitions  originating  from  the  F-2  lower 
hyperfine  level  are  expected  to  contribute  more  to  the  DFWM  signal  for  two 
reasons.  First,  they  carry  more  oscillator  strength  than  those 
originating  from  the  F-1  lower  hyperfine  level  as  can  be  seen  from  the 
values  presented  in  Figure  3-8.  Note  that  the  dominant  contributing 
transition  from  F-2  is  the  2-3  transition,  which  has  the  strongest 
transition  moment  of  all,  while  the  domlnent  transition  from  F-1  is  the 
1-0,  which  has  the  second  smallest  transition  moment.  The  saturation 
intensity  is  proportional  to  l//x^;  the  nonlinear  term  is 

proportional  to  Secondly,  the  F-2  hyperfine  level  has  5/8  of  the 

total  ground  state  population,  so  there  are  more  absorbers  contributing  to 
the  signal  due  to  transitions  from  this  level. 

Hence  the  DFWM  reflectivity  would  be  expected  to  be  largest  on  the 
low  frequency  side  of  the  broad  absorption  line  shown  in  Figure  3-9.  This 
is  in  fact  what  is  observed  under  the  low  pump  power,  low  FWM  reflectivity 
conditions  (R  »  1.6%)  shown  in  Figure  4-4a.  At  even  lower  powers,  the 
reflectivity  due  to  the  F-2  transitions  is  much  larger  than  that  due  to 
the  F-1  transitions  (References  4-2, 4-5). 

Under  high  pump  power  conditions,  the  2-3  peak  is  observed  to  be 
substantially  smaller  than  the  1-0  peak  on  the  high  frequency  side  of  the 
D2  line.  This  is  Illustrated  in  Figure  4-4b  and  4-3  for  cases  where  the 
FWM  reflectivity  is  approximately  75%  and  160%.  Such  a  changeover  in  the 
relative  magnitude  of  reflectivity  arising  from  the  F-2  and  F-1 
transitions  has  been  observed  in  pulsed  experiments  (Reference  4-6).  Peak 
reflectivity  for  high  pump  power  CW  conditions  also  has  been  reported  to 
be  near  the  1-0  transition  rather  than  the  2-3  transition  (Reference  4-7). 

This  apparent  anomaly  of  the  peak  reflectivity  shifting  from  the 
low  frequency  side  of  the  line,  where  the  contribution  to  FWM  reflectivity 


is  expected  to  be  largest,  to  the  high  frequency  side  as  the  pump  power  is 
increased  arises  from  another  nonlinear  effect.  This  nonlinear 
interaction  involves  changes  in  the  propagation  characteristics  of  the 
intense  pump  fields  induced  by  spatial  intensity  variations  in  the  same 
fields.  Two  of  these  self>action  effects,  self -focusing  and 
self-defocusing,  are  responsible  for  the  observed  Intensity  effects  in  the 
FWM  reflectivity  spectrum. 

In  nonlinear  media,  intensity  dependent  changes  in  the  material's 
refractive  index  cause  self -focusing  or  self-defocusing  of  propagating 
laser  beams  that  have  nonuniform  intensity  profiles.  The  equation 
describing  the  nonlinear  index  change  (An)  for  Gaussian  laser  beams 
propagating  in  resonant  media  is: 

-xNIpI^E^ 

An  -  - 2--  (4-1) 

[/I(Wo-Wl)]^ 

where  n  is  the  transition  dipole  moment,  N  is  the  (atomic)  number 
density,  is  the  laser  frequency  and  is  the  transition 
frequency  and  is  the  electric  field  amplitude.  It  can  be  seen  from 
this  equation  that  when  the  laser  is  timed  to  a  frequency  above  the 
transition  frequency,  An  >  0.  This  corresponds  to  a  positive  lens  and 
the  nonlinear  medium  will  have  a  tendency  to  focus  the  laser  beam.  For 
laser  frequencies  below  the  transition  frequency  An  <  0  and  the 
nonlinear  medium  will  act  as  a  negative  lens,  tending  to  defocus  the 
propagating  laser  beams.  Furthermore,  it  has  been  shown  (Reference  4-8) 
that  the  minimum  power  for  self- focusing  to  affect  a  Gaussian  laser  beam 
is  given  by: 


3.174  X  10'^  A^c 


where  A^  is  the  vacuvun  wavelength  of  the  laser,  c  is  the  speed  of 
light,  and  n2  ••  2Ln/\E^\^ .  is  on  the  order  of  100  /iW  for  the 

conditions  in  the  experiments  presented  here. 

The  DFWM  pump  powers  employed  in  TRW's  experiments  were  in  the  range  of  100-600 
mW,  so  self- focusing  (self-defocusing)  is  expected  to  impact  the  observed  spectral 
behavior.  The  1-0  transition  has  a  larger 
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reflectivity  than  expected  which  grows  relative  to  the  2-3  reflectivity 
peak  as  I  increases  because  self- focusing  increases  the  effective  pump 
intensity  and  hence  DFWM  reflectivity  at  this  frequency.  Conversely,  the 
2-3  transition  has  a  smaller  reflectivity  than  expected  because  the 
self-defocusing  on  the  low  frequency  side  of  2-3  transition  decreases  the 
pump  intensity.  Note  that  self-defocussing  is  inherently  a  self-limiting 
process  and  so  is  not  expected  to  play  a  significant  role  in  determining 
the  FWM  reflectivity  near  the  2-3  transition.  Further  evidence  in  support 
of  this  explanation  comes  from  both  previous  observations  of 
self-focusing/self-defocusing  in  sodivim  vapor  (Reference  4-9)  and  from 
observations  of  greatly  increased  beam  divergences  in  our  experiments. 

This  is  discussed  in  more  detail  in  the  next  section. 

4.3.2  Self-action  Effects 

The  self-action  effect  of  self-focusing  and  self-defocusing  of  a 
single  narrowband  CW  dye  laser  beam  in  sodium  vapor  has  been  examined 
previously  in  detail  (Reference  4-9),  The  purpose  of  the  experiments 
described  here  was  to  characterize  the  effects  of  self-focusiftg  and 
self-defocusing  on  the  DFWM  process.  These  often  strong  effects  include 
large  changes  in  reflectivity  and  in  the  divergences  of  all  laser  beams , 
including  the  conjugate  signal  beam.  One  consequence  of  the  latter  is 
that  the  phase  conjugate  fidelity  can  be  degraded  at  high  FWM  reflectivity 
because  of  the  large  pump  powers.  This  may  limit  the  conditions  under 
which  an  atomic  vapor  phase  conjugator  may  be  usefully  operated. 

Both  transmitted  pump  beams  and  the  conjugate  return  were  observed  to 
have  substantially  Increased  divergence  in  high  reflectivity  experiments 
where  the  sodium  density  was  high  and  pump  intensities  were  large.  An 
experiment  to  understand  and  quantify  this  effect  was  performed  using  the 
standard  counterpropagating  pump  geometry  shown  in  Figure  4-6.  In  this 
particular  experiment,  care  was  taken  to  ensure  that  the  focal  points  of 
the  two  counterpropagating  pump  beams  coincided  so  that  and  E^j  were 
phase  conjugates  at  all  points  along  the  pump  axis.  Furthermore,  the  Na 
cell  was  located  at  this  common  focal  point.  The  cell  length  of  1  cm  was 
short  compared  to  the  calculated  Rayleigh  distance  of  7.7  cm,  so  the  pump 
wavefronts  had  a  minimum  of  curvature  in  the 
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Figure  4-6.  Experimental  arrangement  for  measuring  self-action  effects  under  high 
R  conditions  in  Na  vapor.  Spot  size  was  measured  as  a  function  of 
distance  using  a  translating  iris  diaphragm. 


Interaction  region.  The  signal  beao,  incident  at  angle  0^^  with 
respect  to  the  forward  pump,  was  also  focused  at  the  center  of  the  cell. 

The  conjugate  return  beam  diameter  was  measured  as  a  function  of 
distance  from  the  Na  cell  and  found  to  be  diverging  for  experiments  on 
both  the  2-3  and  1-0  sides  of  the  line.  Figure  4-7  presents  these  data 
for  three  different  Na  densities.  Figures  4- 7a  and  4- 7b  show  results  for 
Na  densities  where  self-action  effects  were  observed,  while  4-7c  shows 
observed  divergence  in  the  absence  of  appreciable  self -focusing. 

Geometric  constraints  prevented  measurements  closer  than  about  60  cm  from 
the  cell.  However,  estimated  conjugate  beam  behavior  near  the  Na  cell  was 
obtained  by  fitting  the  data  for  Gaussian  beam  propagation.  Two  different 
Gaussian  beam  profiles  are  shown  in  Figure  4-7a  for  the  1-0  transition 
data:  one  that  matched  beam  dimensions  at  0  and  100  cm  using  the  measured 
probe  spot  with  the  cell  removed,  and  the  other  that  matched  the  observed 
divergence  in  the  60-100  cm  region.  Both  profiles  indicate  that  the  1-0 
conjugate  beam  is  converging  as  it  leaves  the  Na  cell,  then  diverging. 
Figures  4- 7b  and  c,  in  which  the  lines  are  Gaussian  beams  matched  to  the 
observed  divergences  and  propagated  back  to  the  effective  focal  points , 
show  that  the  effective  convergence  imparted  by  the  Na  cell  increased  as 
the  sodium  density  was  increased.  This  observation,  along  with  the  fact 
that  the  2-3  divergence  is  essentially  the  same  for  all  three  conditions, 
suggests  that  the  effect  is  due  to  self-focusing. 

However,  neither  of  the  Gaussian  beam  profiles  accurately  matches 
the  conjugate  beam  propagation  characteristics.  It  is  veiry  unlikely  that 
the  beam  diameter  of  the  probe  in  the  cell  is  that  measured  in  the  absence 
of  sodium  vapor  because  of  self-focusing  induced  by  the  strong  pumps; 
furthermore  the  observed  beam  divergence  is  not  predicted  well.  The 
divergence -matched  profile  requires  inordinately  large  beam  diameters  in 
the  cell,  inconsistent  with  self-focussing  and  filamenting  of  the  pump 
beams.  This  strongly  suggests  that  the  1-0  conjugate  return  is  not  a 
Gaussian  beam,  even  though  the  three  input  waves,  E^,  E^,  and  Eg, 
are  Gaussian.  Such  behavior  was  predicted  by  Grynberg,  et  al.  (Reference 
4-10)  who  showed  that  intensity  overlap  of  comparably-sized  Gaussian  input 
beams  would  yield  a  non-Gaussian  conjugate  wave  because  of  the  intensity 
dependence  of  the  reflectivity. 

Such  an  effect  should  also  impact  the  2-3  conjugate  return.  It  can 
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be  seen  from  Figure  4-7c  that  the  divergence  of  the  2-3  conjugate  beam  is 
significantly  greater  than  that  of  a  freely  propagating  Gaussian  input 
signal  beam  (dotted  line) .  In  these  experiments  the  size  of  the  pump 
beams  was  -  200  pm  (radius  of  1/e  E- field)  whereas  that  of  the 
probe  was  w^  -  150  pm.  For  these  conditions  the  average  pump 
intensity  experienced  by  the  probe  beam  varies  greatly  across  the  beam. 

For  example,  the  central  10%  of  the  probe  beam  will  interact  with  a  pvunp 
intensity  which  is  -6  times  greater  than  the  pump  intensity  that  the  outer 
10%  (r  -  .9  Wq  to  r  -  w^)  of  the  probe  beam  will  interact  with. 

Depending  upon  the  exact  experimental  conditions,  the  phase  conjugate 
reflectivity  can  change  appreciably  across  this  aperture  due  to  variations 
in  intensity.  For  example,  for  strongly  saturating  conditions,  the  outer 
edge  may  be  reflected  more  efficiently.  Thus  when  the  Gaussian  pump  and 
probe  beams  are  comparable  in  size,  the  probe  beam  and  the  phase  conjugate 
beam  will  be  diffracted  in  the  DFWM  process  by  a  gain  aperture  which  is  a 
complicated  function  of  the  pump  transverse  intensity  profile. 

This  effect  can  alter  the  effective  conjugate  spot  size  in  the 
interaction  region,  leading  to  return  beams  that  are  not  true  conjugates 
of  the  input  signal  wave.  For  the  lower  sodium  number  density  condition 
of  Figure  4- 7c,  the  aperturing  effect  appears  to  be  sufficiently  dominant 
to  mask  the  differences  in  the  divergences  of  the  self- focused  and 
self -defocused  phase  conjugate  beams.  At  higher  sodium  number  densities, 
the  self-action  nonlinearity  is  larger  so  differences  in  self-focusing  and 
self-defocuslng  become  evident. 

Even  stronger  lenslng  effects  than  those  shown  in  Figure  4-7  have 
been  observed  in  our  experiments  at  higher  pump  intensities  and  sodium 
number  density.  These  focusing  effects  have  limited  our  maximum 
obtainable  DFWM  reflectivity  in  sodium  vapor.  Despite  this,  we  have 
measured  CW  DFWM  reflectivity  of  up  to  230%,  which  is  about  a  factor  of 
two  larger  than  previously  reported  values  in  the  literature. 

Importantly,  these  lensing  effects  were  not  observed  to  be  significant  for 
reflectivities  below  40-50%.  Good  resolution  fidelity  and  aberration 
correction  was  observed  in  experiments  in  this  lower  reflectivity  regime. 
These  results  are  reported  in  Section  6. 


4.3.3  Saturation  Effects 


The  gain  aperturlng  effect  mentioned  above  depends  strongly  on  the 
saturating  effects  of  the  overlapping  Input  beams.  Saturating  effects 
were  discussed  In  general  In  Section  3;  this  section  will  examine  FWM 
saturation  effects  In  sodium  vapor  In  the  Rabl  regime  where  2r  < 

<  Our  high  reflectivity  measurements  were  In  this  experimental 

regime  where  saturation  characterlsltlcs  were  neither  those  of  a  homo¬ 
geneous  system  (F  or  »  Aq)  nor  those  of  a  system  In  the 
Doppler  limit  (AWp^  «  Awjj) . 

In  homogeneously -broadened  absorption  lines  the  saturation 
Intensity  Increases  by  a  factor  of  (1  +  5  )  as  the  laser  Is  detuned  by 
S  -  Au/r  homogeneous  half -widths  from  line  center.  In  contrast,  the 
saturation  Intensity  In  the  Doppler  limit  Is  Independent  of  the  detuning 
of  the  laser  from  line  center  because  the  laser  Is  on  resonance  with  the 
one  velocity  group  with  which  It  can  Interact.  Thus  5«0  regardless  of 
the  laser  detxinlng  from  line  center  where  v»0  because  S  Is  measured  from 
the  Doppler-shifted  line  center  frequency  of  the  resonant  velocity  group, 
Aw-k«v,  l.e.,  S  -  (A<<^k*v)/r  (see  Equation  3-14). 

Theoretical  analysis  (Reference  4-11)  for  homogeneously-broadened 
lines  Indicates  that  maximum  reflectivity  as  a  function  of  pumping 
Intensity  (I)  occurs  I/lg^t  “  Tkls  condition  can  be  met  for 
arbitrarily  large  I  by  detuning  the  laser  to  Increase  The 

requisite  detuning  Is  approximately 

which  matches  the  condition  derived  In  Section  3.1  for  the  frequency 
offset  where  x  begins  to  saturate  differently  If  the  pumps  are 
significantly  stronger  than  the  signal  wave.  For  detunings  much  smaller 
than  this,  I/lgg^  »  1  and  the  nonlinearity  has  no  "contrast"  between 
light  and  "dark"  regions. 

Since  th^  absorption  line  In  the  Doppler  and  Rabl  regimes  Is  a 
(Gaussian)  distribution  of  homogeneous  absorption  lines,  each  with  a 
slightly  different  line  center  frequency.  It  Is  useful  to  consider  the 
relative  contribution  of  each  of  these  homogeneous  packets  to  the  total 
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DFWM  reflectivity.  Two  laser  parameters  determine  this;  the  laser 
frequency,  which  determines  the  Doppler  shift  k*v  and  the  effects  of 
linear  absorption,  and  the  laser  Intensity,  which  determines  x  as  a 
function  of  S.  The  contribution  of  a  particular  velocity  group  is  then 
determined  by  its  population  fraction  and  its  value  of  S .  Note  that  the 
value  of  S  depends  on  the  homogeneous  width  over  which  v  is  weighted. 

In  the  Doppler  limit,  v  is  essentially  constant  (-Aw/k)  over  IT]  as 
the  power -broadened  homogeneous  width  approaches  <v>*0  and 

S  becomes  Aw/AWpj^j,  the  value  in  the  homogeneous  limit. 

The  discussion  in  Section  3.1.2  showed  that  in  the  case  of  DFWM, 
the  back  pump  "reads"  the  nonlinearity  off-resonance  by  2k*v.  As  a 
result,  one  would  expect  the  highest  reflectivity  when  iF  is  equal  to 
this  offset,  i.e.,  the  peak  reflectivity  would  be  expected  at 

-  1/2  (I/Ig^^(Aw-0))V2r^  (4.4) 

Measured  values  of  Au  for  the  FWM  reflectivity  peaks  on  both  the  2-3  and 
1-0  transitions  were  3-4  times  larger  than  predicted  from  Equation  4-4. 
This  is  expected,  however,  because  Equation  4-4  does  not  take  into  account 
the  impact  of  residual  linear  absorption  or  the  effective  increase  of  I  on 
the  1-0  side  of  the  line  due  to  self -focusing.  Both  of  these  effects 
should  disappear  as  the  Na  density  is  reduced. 

4.3.4  Field  of  View 

Figure  4-8  presents  the  observed  dependence  of  DFWM  reflectivity  on 
^fs’  angle  between  the  forward  pump  and  signal,  under  conditions 
of  high  CW  reflectivity.  The  experimental  arrangement  for  making  these 
measurements  was  the  same  as  for  the  active  tracking  demonstration 
described  elsewhere.  The  scanning  geometry,  reproduced  in  Figure  4-8,  was 
chosen  to  minimize  changes  in  the  interaction  volume  as  the  angle  was 
swept.  The  experiment  was  also  designed  to  keep  any  residual  phase 
mismatch  due  to  misalignment  of  the  pump  beams  constant  as  the  signal  beam 
angle  was  scanned.  It  is  important  to  minimize  the  impact  of  these  other 
effects  that  can  reduce  FWM  reflectivity  so  that  the  true  angular 
dependence  under  high  R  conditions  can  be  determined. 

The  angular  dependence  of  the  reflectivity  can  be  understood 
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Figure  4-8.  Observed  angular  dependence  of  DFWM  reflectivity  under 
high  reflectivity  (R  >  150%)  conditions.  Shown  is  the 
predicted  sin“^  6  dependence. 


qualitatively  in  terms  of  the  angular  dependence  of  the  population  of 
sodium  atoms  comprising  the  refractive  index  grating.  When  the 

population  discs  (c.f.  Figure  3-4)  resonant  with  E£  and  Eg  coincide, 
yielding  the  largest  possible  sodium  atom  population  for  the  nonlinear 
interaction.  For  large  6^^  in  the  Doppler  limit,  the  velocity  is 
essentially  constant  over  the  population  diamond,  and  so  the  sodium 
density  available  for  forming  the  grating  is  just 

Nex=Ntot(®”/^'^''‘'^^®^PL'W/2RTk2cos2(e/2))](r2/k2sine)  (4-5) 

which  for  Aw— 0  varies  as  1/sin  B.  This  is  the  population  fraction 
that  interacts  with  the  optical  fields,  but  because  of  grating  washout  it 
is  not  the  effective  population  comprising  the  grating.  At  most  angles, 
the  grating  lifetime  is  determined  by  washout,  which  in  turn  is  determined 
by  the  grating  spacing  given  in  Equation  3-8.  Thus  the  effective  N 
for  large  8  varies  as  (sin  B)~  ,  and  the  expected  angular  dependence 
of  FWM  reflectivity  varies  as  (sin  8)~^  since  R  a  N^.  This 
simple  derivation,  based  on  the  physical  model  discussed  in  Section  3, 
matches  the  large  angle  predictions  from  more  complex  formulations 
published  by  Vandzura  (Reference  4-12)  and  Nilsen  and  Yariv  (Reference 
4-13).  Equation  4-5  is  not  valid  for  the  small  angles  of  our  data  because 
the  population  diamond  extends  over  a  wide  range  of  velocities 
perpendicular  to  E£  and  Eg.  In  this  case,  is  given  by 


Aw+r/k  VL+2r/ksin5 


Nex  -  Ntot(M/2’rRT) 


/nurri  /  is, 

/ 

Aw-r/k  Vt 


y  X 


where 


=  j  (Aa)-r)/k[tan(0/2)  -  {  v^-(Aw-r)/k[cot  e 


(4-6) 


Note  that  F  should  be  replaced  by  l/2AWpjj  if  power-broadening  is 
contributing  to  the  homogemeous  width.  calculated  from  Equation  4-6 

must  still  be  reduced  by  the  effect  of  grating  washout.  The  washout 


effect  is  also  reduced  substantially  because  the  grating  spacing  is 
increased  and  because  the  grating  lifetime  can  now  be  controlled  by  the 
radiative  lifetime  of  the  excited  state.  In  sodium,  radiative  decay 
(upper  state  lifetime  t  -  16  ns)  is  a  stronger  grating  loss  mechanism 
than  grating  washout  for  angles  smaller  than  about  50  mrad;  the  grating 
steady  state  population  is  determined  primarily  by  radiative  lifetime  for 
angles  less  than  about  15  mrad. 

The  above  model  would  predict  that  reflectivity  in  homogeneously 

-  0 

broadened  lines  would  fall  off  as  (sin  $)~  due  only  to  grating 
washout  because  there  would  be  no  angular  dependence  to  the  participating 
population,  Calculations  using  Equation  4-6  would  indicate  that 

similar  behavior  would  be  expected  in  the  experiment  reported  here  because 
dNex/d^  is  essentially  zero  if  <  hWpjj/Awjj,  the  ratio  of  the 
power -broadened  homogeneous  width  to  the  inhomogeneous  width.  This 
conclusion  is  supported  by  the  numerical  computation  of  Nilsen  and  Yariv, 
presented  as  Figure  4-9,  which  shows  an  approximate  (sin 

dependence  for  angles  in  our  data  range.  This  calculation,  while  generic, 
does  fit  closely  parameters  for  sodium  at  temperature  used  in  our 
experiments.  It  clearly  shows  the  decreasing  importance  of  grating 
washout  at  small  angles  as  the  grating  spacing  increases  and  radiative 
processes  begin  to  control  the  grating  lifetime  and  steady- state 
population.  It  also  begins  to  show  the  effects  of  nonzero  dN^^^/d^  at 
the  larger  angles,  where  the  calculation  begins  to  fall  off  faster  than 
(sin  sy^. 

Our  experiments  were  performed  in  the  Rabi  regime,  so  the  pertinent 
homogeneous  width  was  Aw  .  ,  which  was  typically  0.2Awjj.  Thus  the 
observed  (sin  8)  for  angles  larger  than  “15  mrad  is  as  expected. 

For  smaller  angles,  the  reflectivity  is  observed  to  exhibit  a  stronger 
angular  dependence  rather  than  rolling  off  as  grating  washout  becomes 
unimportant.  This  may  be  an  experimental  artifact  arising  from  changes  in 
the  effective  interaction  volume  defined  by  the  overlap  of  Gaussian  signal 
and  pump  beams,  or  may  be  due  to  residual  pump  misalignment  giving  rise  to 
an  angular  phase  mismatch.  It  can  be  easily  shown  that  even  with  pump 
misalignment,  the  phase  mismatch  in  a  DFWM  experiment  goes  to  zero  when 
Ej  and  E^  are  colinear.  The  impact  on  an  experimental  measurement 
depends  on  the  acceptance  angle  of  the  conjugate  beam  detector. 
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Figure  4-9. 


-2 

Predicted  sin  6  dependence  compared  to  numerical 
calculation  of  FWM  reflectivity  in  the  Doppler  limit. 
The  numerical  calculation  follows  predicted  behavior 
at  small  angles. 
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The  simple  model  presented  above  must  be  considered  In  more  detail  In 
light  of  the  discussion  of  saturation  and  off-resonance  contributions 
(Section  3.1.1  and  4.3.3),  as  well  as  the  effect  of  optical  pvunplng  when 
there  Is  strong  overlap  of  neighboring  hyperflne  transitions  (Section 
3.1.2).  These  effects  are  Important  In  our  experiments  and  will 
substantially  alter  the  velocity  distributions  that  must  be  Integrated 
over  In  Equation  4-6.  Optical  pumping  modification  of  the  velocity 
profile  has  been  considered  by  Humphrey,  et  al.  (Reference  4-5),  whose 
theoretical  treatment  shows  that  optical  pumping  can  account  for  the  (sin 
dependence  they  observed  for  angles  up  to  1  rad  in  a  low-power, 
Doppler -limit  experiment.  However,  those  experiments  were  performed  using 
a  magnetic  field  to  split  the  magnetic  sub-level  degeneracies;  in  the 
absence  of  the  field  they  report  a  (sin  tf)'^  dependence  at  these  large 
angles.  This  change  In  angular  dependence  with  the  application  of  a 
magnetic  field  of  a  few  Gauss  was  unexplained. 

The  above  model,  in  conjunction  with  Equations  3-9,  3-11,  3-12,  and 
4-6,  can  be  used  to  estimate  the  angular  dependence  of  R  for  different 
values  of  I  and  Aw,  once  the  effects  of  optical  pumping  and  strong 
saturation  have  been  considered.  A  few  basic  observations  can  be  made, 
however,  concerning  the  field-of-view  in  the  high  reflectivity  regime  of 
our  experiments.  First,  large  Aw  is  necessary,  which  effectively 
eliminates  the  contribution  of  the  back  grating  formed  by  Ej^  and  E^. 

Under  low  power,  on- resonance  conditions,  this  grating  becomes 
Increasingly  Important  as  Increases,  mitigating  the  loss  of 
reflectivity  from  the  forward  grating.  Secondly,  the  Increased 
homogeneous  width  due  to  power -broadening  leads  to  Increased  and 
larger  6  before  becomes  Important,  but  also  leads  to  more 

rapid  grating  decay  because  the  average  cross-grating  velocity,  <Vgjj>, 

Is  Increased  (see  Equation  3-9  and  accompanying  discussion).  For  small 
9,  <Vgd>  “  ■'^thermal  independent  of  F  (or  AWp^)  so  one  would 
expect  that  the  field  of  view  would  increase  with  F,  as  predicted  by 
Nilsen  and  Yariv.  However,  strong  overlap  of  the  hjrperfine  components 
coupled  with  optical  pumping  (c.f.  Figure  3-lOb)  can  effectively  reduce 
and  Increase  angular  sensitivity  at  the  higher  intensities  required 
for  good  FWM  reflectivity. 


59 


4.4  BROADBAND  REFLECTIVITY 

The  reflectivity  experiments  described  so  far  were  performed  with  a 
laser  bandwidth  of  »  1  MHz ,  narrow  compared  to  any  of  the  characteristic 
frequencies  of  the  sodium  medium.  Such  a  source  bandwidth  is  not  easily 
achieved  under  laboratory  conditions ;  a  requirement  for  such  a  bandwidth 
would  have  significant  impact  on  the  practicality  of  fielding  a  sodium 
vapor  FWM  for  mission  applications.  To  ascertain  the  effects  of  wider 
source  bandwidths.  FWM  reflectivity  was  measured  using  the  broadest 
bandwidth  available  from  our  laser  source. 

The  ring  dye  laser  when  operated  in  the  broadband  configuration  has 
a  frequency  bandwidth  of  2  GHz  consisting  of  a  nvimber  of  longitudinal  loop 
modes  spaced  by  -90  MHz.  The  coherence  length  of  15  cm  was  long  compared 
to  the  interaction  lengths  in  our  experiments  so  the  modes  were  mutually 
coherent.  Even  though  the  broadband  laser  power  (1.25W)  is  greater  than 
the  narrowband  power  (800  mW) ,  the  broadband  power  is  distributed  over  -22 
longitudinal  modes,  resulting  in  a  power  per  mode  that  is  -14  times  less 
than  for  single  frequency  operation.  Changes  in  the  output  beam  diameter 
result  in  an  intensity  -20  times  lower  per  optical  mode.  The  experimental 
setup  was  the  same  as  for  the  narrowband  reflectivity  measurement,  except 
that  meaningful  frequency  measurements  of  the  broadband  laser  were  not 
possible  because  of  inherently  unstable  operation  of  the  laser.  Because 
of  this  instability,  it  was  not  possible  to  measure  the  spectral 
dependence  of  broadband  reflectivity. 

Surprisingly  high  reflectivity  was  observed.  Reflectivity  was 
measured  to  be  -15-20%  in  a  3  cm  cell  containing  sodium  in  a  1  torr  buffer 
of  Ar.  A  6%  reflectivity  was  measured  in  a  1  cm  cell  with  no  added  buffer 
gas.  To  first  order,  R  «  L  so  the  reflectivity  in  the  3  cm  cell 
appears  to  be  diminished  by  the  buffer  gas;  such  a  result  is  consistent 
with  our  observations  in  narrowband  experiments  and  the  results  of  other 
workers  (Reference  4-7).  An  interesting  and  unexplained  phenomenon  was 
that  the  best  reflectivity  in  the  case  of  the  3  cm  cell  was  observed  to 
occur  when  the  two  counterpropagating  pump  beams  were  deliberately 
misaligned.  These  results  can  be  explained  qualitatively  in  terms  of  the 
differences  in  broadband  and  narrowband  dye  laser  operation.  These 
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differences  affect  relative  contributions  to  the  DFWM  signal  by  various 
velocity  groups  of  the  Doppler -broadened  absorption  line. 

One  of  the  main  differences  between  broadband  and  narrowband  DFUM 
Is  that  the  level  of  saturation  per  optical  mode  will  be  smaller  for 
broadband  operation  by  a  factor  of  -20.  Since  the  pump  Intensity/satur¬ 
ation  dependence  of  the  DFWM  signal  (Section  4.3.1)  can  be  quite  strong, 

It  would  be  expected  that  the  broadband  reflectivity  will  be  much  less 
than  the  narrowband  reflectivity.  I/lgat  each  of  the  broadband 
modes  Is  a  factor  of  20  less  than  the  narrowband  single  mode  but 

there  are  —22  of  these  broadband  modes.  These  two  factors  do  not  cancel 
one  another,  however,  because  the  FVM  signal  Is  proportional  to  the  square 
of  the  polarization  and  the  sum  of  the  squares  of  each  of  the  22  broadband 
polarization  components  (which  assumes  mutual  coherence  among  the  modes) 
will  be  less  than  the  square  of  the  single  (large)  narrowband  polarization 
component.  For  the  numbers  In  our  experiments,  this  simplistic  argument 
would  Indicate  that  the  broadband  reflectivity  would  be  about  6%  of  the 
narrowband  reflectivity,  remarkably  within  a  factor  of  two  of  the  observed 
ratio  of  about  3%  (6%/230%)  tinder  high  reflectivity  conditions  In  the  1  cm 
cell.  Stronger  self- focussing  In  the  narrowband  experiments  could  account 
for  this  difference.  A  more  complete  analysis  would  require  Incorporating 
the  effects  of  a  nonuniform  distribution  of  optical  energy  over  the 
longitudinal  modes  and  consideration  of  saturation  effects.  These 
considerations  would  require  detailed  Information  on  the  actual  frequency 
spectrum  of  the  laser  line  with  respect  to  the  Doppler-broadened  D2  line 
before  a  more  careful  comparison  between  broadband  and  narrowband 
reflectivities  can  be  made.  We  were  not  able  to  collect  this  Information 
with  our  present  experimental  configuration. 

Another  consequence  of  the  broadband  power  distribution  is  that 
contributions  to  the  signal  will  now  be  possible  from  non- degenerate 
four -wave  mixing  (NDFM)  processes.  Because  of  the  range  of  optical 
frequencies,  the  resonance  conditions  described  by  Equations  3-3 

and  3-4  In  Section  3.1.2  can  be  satisfied  over  a  wide  range  of  velocities 
for  both  the  forward  and  backward  gratings.  The  nearly  co- linear  E£  and 
Eg  beams  interact  with  the  same  sodium  velocity  group  to  write  an 
efficient  refractive  index  grating  at  a  detuning  A(i)  from  the 
Doppler- line  center.  For  very  small  6,  these  two  fields  will  generally 
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correspond  to  the  same  longitudinal  optical  mode  in  both  the  signal  and 
forward  piimp  fields  because  they  must  be  within  a  homogeneous  width  of 
each  other  to  write  the  interference  grating.  For  our  experiments,  the 
power  broadening  under  broadband  conditions  was  such  that  only  one  or  two 
modes  could  Interact  to  form  a  grating.  (This  limitation  in  the  number  of 
modes  writing  the  grating  holds  at  larger  angles  also,  but  the  mode  may 
not  be  the  same  in  the  signal  and  forward  pump  fields  because  they 
experience  different  Doppler  shifts.  This  "mode-hop” ,  where  it  is  the  n+1 
mode  of  the  signal  field  writing  a  grating  with  the  n^^  mode  of  the 
forward  pvuop,  occurs  at  such  large  angles  in  our  experiments  that  such 
interactions  do  not  need  to  be  considered  here.) 

Unlike  the  case  of  DFUM  discussed  in  Section  3.1.2,  there  are  now  ' 
many  optical  modes  in  for  reading  out  the  forward  grating.  Consider 
a  case  where  the  homogeneous  width  allows  interaction  of  only  one  optical 
mode  from  each  field  with  a  particular  velocity  group.  At  small  angles 
when  the  grating  is  written  by  the  same  optical  mode  in  Ej  and 

E^,  a  phase-matched  DFWM  experiment  occurs  for  the  mode  in  E^^. 

This  mode  sees  the  nonlinearity  Doppler- shifted  by  2(Wjj-wq).  All 
other  modes  in  Ejj  participate  in  NDFWM  experiments.  Figure  4-10  shows 
that  these  NDFWM  interactions  are  also  phase  matched  if  For 

nonzero  the  phase  mismatch  for  mode  ii  «  is 

Ak  «  (|/c)[(l-cosfljj)/costfm]  (^*7) 

where  8^-  arcsin  {WjjSintf£g/((i>jj+|) )  »  and  $  is  the  frequency 

difference  of  the  two  modes  in  the  NDFWM  experiment.  In  our  experiments, 
the  maximum  value  of  i  is  2x(2  GHz),  which  for  a  8^^  of  8  mrad, 
results  in  excellent  phase  matching  for  all  modes  over  the  1  cm  inter¬ 
action  length.  Thus  all  modes  in  E^  are  backscattered  in  phase-matched 
FWM  experiments.  The  efficiency  for  a  particular  mode,  £,  is  deter¬ 
mined  by  how  much  off  resonance  it  reads  the  x  grating  generated  by  E^ 
and  Eg  at  the  linear  absorption  at  w^,  and  the  effects  of 
contrast  and  saturation  discussed  in  Section  3.1.  For  example,  given  that 
E£  and  Eg  are  detuned  by  6u,  the  back  pump  mode  at  Wjj-2Aw 
reads  the  x  grating  on  resonance,  but  because  of  the  lineshape  of  the 
D2  line  tends  to  be  strongly  absorbed. 
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Figure  4-10.  Frequency  and  phase-matching  conditions  in  NDFWM  experiments  occurlng  during 
broadband  laser  pumping.  The  experiments  are  strictly  phase  matched  only  for 

0fs  =  0,  but  the  mismatch  at  0^^  +  0  is  negligible  for  the  2  GHz  bandwidth  used 
in  these  experiments. 


At  large  angles  where  "mode-hop"  has  occurred,  the  signal  now  arises 
purely  from  NDFWM  Interactions.  Letting  mode  be  that  of  Ef,  and 
considering  to  be  mode  a  conjugate  beam  at  will  be 

formed  by  scattering  of  the  (back  pump)  mode  In  E^.  The 
appropriate  shift  In  conjugate  wave  frequency  about  Is  a  result  of 
the  Doppler  shift  from  the  Ef/E^  grating,  which  In  this  nondegeneate 
experiment  Is  produced  by  the  traveling  Interference  wave  at 

.  Phase  matching  arguments  are  essentially  the  same  as 
those  presented  above.  Efficiency  considerations  are  similar  to  the  case 
where  the  same  mode  In  E£  and  E^  Interact,  except  that  the  finite 
response  time  of  sodium  becomes  Important  when  nondegenerate  fields  write 
the  gratings.  The  time  dependence  of  the  traveling  grating  must  be  slow 
compared  to  the  response  time  of  sodltjua,  l.e.,  |u>£-(<>g|  <  1/T^  (or 
the  Rabl  frequency),  otherwise  the  grating  contrast  will  become  blurred. 

An  Important  difference  between  hl^  reflectivity  experiments  using 
narrowband  and  broadband  sources  Is  the  contribution  of  gratings  formed  by 
signal  and  back  pump  fields  (  the  "back"  grating) .  The  discussion  In 
Section  3.1.2  and  Figure  3-5  Indicate  that  back  gratings  are  not  Important 
In  narrowband  DFWM  experiments.  This  Is  not  the  case  In  broadband 
experiments,  where  mode  In  the  signal  field  can  Interact  with  the 
same  velocity  group  as  mode  (i)^-2Aw  In  the  back  pump,  producing  a 
grating  that  Is  well  within  Such  gratings  are  "read"  by  all 

modes  In  E£,  and  are  Inherently  NDFWM  Interactions.  Efficiency 
arguments  are  similar  to  those  presented  above,  with  some  modifications. 
First,  because  of  the  finite  response  time  of  sodium,  the  traveling 
grating  written  by  Eg  and  E|j  will  'imear  unless 

2A(j  <  1/T]^.  Thus  only  one  or  two  modes  near  wq  can  result  In 
Important  back  gratings.  These  few  back  gratings  are  still  "read"  by  all 
modes  In  E£,  with  phase  mismatches  that  are  negligible  In  our 
experiments . 

If  sufficient  overlap  of  the  laser  line  with  the  sodium  absorption 
line  exists  so  that  Aw  approaches  zero,  back  gratings  and  their  NDFWM 
contribution  will  be  Important  In  broadband  experiments .  An  example  Is 
fleld-of-vlew.  The  back  grating  has  Its  minimum  spacing  at  ^£g=0, 
where  the  ravages  of  thermal  washout  minimize  Its  contribution  to  the  FWM 


64 


reflectivity.  As  8^^  is  increased,  the  back  grating  contribution  to 
FWM  reflectivity  increases,  offsetting  the  loss  of  reflectivity  from  the 
forward  grating.  Thus  one  would  expect  a  larger  field-of-view  for 
broadband  operation.  Experimental  observation  of  this  is  presented  in 
Figure  4-11,  which  shows  the  angular  dependence  of  broadband  reflectivity 
compared  to  that  of  normalized  narrowband  reflectivity. 

4.5  ACTIVE  TRACKING  DEMONSTRATION 

An  active  tracking  demonstration  was  performed  to  show  the  automatic 
tracking  capabilities  of  a  phase-conjugating  FWM.  Figure  4-12  shows  the 
experimental  setup,  which  consisted  of  the  standard  counterpropagating 
pump  beam  geometry,  except  that  the  back  pvunp  was  generated  by  a 
retroreflectlng  mirror.  The  signal  beam,  which  essentially  constituted  a 
very  bright  glint,  was  split  off  of  the  pump  leg  by  a  beamsplitter  and 
directed  to  a  galvometrically  scanned  rotating  mirror  (RM)  located  at  the 
radius  of  curvature  of  a  spherical  mirror  (SM) .  The  sodium  cell  FWM  was 
also  located  at  the  radius  of  curvature  of  the  spherical  mirror  and  placed 
such  that  the  cell  and  the  rotating  mirror  were  on  opposite  sides  and 
equally  spaced  from  the  optical  axis  of  the  spherical  mirror.  This 
arrangement  permitted  a  substantial  angular  sweep  of  the  incoming  signal 
beam  over  the  surface  of  the  spherical  mirror  while  maintaining 
approximately  the  same  interaction  length  in  the  sodium  cell.  The 
conjugate  return  reflected  off  of  the  spherical  and  rotating  mirrors  and 
was  separated  from  the  incoming  signal  beam  by  a  beam  splitter.  The 
conjugate  was  monitored  by  passing  it  through  an  aperture  to  measure 
divlation  from  the  signal  path,  and  displayed  on  a  white  card.  If  the  FWM 
cell  provided  accurate  automatic  tracking,  the  conjugate  beam  would 
retrace  the  signal  beam  path  precisely,  and  since  the  rotating  mirror  was 
essentially  stationary  on  the  timescale  of  the  optical  roundtrip  from  it 
to  the  sodium  cell,  the  conjugate  beam  would  pass  through  the  aperture  and 
be  visible  on  the  display  card. 

Automatic  tracking  of  the  input  signal  beam  was  observed  for  angular 
scans  of  50  mrad,  a  llmication  Imposed  by  the  size  of  the  spherical 
mirror.  The  rotating  mirror  was  driven  at  1000  Hz,  providing  an  angular 
scan  rate  in  excess  of  100  radians/s.  Under  all  conditions,  the  conjugate 
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Figure  4-11.  Observed  angular  dependence  of  DFWM  reflectivity  using 

broadband  (2  GHz)  laser  irradiation.  Data  was  collected 
in  the  same  manner  as  narrowband  reflectivity  (also  shown) , 
except  that  data  was  recorded  photographically  to  overcome 
frequency  instability  of  laser  during  broadband  operation. 


return  was  observed  to  pass  through  the  aperture.  Angular  deviation  was 
less  than  1  prad.  which  was  the  precision  of  our  measurement.  This 
experiment  demonstrates  that  while  FWM  reflectivity  drops  off  very  rapidly 
in  sodium  vapor,  the  effective  field-of-view  for  automatic  tracking  can  be 
substantially  larger  if  the  FWH  reflectivity  is  adequate. 
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5.0  CROSSTALK  ANALYSIS  AND  EXPERIMENTS 


Crosstalk  between  a  phase  conjugate  wave  and  noise  waves  was  investi¬ 
gated  experimentally  to  evaluate  detectivity  limits  intrinsic  to  the  FWM 
process.  Crosstalk  was  observed,  but  only  under  high  probe/noise  condi¬ 
tions  where  pump  depletion  occurs .  These  measurements ,  which  fully 
satisfy  the  contract  SOW,  thus  show  that  crosstalk  is  insignificant  under 
the  weak  probe  conditions  likely  to  be  encountered  in  phase  conjugated 
optical  links.  It  was  determined  that  analysis  of  FWM  physics  beyond  that 
required  by  the  SOW  was  desireable,  however,  so  a  theoretical  evaluation 
of  crosstalk  was  undertaken  on  TRW  IR&D.  This  evaluation  is  included 
because  it  substantiates  our  experimental  conclusions. 

5.1  OBJECTIVES  AND  BACKGROUND 

The  detectivity  of  a  signal  on  a  phase  conjugate  wave  can  be  limited 
by  crosstalk  between  the  phase  conjugate  wave  and  one  or  more  noise  waves 
in  the  FWM  cell.  In  this  section  we  shall  examine 

a.  intensity  crosstalk,  i.e.,  unwanted  temporal  intensity  varia¬ 
tions  ,  and 

b.  phase  crosstalk,  i.e.,  unwanted  effects  of  a  noise  wave's 
transverse  phase  distribution. 

In  Section  6,  spatial  crosstalk,  i.e.,  unwanted  transverse  Intensity 
patterns,  is  examined  and  shown  to  derive  from  the  intensity  crosstalk 
mechanisms  identified  in  this  section. 

We  start  with  a  survey  of  possible  crosstalk  mechanisms  such  as  pump 
depletion,  probe  depletion,  etc.,  as  well  as  some  observations  about  what 
might  be  expected  in  the  Na  vapor  experiments  to  be  discussed  next.  For 
each  experiment  there  is  a  description  of  its  design,  followed  by  a 
presentation  of  the  results  and  conclusions. 

It  will  be  seen  that,  for  combined  probe -noise  beam  intensities 
that  are  low  to  moderate  compared  to  the  total  pump  intensity,  the 
crosstalk  appears  to  be  due  principally  to  a  combination  of  pump  depletion 
and  intensity  saturation  of  the  FWM  interaction.  At  higher  probe-noise 
Intensities,  there  is  evidence  of  contributions  from  other  mechanisms.  An 
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empirical  expression  that  fits  the  data  on  probe  reflectivity  vs.  probe 
intensity  can  be  used  to  predict  the  intensity  crosstalk  for  arbitrary 
probe/noise  intensity  ratios.  It  will  be  shown  experimentally  that  the 
predictions  are  quite  accurate.  This  same  expression  will  be  used  to 
predict  the  tolerable  noise  intensity  (or,  conversely  the  necessary  probe 
intensity)  in  a  system  for  communicating  via  intensity  modulation. 

5.2  CROSSTALK  MECHANISMS 

To  describe  FWM,  we  start  with  Maxwell's  equations  and  eliminate  the 
B  field  to  obtain  the  wave  equation  for  the  electric  field.  If  all 
electric  fields  are  copolarized  the  result  in  cgs  units  is 

«  4ir 

V^E  - - . 

c^  at2  dt^ 

where  Pjjl  is  the  nonlinear  part  of  the  induced  polarization.  There  are  at 
least  four  waves  present  (the  two  pumps,  the  probe  and  the  probe 
conjugate) ,  so  we  make  this  explicit  by  writing 

In  the  slowly  varying  wave  approximation,  the  expression  for  each  field  is 

Em  -  1/2  Eni(r,t)  exp[i(^*r  -  Wn,t)]  +  c.c., 

where  c.c.  is  the  complex  conjugate  and  Em(r,t)  is  assumed  to  vary  slowly 
with  respect  to  r  and  t  compared  with  the  rapid  variation  of  the  exponen¬ 
tial.  For  FWM  in  a  Kerr-like  medium,  the  nonlinear  polarization  has  the 
form 

PnL  -  e3, 

where  Is  the  third  order  susceptibility.  In  steady  state,  the  slowly 

varying  envelope  approximation  to  the  wave  equation  is 

^Eni(r)  iJTWjj  __  _ 

— di - 2^^^  ^  Ej(r)Ei,(r)E^(r)  exp[i(K-k^).r]  (5-la) 

where  wj  +  -  Wjj,,  (5 -lb) 

and  K  ■  icj  +  -  k]^,  (5-lc) 


and  where  the  z-axls  is  parallel  to  kjj.  Strictly  speaking,  all  combina- 


tions  Ej ,  E^,  and  Ej^,  which  are  any  three  of  the  waves  composing  the  total 
field,  must  be  included  as  source  terms  on  the  right  hand  side  of  Eqn. 

5 -la,  but  very  few  of  the  possible  triplets  make  significant  contribu¬ 
tions.^  Unless  kjj  is  very  close  to  K,  the  rapidly  varying  exponential 
causes  E^,  to  be  very  small.  Thus  k^,  must  be  parallel  to  K,  and  we  must 
have  |K|  »  w^/nc,  ie.,  there  must  be  very  little  phase  mismatch.  For 
degenerate  FWM,  —  w,  the  input  frequency,  for  all  m. 

Eqn.  5 -la  shows  that  the  element  of  wave  dEjj(r)  created  in  the 
interval  dz  is  proportional  to  the  three  contributing  waves,  Ej ,  Ej^,  and 
E]^.  Furthermore,  at  any  point  (x,y)  in  the  plane  normal  to  the  z-axis, 
dEjj(x,y,z)  is  proportional  to  the  magnitudes  of  these  three  waves  at  that 
point.  This  behavior  will  be  important  in  the  understanding  of  transverse 
spatial  crosstalk. 

Eqn.  5-la  also  says  that  any  trio  of  existing  waves,  E j  ,  E]^, 

whose  new  wave  vector  |K|  =  w/nc,  may  give  rise  to  a  new  wave;  if  K  is 
equal  to  the  wave  vector  of  an  existing  wave,  then  dE^  will  be  an  incre¬ 
ment  to  that  wave.  The  crosstalk  experiments  were  done  near  absorption 
resonances  in  Na  vapor,  so  the  medium  was  not  Kerr-like.  Instead  the 
susceptibility  is  intensity  dependent  and  the  induced  polarization  at  any 
point  vanishes  as  the  total  intensity  at  that  point  becomes  arbitrarily 
large  (Ref.  5-1,  5-2,  5-3  and  discussion  in  Section  3.1).  Consequently, 
the  fields  radiated  by  the  oscillating  polarization  decrease  for  suffi¬ 
ciently  large  input  intensity.  This  is  reflected  in  the  replacement  of 
the  constant  (iww/2nc)x^^^  in  Eqn.  5-1  by  Wjj^]^(I),  a  function  of  the  total 
intensity  in  the  interaction  region,  ^jkl^^^  decreases  for  very  large 
total  intensity.  This  saturation  behavior  will  later  be  seen  to  be  a 
significant  factor  contributing  to  crosstalk. 

The  FWM  generation  of  the  probe -conjugate  wave  Ep^,,  and  the  amplifica¬ 
tion  of  the  transmitted  probe  Ep  are  described  by 

dEpc  ”  ^fbp^^^  ^f^b^p  (5-2a) 


^^p  ”  ^fbpc^^^  ^f^b^'pc 


(5-2b) 


^To  simplify  the  notation,  £  will  be  represented  as  E  from  here  on. 


where  Ef  and  are  the  forward  and  backward  pump  waves .  Now  we  can 
anticipate  that  if  the  total  intensity  I  is  already  in  the  saturation 
regime,  then  a  noise  wave,  strong  enough  to  significantly  increase  I,  will 
decrease  the  probe -conjugate  intensity,  i.e.,  it  will  cause  crosstalk. 

Terms  describing  absorption  or  gain  should  be  included  in  Eqn.  5-2: 
they  are  ignored  here  because  we  are  primarily  interested  in  terms  that 
create  new  waves.  Also  ignored  are  focusing  terms  that  are  proportional  to 

f  /I?  \  T?  /  T?  \  T-  /  m  \  f-  _ 3  /  T-  \  T- 


coupling  triplets  of  the  form  (EfEf)Ep,  (EfEf)Epj.,  (EbEb)Ep,  and  (EbEb)Epc. 


5.2.1  Intensity  Crosstalk 


First  we  consider  FWM  waves  that  are  generated  in  the  absence  of  noise 
beams.  Two  waves,  dEj^  -  EpE^jEf  with  wave  vector  kp-2kf  and  its  conjugate 
dEj^j.  -  EpjjEjE^  with  wave  vector  2kf-kp,  are  created  with  strong  coupling 
since  each  triplet  contains  two  strong  pump  waves.  There  is,  however, 
some  phase  mismatch  because  our  probe  beam  was  at  an  angle  of  the  order  of 
10  mrad  from  the  pump  beams ,  so  that  the  above  wave  vector  magnitudes  were 
not  exactly  equal  to  w/nc.  Thus  the  intensities  of  these  waves  are 
expected  to  be  considerably  smaller  than  those  of  the  probe  and  probe - 
conjugate. 

Other  weaker  waves,  created  by  triplets  with  only  one  strong  pump 
term,  are  dE2  -  EpEjjE^^,  with  wave  vector  2kp-kf  and  its  quasi -conjugate^ 
dE2c  ~  Epc^f^p  (wave  vector:  k£-2kp);  again  there  is  an  angular  phase 
mismatch.  These  two  now  give  rise  to  a  two-step  infinite  series  of  waves 
and  quasi-conjugates : 


Wave  Vector 


dE„.  -  Em-lEbEf 

‘^^m  ~  ^-l^p^pc 
‘^Ejiic  ~  ^(m-l)c^b^f 
^^c  ~  ^(m-l)c^pc^p 


(m-l)kp-mk£ 
mkp- (m-l)kf 
mkf- (m-l)kp 
(m-l)kf-mkp 


m  odd  3 . . . « 


m  even  4.  .  .<» 


m  odd  3 . . . « 


m  even  4. . .« 


The  phase  mismatch  for  each  succeeding  generation  increases,  thus  reducing 
the  power  fed  from  one  generation  to  the  next.  In  addition,  this  phase 
mismatch- induced  decrease  in  |E„,[  relative  to  lE^,.]^!  itself  enhances  the 
relative  decrease  of  |  |  .  The  result  is  a  rapid  decrease  in  the  power 


-If  the  pump  waves  are  exactly  plane,  these  will  be  true  conjugates. 


channeled  into  higher  generation  waves,  Ej^.  The  above  considerations 

indicate  that  at  fixed  input  pump  intensities,  the  probe -conjugate 

# 

intensity  will  first  Increase,  and  then  saturate  as  the  input  probe 
intensity  increases.  There  are  three  causes  of  saturation: 

1.  There  is  an  intrinsic  saturation  built  into  the  coupling 

coefficient  . 

2.  The  pump  depletion  necessary  to  provide  energy  to  the  conjugate 
beam  decreases  the  coupling  triplet. 

3.  The  pump  depletion  necessary  to  create  the  higher  generation 
waves  also  decreases  the  coupling. 

It  will  be  seen  later  that  a  quantitative  description  of  the  saturation  of 
the  probe -conjugate  intensity  allows  us  to  predict  the  Intensity  crosstalk 
accurately. 

With  the  introduction  of  a  noise  wave  E^,  many  new  interactions  occur 
that  can  in  principle  contribute  to  crosstalk  because  they  lead  to 
depletion  of  the  pvimp,  probe,  and/or  conjugate  beams.  First,  in  a  manner 
analogous  to  that  for  the  probe  beam,  the  pump  beam  is  depleted  by  the 
creation  of  a  noise -conjugate  beam  and  amplification  of  the  transmit¬ 

ted  noise.  Second,  there  is  creation  of  waves  analogous  to  E]^  and  Ej^j.- 
e{  and  its  conjugate  E^^.,  where  Ep  and  Epg  are  replaced  in  the  coupling 
triplets  by  Ejj  and  E^^,  respectively.  Third,  there  are  completely  new 
interactions  capable  of  producing  crosstalk  that  are  not  present  in  the 
absence  of  a  noise  beam.  For  example,  both  the  pump  and  either  the  probe 
or  the  probe  conjugate  are  depleted  by  the  creation  of  new  waves  via  the 
following  coupling  triplets : 


Triplet 

Wave  Vector' 

EbEpE^ 

kp"kf  *kjj 

EfEpE^ 

kp+kf-kjj 

EbEpEnc 

kp'kf+kjj 

EfEpEjic 

tp+kf+i^n^ 

EbEpcEn 

kp*kf-kjj 

EfEpcE?; 

-kp+kf-kjj 

EbEpcEnc 

■kp'kf-kjjt 

EfEpcEnc 

-kp+kf+kjj 

The  coupling  for  this  third  category  of  waves  is  weaker  because  only 
one  factor  in  the  triplet  is  a  pump  wave;  also  two  of  the  8  terns  have  a 
large  phase  mismatch  and  the  remainder  have  small  mismatches.  Therefore, 
for  probe  and  noise  Intensities  small  compared  to  that  of  the  pump,  they 
cannot  be  as  significant  in  crosstalk  as  the  pump  depletion  due  to 
generation  of  the  noise  conjugate,  dEjjg  ~  and  amplification  of 

the  noise  beam,  dEjj  ~  EjE^E^j,. 

The  probe  and  its  conjugate  can  also  be  directly  depleted  through 
coupling  triplets  that  lead  to  growth  of  the  noise  and  noise-conjugate 
waves:  dEjj  -  EpEp^E^J^,  (wave  vector:  icjj)  and  dEjj^  ~  EpEp(.E^  (wave  vector: 
-icjj)  .  This  coupling  is  weak,  because  there  are  no  pump  waves  in  the 
triplet,  but  the  phase  match  is  perfect.  It  will  be  seen  later  that  our 
experimental  results  eliminate  these  two  temns  as  significant  sources  of 
crosstalk  in  the  range  of  probe-plus-nolse  intensities  we  studied. 

There  are  also  Infinite  sequences  of  waves  created  in  the  presence  of 
the  noise  wave  such  as: 

dEj  -  EbEpE*  kp-kf-i^ 

-  E(ni.i)EpEf  m[kp-kf]-^  m  integer  2...® 

In  general,  some  higher  generation  waves  are  simply  higher  order  contribu¬ 
tions  to  an  earlier  generation  wave,  and  therefore  end  the  sequence. 

Some,  like  those  just  listed,  are  part  of  an  infinte  sequence  in  which  the 
phase  mismatch  gets  increasingly  large,  and  the  coupling  terms  contain 
increasingly  smaller  fields. 

^Wave  vectors  with  three  terms  of  the  same  sign,  marked  with  f.  have  a 
momentum  magnitude  of  3u/nc  and  a  very  large  phase  mismatch. 


Ve  will  show  latter  in  this  section  that  most  of  our  crosstalk 
results  can  be  explained  by  Invoking  only  pump  depletion  caused  by  growth 
of  the  probe,  probe • conj ugate ,  noise,  and  nolse-conjugate  beams;  and  by 
the  susceptibility  saturation  Inherent  In  resonant  media.  At  very  high 
combined  probe  and  noise  Intensities,  however.  It  may  be  necessary  to 
attribute  part  of  the  crosstalk  to  the  creation  of  new  waves  that  can 
exist  only  In  the  presence  of  a  noise  wave. 


5.2.2 


Phase  crosstalk,  that  is,  phase  variations  on  that  are  imparted  to 
Epg,  can  be  studied  by  examining  the  coupling  terms  that  contribute  to  the 
creation  of  Ep^;  the  conclusion  below  will  be  that  no  phase  noise  is 
transmitted. 


In  the  first  generation,  we  have  the  pumps,  probe,  and  noise  beams. 
Assume  that  an  element  of  some  second  generation  wave  E^  results  from 
the  coupling  of  three  first  generation  waves;  dE4  -  E]^E2E3.  Then  an 
element  of  a  third  generation  wave  results  from  the  coupling  of  three 
second  generation  waves  Including  E4:  dE7  ~  E4E5E5.  Substituting  for  E4. 
we  obtain  dEy  -  E]^E2E3E5Eg.  In  general,  an  element  of  any  higher  genera¬ 
tion  wave  is  proportional  to  all  of  the  lower  generation  waves  that  led  to 
its  creation.  At  each  stage  the  number  of  factors  In  the  coupling  term  is 
three,  so  the  final  number  of  factors  must  be  odd,  as  in  the  expression 
for  dEy  above.  Further  the  array  of  factors  in  any  contribution  to  dEp^, 
must  be  of  the  form, 

(EfE*)  (EjEJ)  ...  (EkEj^)ExE.xEj 

where  E^  and  E.^  are  waves  travelling  antiparallel  to  each  other,  such  as 
EfEb,  or  EpEp^,  or  EjjEjj^..  First,  the  nvunber  of  luiconjugated  (unstarred) 
factors  must  be  one  greater  than  the  number  of  conjugated  (starred) 
factors  for  the  wave  element  to  have  the  frequency  u  appropriate  to  Ep^ . 
Second,  each  wave  other  than  E^  or  E.^  must  be  paired  with  its  conjugate 
in  order  that  the  resultant  wave  vector  be  that  of  Ep^,;  the  result  is  that 
all  wave  vectors  cancel  out,  except  that  of  E^.  In  the  above  expression, 
then,  we  see  that  if  the  noise  beam  appears,  it  does  so  only  in  the  pairs 
^^nc‘  either  case,  its  phase  distribution  cancels  out  and  is 
not  Imparted  to  the  probe -conj ugate  wave. 


In  Na  vapor  the  angular  response  of  FWM  is  limited  by  thermal  motion 
of  the  atoms.  Particularly  in  cw  experiments,  atomic  motion  washes  out 
the  volume  diffraction  gratings  formed  by  the  interference  of  the  probe 
and  pump  waves ,  unless  there  is  only  a  very  small  angle  between  the  two 
corresponding  wave  vectors.  This  is  discussed  in  Section  4.3.4.  The 
result  is  that  there  is  essentially  no  reflectivity  outside  some  small 
angle.  Similarly,  a  noise  beam  cannot  reflect  outside  of  this  angle,  and 
therefore  cannot  cause  pump  depletion  at  larger  angles.  Thus  we  expect  a 
rapid  decrease  in  crosstalk  as  angular  separation  between  the  noise  and 
pumps  Increases;  at  larger  angles,  however,  we  expect  some  residual 
crosstalk  due  to  the  intrinsic  saturation  of  the  susceptibility  in 
resonant  interactions. 

5.2.4  Crosstalk  Frequency  Dependence 

An  off-resonant  noise  beam  whose  frequency  is  outside  the  Na  linewldth 
will  not  reflect  in  the  Intensity  range  considered  here.  Furthermore,  if 
the  noise  beam  frequency  la  offset  by  more  than  the  Na  vapor  homogeneous 
linewldth  from  the  pximp  frequency,  it  will  be  poorly  reflected  (See 
Section  3.1).  Therefore  it  cannot  deplete  the  pump  and  cause  crosstalk. 

We  anticipate  that  the  crosstalk  caused  by  a  noise  beam  at  a  given 
frequency  will  be  correlated  with  its  phase  conjugate  reflectivity  at  that 
frequency . 

5.3  INTENSITY  CROSSTALK  EXPERIMENTAL  DESIGN 

The  effect  of  a  controllable  noise  source  on  the  probe  conjugate 
intensity  was  initially  studied  with  the  experimental  arrangement  in 
Figure  5-1.  Two  counterpropagating  beams  of  equal  intensity  were  produced 
by  beam  splitter  BSl  and  overlapped  in  the  Na  cell.  A  slight  misalignment 
in  the  pump  beams  prevented  feedback  to  the  ring  dye  laser  and  was 
essential  in  maintaining  a  stable  single  frequency  output.  A  1.0  m 
focal  length  lens  LI  was  placed  before  BSl  to  impart  equal  convergence  to 
both  pump  legs.  Pump  intensities  at  the  Na  cell  were  typically  60-70 
W/cm^  after  correction  for  cell  window  transmittance.  The  wedged  beam¬ 
splitter  UBS  produced  two  separate  beams;  the  first  surface  reflection  was 
used  to  diagnose  the  frequency  output  of  the  dye  laser  by  directing  it 
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Figure  5-1.  Intensity  crosstalk  experiment. 


Into  a  Burleigh  wavemeter  (»  1  GHz  frequency  resolution) ,  or  when  greater 
frequency  resolution  was  required,  Into  the  Lamb  dip  diagnostic  described 
In  Section  4.1.  The  stronger  second  surface  reflection  was  used  to 
generate  both  the  probe  and  noise  beams.  With  a  fixed  total  power  of  both 
beams  established  by  WBS,  the  setting  of  the  variable  beamsplitter  VBS 
determined  the  nolse*probe  Intensity  ratio. 

The  probe  beam  was  passed  through  beamsplitter  BS3,  focusing  lens  L3, 
and  directed  Into  the  Na  cell  at  a  fixed  angle  tp  with  respect  to  the  pump 
beam  axis.  was  kept  at  the  minimum  allowed  by  finite  beam  diameters  and 
was  measured  to  be  6.6  mrad  for  the  arrangement  diagrammed  In  Figure 
5-1.  Variable  noise  beam  angles  were  obtained  by  the  use  of  a  precision 
spherical  mirror  SM  and  an  electromechanlcally  driven  rotatable  mirror  RM. 
Both  RM  and  the  Na  cell  were  located  at  radius  of  curvature  distances  from 
SM  such  that  the  noise  beam  overlap  with  the  pump  and  probe  beams  remained 
constant  at  all  noise  beam  Input  angles  0.  $  was  varied  from  a  minimum  of 

13  mrad  to  24  mrad  with  respect  to  the  pump  beam  axis  while  the  noise  beam 
Itself  was  moved  In  a  direction  perpendicular  to  the  plane  containing  the 
pump  and  probe  beams  (See  Figure  5-2).  The  noise  beam  was  temporally 
chopped  to  detect  small  changes  In  the  probe -conjugate  Intensity  level  In 
Its  presence.  The  probe-  and  noise-beam-conjugate  Intensities  were 
measured  by  photodiodes  PD2  and  PD3  respectively,  and  their  outputs  were 
displayed  on  an  oscilloscope. 

5.4  INTENSITY  CROSSTALK  RESULTS 

Typical  crosstalk  data  appears  In  Figure  5-3  where  S  Is  the  measured 
probe -conjugate  Intensity  level  In  the  absence  of  the  noise  beam  and  AS  Is 
the  decrease  In  the  conjugate  Intensity  level  In  the  presence  of  the  noise 
beam.  The  ratio  AS/S  Is  defined  to  be  the  fractional  crosstalk  experi¬ 
enced  by  the  probe  conjugate  due  to  the  presence  of  the  noise  beam. 

The  lower  trace  In  Figure  5-3  displays  the  conjugate  intensity  of  the 
chopped  noise  beam  which  was  found  to  be  correlated  to  the  decrease 
observed  In  the  probe  conjugate  level.  Only  decreases  In  probe -conjugate 
Intensities  were  observed  In  the  presence  of  noise  beams ;  no  Increases  in 
probe -conjugate  levels  were  discerned  with  the  addition  of  noise  sources. 
This  Is  consistent  with  the  crosstalk  mechanisms  listed  In  Section  5.2, 
l.e.,  saturation  of  the  resonant  susceptibility  and  pump  depletion. 
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Na  CELL  LOCATION 


<t>  -  6.6  mrad 
24  mrad  >d>  13  mrad 


SPHERICAL 

MIRROR 


DIRECTION  OF  NOISE 
BEAM  MOVEMENT 


Figure  5-2.  Geometry  for  Intensity  crosstalk  experiments. 
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Crosstalk  Experiment 


Figure  5-3.  Typical  intensity  crosstalk  data  showing  phase  conjugate  returns  of  the 
probe  and  noise  beams.  The  noise  beam  wds  mechanically  chopped. 


;«y‘. 


Noise-probe  intensity  ratios  of  3:1  and  1:1  were  used  in  the  study  of 
the  effect  of  noise  beam  intensity  and  input  angle  on  probe -conjugate 
Intensity  (See  Figures  5-4  and  5-5).  At  noise-probe  intensity  ratios 
below  1:1,  crosstalk  was  observed  but  was  difficult  to  quantify  due  to  the 
small  change  in  the  probe  conjugate  intensity.  Data  from  Figures  5-4 
and  5-5  are  plotted  in  Figure  5-6  as  percent  crosstalk  versus  input 
noise  beam  angle  from  the  minimum  9  of  13  mrad.  For  both  sets  of  data, 
increases  in  the  noise  input  angle  reduce  the  observed  percent  crosstalk. 
The  3:1  noise-probe  intensity  ratio  produced  consistently  higher  crosstalk 
values  at  all  noise  input  angles  when  compared  to  the  1:1  data  set; 
however,  when  normalized  to  their  respective  maxima  in  Figure  5-7,  both 
data  sets  display  similar  angular  dependences.  A  correlation  between 
noise  beam  conjugate  intensity  and  percent  crosstalk  can  be  seen  in  Figures 
5-8  and  5-9  where  the  normalized  noise -conjugate  Intensity  and  normalized 
percent  crosstalk  have  been  plotted  against  noise  beam  input  angle.  For 
the  3:1  intensity  ratio  data  set  (Figure  5-8),  the  noise  conjugate 
intensity  is  seen  to  decrease  more  rapidly  than  the  percent  crosstalk  with 
Increases  in  noise  input  angle.  A  more  similar  angular  dependence  was 
demonstrated  by  the  1:1  intensity  ratio  data  set  in  Figure  5-9  where 
relative  noise -conjugate  intensity  and  crosstalk  differ  at  most  by  6-7%. 

These  results  suggest  that  the  observed  intensity  crosstalk  is  caused 
by  two  terms:  one  angle  dependent  and  one  angle  independent.  In  Section 
5.2  we  indicated  that  crosstalk  induced  by  pvimp  depletion  should  occur  only 
within  a  small  angular  range  between  the  noise  and  pump  beams ,  namely ,  the 
range  over  which  the  noise -conjugate  Intensity  is  large.  Figures  5-8 
and  5-9  show  that  this  range  is  about  8-10  mrad,  consistent  with  the 
angular  dependence  data  presented  in  Section  4.3.4,  and  Figure  5-6 
confirms  a  sharp  fall-off  in  crosstalk  out  to  this  angle.  Section  5.2 
also  indicated  that  angle  independent  crosstalk  due  to  susceptibility 
saturation  is  expected.  The  data  does  not  extend  to  large  enough  angles 
to  verify  this  expectation  completely,  but  the  curves  in  Figure  5-6  show 
signs  of  flattening  at  about  10  mrad. 

Angle  independent  effects  can  be  inferred  from  another  viewpoint  in 
Figure  5-8,  where  the  noise -conjugate  intensity  and  the  crosstalk  are 
both  plotted  as  percentages  of  their  maximum  values  versus  the  pump- 
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Figure  5-5.  Intensity  crosstalk  results  when  probe  and  noise 
powers  are  equal. 
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Figure  5- 


.  Intensity  crosstalk  data  plotted  in  a  function  of 
noise-pump  angle. 
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5-8.  Noise  conjugate  intensity  and  observed  ''.rosstalk  as  a 

function  of  noise-pump  angle  for  noise/probe  ratio  of  3.0 


noise  angle.  Here  the  noise -probe  Intensity  ratio  is  3:1,  and  we  see  the 
noise-conjugate  intensity  decreasing  nuch  more  rapidly  than  the  cross¬ 
talk.  This  should  be  coopered  with  the  behavior  seen  in  Figure  5-9, 
where  the  same  quantities  are  plotted  for  a  noise-probe  intensity  ratio  of 
1:1.  Here  the  noise -conjugate  intensity  falls  off  only  slightly  faster 
than  the  crosstalk.  The  relative  behavior  in  these  two  figures  can  be 
explained  by  assuming  that  the  crosstalk  is  dominated  at  smaller  angles  by 
pump  depletion,  and  at  larger  angles,  where  the  noise -conjugate  Intensity 
is  small,  by  susceptibility  saturation.  The  greater  large-angle  crosstalk 
in  Figure  5-6  is  expected  because  the  higher  noise  intensity  produced  more 
saturation. 

5.5  POLARIZATION  CROSSTALK  EXPERIMENTAL  DESIGN 

Intensity  crosstalk  studies  were  carried  out  where  the  polarization 
states  of  the  probe  and  noise  beams  were  varied.  The  experimental 
arrangement  diagrammed  in  Figure  5-10  differs  from  the  first  crosstalk 
experiment  (Figure  5-1)  by  the  placement  of  the  noise  and  probe  beams  on 
opposite  sides  of  the  pump  beam  axis.  The  data  collected  here,  when  all 
beams  beams  were  copolarized,  could  then  be  compared  with  early  crosstalk 
results  to  determine  whether  noise-probe  beam  angle  and/or  pump-noise/pump 
probe  beam  angles  were  significant  parameters  in  the  determination  of  the 
extent  of  crosstalk. 

Both  the  probe  and  noise  beams  were  focused  at  the  Na  cell  by  a 
single  high  quality  achromat  L2  (focal  length  45  cm)  and  their  overlap 
with  the  pump  beams  carefully  checked  by  replacement  of  the  Na  cell  with  a 
precision  pinhole.  The  probe-pump  angle  tp  and  the  noise-pump  angle  6  were 
kept  to  the  minimum  allowed  by  finite  beam  diameters  and  were  measured  to 
be  8.5  mrad  and  10.5  mrad,  respectively.  The  linear  polarization  of  the 
noise  and  probe  beams  was  manipulated  by  half  waveplate  polarization 
rotators  WP  placed  behind  L2.  The  waveplates  were  deliberately  arranged 
to  be  the  last  optical  elemnt  in  the  noise  and  probe  beam  paths  in  order 
to  maintain  the  direction  of  linear  polarization  of  both  beams.  Polariza¬ 
tion  directions  were  determined  with  a  Karl  Lambrecht  laser  beam  polarizer 
(MSCAIO)  and  the  degree  of  polarization,  defined  by 

^  "  ^parallel/^^perpendicular'*‘^parallel)  > 
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Noise  conjugate  intensity  and  observed  crosstalk 
as  function  of  noise-pump  angle  for  noise/probe 
ratio  of  0.8. 


was  measured  to  be  greater  than  0.98  for  all  beams  at  the  Na  cell. 

Relative  probe  and  noise  beam  Intensities  were  adjusted  by  variable 
beamsplitter  VBS2  and  the  total  Input  Intensity 

^total  “  ^nolse  ^probe 

was  controlled  by  Insertion  of  neutral  density  filters  NDF  Into  the  beam 
generated  by  variable  beamsplitter  VBSl.  As  before,  both  probe  and  noise 
conjugates  were  monitored  with  matched  photodiodes  PDl  and  FD2  and  their 
outputs  displayed  on  an  oscilloscope  triggered  by  the  noise  beam  chopper 
CH. 

5.6  POLARIZATION  CROSSTALK  RESULTS 

Intensity  crosstalk  was  measured  In  three  distinct  probe -pump -noise 
polarization  states.  The  copolarized  state,  where  all  planes  of  polariza¬ 
tion  were  vertical  and  parallel,  served  as  a  reference  state  for  Intensity 
crosstalk.  The  second  state  Is  referred  to  as  the  parallel  state  and  Is 
established  when  the  plane  of  polarization  of  the  probe  and  noise  beams 
are  parallel,  but  both  are  angled  at  45*  with  respect  to  the  pump  polariza¬ 
tion  plane.  The  third  state,  wh^re  the  probe  and  noise  polarization 
planes  are  orthogonal  to  each  other  but  at  45*  with  respect  to  the  pump 
plane.  Is  referred  to  as  the  perpendicular  state.  The  results  of  two 
Independent  experiments  are  presented  In  Figure  5-11  where  the  crosstalk 


obtained  In  the  parallel  and  perpendicular  states  Iti  compared  to  that  of 
the  reference  copolarized  state.  Although  the  measured  crosstalk  for  the 
copolarized  state  differed  slightly  In  the  two  experlmencs.  In  both  cases 
the  parallel  and  perpendicular  states  yielded  essentially  the  same 
percentage  crosstalk  values  as  the  reference  state.  Hence,  Intensity 
crosstalk  does  not  appear  to  depend  on  the  orientations  of  the  pump, 


probe,  and  noise  beam  polarization  planes. 
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Figure 

5-11.  Polarization  Crosstalk 

Results 

Case 

Na  Cell 

Polarization 

Pump 

Probe 

Noise 

Percent 

Temperature 

State 

Power 

Power 

Power 

Crosstalk 

CC) 

copolarized 

(mW) 

(mW) 

(mW) 

31.8 

1 

250 

perpendicular 

400 

9.9 

8.4 

29.9 

copolarized 

26.9 

2 

248 

perpendicular 

380 

8.5 

7.6 

26.7 

*Average 

of  six  measurements 

5.7  CROSSTALK  MODEL  VERIFICATION 

If  pump  depletion  and  susceptibility  saturation  are  the  principal 
causes  of  noise -Induced  crosstalk,  then  in  a  noise -free  interaction  they 
should  manifest  themselves  as  saturation  of  the  probe -conjugate  intensity 
at  higji  probe  intensities.  The  experimental  approach  used  to  test  this 
assumption  is  outlined  in  Figure  5-12.  Two  counterpropagating  pump  beams 
of  equal  powers  were  formed  by  BSl  and  focused  to  a  400  /xm  1/e^  diameter 
at  the  Na  cell.  Pump  intensities  were  typically  125-132  W/cm^.  A  single 
probe  beam  was  directed  into  the  cell  by  BS2  and  L2  to  form  a  200  ^m  1/e^ 
diameter  spot  at  the  cell.  The  pump-probe  angle  was  8.5  mrad.  Probe 
intensity  was  varied  from  2-40  W/cm^  by  insertion  of  neutral  density 
filters  NDF  into  the  probe  beam  path.  Care  was  taken  to  precisely  locate 
the  neutral  density  filters  to  insure  constant  pump-probe  overlap  at  the 
cell.  At  each  probe  input  intensity,  reflected  power  and  gain  were 
measured  by  PDl  and  PD2  respectively.  The  results  of  two  independent 
trials  are  listed  in  Figures  5-13  and  5-14  and  displayed  in  Figure  5-15, 
where  reflected  power  as  a  percentage  of  pump  power  is  plotted  against 
input  probe  intensity  for  both  sets  of  data;  the  two  trials  are  distin¬ 
guished  by  circles  or  triangles  around  the  plotted  points.  It  is  evident 
that  the  reflected  power  saturates  as  the  probe  intensity  is  increased; 
this  will  be  shown  to  account  for  the  observed  intensity  crosstalk. 
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The  experimental  points  in  Figure  5-15  can  be  fit  very  well  with  the 
expression, 

Ipc  Alp 

Ipump  ^  ®Ip 

where  Ipump  is  the  total  pump  intensity,  Ip  and  Ip^,  are  the  probe  and 
probe -conjugate  intensities,  all  in  W/cm^,  and  A  and  B  are 
A  -  1.71  X  10-3  cm^/W, 

B  -  7.27  X  10*2  cm2/W 


(5-3b) 


The  points  inside  the  squares  in  Figure  5-15  are  calculated  from  Eqn 
5-3.  It  should  be  remembered  that  the  values  of  A  and  B  apply  to  a 
particular  pump-probe  angle  (8.5  mrad)  and  pump  Intensity  (125-132  W/cm2) , 
and  presumably  are  dependent  upon  these  variables. 

It  has  been  shown  that  numerical,  theoretical  predictions  of  probe- 
conjugate  intensity  in  a  Kerr-like  medivtm  can  be  fit  quite  accurately  with 
an  expression  of  the  form  (Ref  5-4) , 

^pc 

^pump  ^  +  alp 

Since  B/A  »  1  in  Eqn.  5-3,  it  follows  that  the  curve  in  Figure  5-15 
saturates  much  more  strongly  than  that  predicted  for  a  Kerr- like  medium, 
although  it  must  be  remembered  that  Eqn.  5-3  has  not  been  verified  beyond 
a  probe  intensity  of  about  40  W/cm2.  Nevertheless,  the  comparison  shows 
that  in  Na,  if  the  theory  is  correct,  there  must  be  an  additional  saturat¬ 
ing  mechanism  at  least  as  significant  as  pump  depletion.  The  most  likely 
candidate  is  coupling  saturation,  as  suggested  by  the  experimental  results 
on  crosstalk  vs.  noise  input  angle,  and  crosstalk  vs.  probernoise  intensi¬ 
ty  ratio. 


With  the  aid  of  two  assumptions,  Eqn.  5-3  can  be  used  to  predict 
the  crosstalk  for  an  experiment  in  which  the  noise  and  probe  beams  make 
the  same  angle  with  the  pump  beam  as  the  probe  beam  did  in  the  previous 
experiment,  and  for  which  the  pump  intensity  is  also  the  same  as  before. 
The  first  assumption  is  that  the  noise  and  probe  together  can  be  regarded 
as  a  single  beam,  and  that  the  noise-  and  probe -conjugates  can  similarly 
be  regarded  as  a  single  beam.  That  is. 
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lyyr. 


PROBE  INTENSITY  (W/cm^) 


Fig.  5-15. 
Reflected  power  as 


Probe  conjugate  power  as  function  of  input  intensity. 

%  of  pump  power  saturates  at  higher  power  input  intensitie 


+ 


^pc  ^nc  A(Ip  +  Ijj) 

Ipump  ^  ®^^p 


AI, 


A1 


n 


1  +  B(Ip  +  In) 


1  +  B(Ip  +  In)  • 


The  second  assumption  is  that  the  ratio  of  probe -conjugate  to  noise - 
conjugate  Intensities  is  equal  to  the  ratio  of  probe  to  noise  intensities, 


'•pc 


AI, 


Ipump  ^  ®^^p  ^n^ 

Defining  the  crosstalk  C  by 

Ipp(no  noise)  -  Ipn(noise) 

^  B  ■  - 

Ipp(no  noise) 

and  using  Eqn.  5-4,  we  obtain 


(5-4) 


C 


BI 


n 


1  +  B(Ip  +  In) 


(5-5) 


Figure  5-16  is  a  graphical  explanation  of  Eqn.  5-5.  Probe  and 
noise  intensities  are  summed  on  the  x-axis  to  determine  the  total  conju¬ 
gate  power  expected.  This  is  divided  into  noise  and  probe  components  in 
the  ratio  of  the  noise  and  probe  input  intensities.  For  equal  noise  and 
probe  Intensities,  the  probe  conjugate  becomes  Cp+jj/2  as  shown  in  the 
bottom  portion  of  Figure  5-16,  and  that  is  less  than  the  intensity  Cp  that 
would  have  been  obtained  in  the  absence  of  a  noise  beam. 

A  test  of  the  predictive  capability  of  Eqn.  5-5  was  carried  out 
with  the  experimental  arrangement  of  Figure  5-10.  A  noise; probe  inten¬ 
sity  ratio  of  1.16:1.00  was  used  with  a  pump-probe  angle  of  8.5  mrad  and  a 
pump-noise  angle  of  10.5  mrad;  one  can  see  from  Figures  5-8  and  5-9  that 
there  is  a  relatively  small  difference  in  noise  beam  reflectivity  at 
these  two  angles,  so  that  Eqn.  5-5  should  be  applicable.  Total  input 
inteiisity  (Ipj-  +  I^)  tanged  from  1.5-28  W/cm^  and  was  varied  by  the 
insertion  of  neutral  density  filters  into  the  beam  which  generated  both 
the  probe  and  noise  beams.  The  effect  of  the  insertion  of  each  neutral 
density  filter  on  the  spatial  overlap  of  pump,  noise,  and  probe  beams  was 
carefully  checked  in  order  maintain  a  constant  overlap  in  the  Na  cell. 

Experimental  conditions  were  otherwise  Identical  to  those  of  Figure  5-15. 

The  results  of  this  intensity  crosstalk  experiment  are  shown  in 
Figure  5-17  together  with  the  predicted  crosstalk  values  derived  from 
Eqn.  5-5.  There  Is  good  agreement  between  predicted  and  measured  values 


TOTAL  INPUT  INTENSITY  (Ip  +  In) 


NOTE;  Cp>n  EQUALLY  OISRIBUTED  TO  NOISE  AND  PROBE  CONJUGATES 


Diagramatic  representation  of  intensity  crosstalk 
mechanism.  It  is  the  nonlinear  dependence  of  conjugate 
reflectivity  on  input  probe/noise  intensity  that  leads 
to  crosstalk. 


NOISE:  PROBE  -  1.16:1.00  (200  /i  l/a^  DMMETER):  PUMP-PROBE  ANGLE  -  S.1-9.3  nrid 
PUMP  INTENSITY  -  132  W/cm^  (400  /i  l/a^  DIAMETER):  PUMP-NOISE  ANGLE  •  10.S  Mnd 


Predicted  and  observed  intensity  crosstalk  measurements 


at  all  but  the  highest  input  intensity.  The  agreement  at  the  lower 
intensities  indicates  that  the  same  mechanisms  that  cause  saturation  of 
the  probe -conjugate  intensity  in  the  noise- free  experiment  -pump  depletion 
and  coupling  saturation  being  the  prime  candidates-  are  principally 
responsible  for  crosstalk.  Therefore,  at  the  smaller  input  intensities, 
those  mechanisms  discussed  in  Section  5.2  that  are  operative  only  in  the 
presence  of  noise  are  not  significant  in  causing  crosstalk.  In  partic¬ 
ular,  this  excludes  direct  probe  and  probe -conjugate  depletion  as  impor¬ 
tant  sources  of  crosstalk  at  the  lower  noise -plus -probe  intensities. 

At  the  high  probe-plus-noise  Intensity  of  28  W/cm^,  however,  the 
sltviation  is  somewhat  different.  Here  the  measured  crosstalk  is  signifi¬ 
cantly  greater  than  that  predicted,  and  the  latter  is  well  outside  the 
experimental  error  bars.  A  possible  explanation  is  that,  at  sufficiently 
high  probe  and  noise  intensities,  additional  coupling  mechanisms  such  as 
that  described  by  the  triplet  E^EpE^  come  into  play. 

An  indication  of  the  effect  of  increased  probe  intensity  on  extracta¬ 
ble  pump  power  can  be  seen  when  the  net  extracted  power  (NEP)  defined  as : 

NEP  -  (transmitted  power  +  reflected  power)  -  input  power 
is  plotted  against  input  intensity  as  in  Figure  5-18.  For  probe  intensi¬ 
ties  up  to  approximately  15  W/cm^,  a  steady  increase  in  extracted  pump 
power  is  recorded.  However,  with  greater  intensities  the  extracted  pump 
power  decreases  until  at  AO  W/cm^  the  NEP  is  approximately  the  value 
obtained  at  2  W/cm^.  One  explanation  for  the  observed  roll-over  is  that  at 
probe  Intensities  >  15  W/cm^,  the  probe  is  strong  enough  to  act  as  a 
"pump"  beam  and  scatter  power  back  into  the  true  pump  beams,  thus  decreas¬ 
ing  its  intensity  and  its  ability  to  extract  power  from  the  true  pumps. 

5.8  TOLERABLE  NOISE  INTENSITY;  An  Example 

To  illustrate  the  potential  practicability  of  our  crosstalk  results, 
we  use  them  to  estimate  the  tolerable  noise  intensity  in  an  intensity 
modulated  communication  system.  We  assume  that  a  constant  intensity  probe 
beam  is  sent  through  an  aberrating  medium  to  a  FWM,  where  a  time -varying 
intensity  signal  is  Imposed  on  the  probe -conjugate  beam;  the  latter  than 
returns  through  the  aberrating  medium  to  a  detector.  The  information  is 
encoded  in  binary  form;  the  probe -conjugate  Intensity  Ip,,  is  either 
greater  than  or  equal  to  a  preset  switch  intensity  I^^  (Ip,,  on,  binary 


Figure  5-18.  Net  extracted  power  as  a  function  of  probe  input  intensity,  indicating 
the  possible  contributions  of  other  FWM  processes  when  the  probe  input 
field  is  strong. 


one),  or  less  than  Ig^  (Ipc  binary  zero).  For  a  valid  "on"  signal  to 

be  received,  crosstalk  noise  in  the  FUM  must  not  reduce  Ip^  below  Ig^.  To 
examine  the  noise  problem,  we  start  with  Eqn.  S-4.  Remember  that  this 
equation  applies  to  a  specific  situation  where  the  total  pump  intensity  is 
132  U/cm^,  the  noise  beam  angle  with  the  pvimp  beam  is  8.5  mrad,  and  the 
noise  frequency  is  degenerate  with  the  pump  frequency.  This  estimate, 
therefore,  is  not  general  but,  because  of  the  small  angle  and  frequency 
range  retrlctions,  it  may  be  regarded  as  a  worst  case  calculation. 

It  is  convenient  to  rewrite  Eqn.  5-4  in  the  form. 


Ipc  “  ^a 


Ip/Ib 

1  +  (Ip  +  In)/Ib  • 


where 


Ig  -  (A/B)  Ipump  -  31.0  W/cm2,  and 
1^  -  1/B  -  13.8  W/cm2. 

From  this  expression,  we  solve  for  Ipt 

^pc^^n  ^b^ 

I - ^ , 

J-a  ■'■pc 


(5-6) 


Ue  now  ask  how  much  the  probe  intensity  must  be  increased  to  maintain 
an  "on"  signal  in  the  presence  of  crosstalk  noise.  Without  noise,  the 
minimum  probe  intensity  for  an  "on"  signal  is,  from  Eqn.  5-6, 


lp(0)  “  I  .  I  ’ 

■‘■a  ^sw 

and  in  the  presence  of  noise  it  is 


Ipdn)  - 

The  ratio  is  then 


^sw^^n  ^b^ 
^a  '  ^sw 


Ipdn) 

lp(0) 


-  1  -f 


Ipdn) 

lp(0) 


-  1  + 


In(W/cm2) 


Ue  see  that  a  noise  intensity  as  large  as  1.38  V/cm^  would  require  a 
relatively  insignificant  10%  increase  in  the  probe  intensity  to  overcome 


the  intensity  crosstalk.  In  addition,  is  tightly  restricted  to  be 
within  about  8.5  mrad  of  the  pump  beam  and  degenerate  with  the-FWM  pump 
frequency,  so  the  overall  result  is  that  intensity  crosstalk  will  have 
little  impact  on  field  applications. 

5.9  INTENSITY  CROSSTALK  SUMMARY 

Simultaneous  measurements  of  intensity  crosstalk  and  noise-conjugate 
intensity  versus  the  noise-pump  angle  indicate  two  principal  crosstalk 
mechanisms,  one  having  a  sharp  angular  dependence  and  vanishing  for  noise - 
pump  angles  greater  than  about  10  mrad,  the  other  being  independent  of 
angle.  A  strong  correlation  between  the  sharp,  small-angle  fall-offs  of 
the  crosstalk  and  the  noise-conjugate  intensity  indicates  that  the  angle- 
dependent  crosstalk  source  is  depletion  of  the  pump  beams  to  provide  power 
for  the  generation  of  the  noise -conjugate  beam.  The  angle -independent 
source  appears  to  be  intensity  saturation  of  the  nonlinear  susceptibility, 
that  is,  of  the  coupling  between  the  pumps  and  the  probe  and  probe - 
conjugate  beams. 

Experiments  to  determine  the  sensitivity  of  intensity  crosstalk  to 
the  relative  polarization  directions  of  the  pump,  probe,  and  noise  waves 
showed  no  discemable  effect. 

The  probe • conj ugate  intensity  was  measured  as  a  function  of  input 
probe  intensity  for  fixed  pump  intensity.  The  curve  was  found  to  increase 
linearly  for  small  probe  intensity,  and  then  to  saturate  for  larger  probe 
intensity.  This  behavior  can  also  be  explained  as  being  due  to  a  combina¬ 
tion  of  pump  depletion  and  intensity  saturation  of  the  nonlinear  coupling. 
With  two  assumptions ,  this  data  was  used  to  predict  the  results  of  a 
measurement  of  crosstalk  vs.  combined  probe  and  noise  intensities.  The 
first  assumption  was  that  the  probe  and  noise  beams  can  be  treated  as 
parts  of  the  same  beam,  l.e.,  that  in  the  presence  of  noise,  the  probe- 
conjugate  intensity  bears  the  same  relationship  to  the  total  (probe  plus 
noise)  input  intensity  that  it  does,  in  the  absence  of  noise,  to  the  input 
probe  intensity  alone.  The  second  assumption  was  that  the  probe-  and 
noise -conj ugate  beams  can  also  be  treated  as  parts  of  the  same  beam,  i.e., 
that  the  ratio  of  probe -conj ugate  to  noise -conj ugate  intensities  is  equal 
to  the  ratio  of  input  probe  to  input  noise  intensities.  The  prediction 
and  the  experimental  results  agreed  well,  demonstrating  that  the  same 


103 


mechanisms  responsible  for  saturation  of  the  probe - conj ugate  Intensity  vs. 
probe  Intensity  curve  are  also  responsible  for  Intensity  crosstalk.  This 
agreement  lends  strong  support  to  the  hypothesis  that  crosstalk  results 
principally  from  pump  depletion  and  Intensity  saturation  of  the  nonlinear 
coupling. 
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6.0  SPATIAL  CROSSTALK  AND  IMAGE  FIDELITY 


Decrlbed  in  this  section  are  experiments  to  establish  the  Importance 
of  another  type  of  crosstalk  involving  transfer  of  spatial  information 
from  one  beam  to  another.  Such  crosstalk  would  be  deleterious  in 
applications  where  image  fidelity  was  important.  The  results  presented 
here  Indicate  that  spatial  crosstalk  arises  only  through  the  intensity 
crosstalk  mechanism  described  in  Section  5,  i.e.,  as  predicted  there  is  no 
transfer  of  phase  information  from  a  noise  beam  to  a  signal  beam  if  their 
intensities  are  both  small  compared  to  the  pump  beams.  The  imaging  FWM 
experiments  reported  here  also  demonstrate  excellent  image  fidelity.  We 
report  a  resolution  of  better  than  14.3  line  pairs  per  millimeter,  the 
highest  line  pair  resolution  obtained  in  Na  vapor,  and  have  demonstrated 
faithful  correction  of  a  severe  optical  aberration. 

6 . 1  Experimental  Description 

To  establish  the  effects  of  spatial  crosstalk  on  the  contrast  and 
resolution  of  the  probe  conjugate  image,  initial  imaging  efforts 
concentrated  on  phase  conjugation  of  high  frequency  spatial  images.  The 
arrangement  in  Figure  6-1  produced  conjugate  images  with  up  to  14.3  line 
pairs  per  millimeter  resolution,  limited  to  this  figure  by  camera  pixel 
resolution  and  the  available  cw  laser  power.  A  90.0  cm  focal  length 
achromat  (LI)  was  placed  in  the  pump  be£uii  path  to  produce  an  1100  fim 
1/e  diameter  pump  cross  section  at  the  Na  cell.  Spatial  information 
was  imprinted  on  the  probe  beam  by  introducing  a  standard  USAF  resolution 
target  (RT)  into  the  probe  beam  path  upstream  of  BS2.  The  spatially 
Imprinted  probe  beam  was  then  focused  into  the  Na  cell  at  a  pump -probe 
angle  of  25.3  mrad  by  a  20.0  cm  focal  length  achromat  (L2) .  The  conjugate 
image  retraced  the  probe  beam  path  (dotted  line)  and  was  reflected  from 
BS2  into  the  camera  (C) .  Care  was  taken  to  optimize  the  positions  of  L2 , 
RT,  and  C  to  produce  the  highest  resolution  images. 

6.2  Image  Fidelity  Results 

Figure  6 -2a  is  a  photograph  of  a  video  monitor  which  displayed  a 
portion  of  the  resolution  target  used  in  the  image  resolution  experiment. 


The  image  for  Figure  6-2a  was  transferred  from  the  resolution  target  to 
the  camera  detector  plane  by  1:1  Imaging  optics  to  test  the  spatial 
frequency  resolution  limit  Imposed  by  the  camera  pixel  spacing.  The 
resolution  portrayed  in  the  photograph  is  unavoidably  degraded  compared  to 
the  actual  resolution  obtained  due  to  the  display  resolution  of  the  video 
monitor.  However,  video  recording  (as  opposed  to  direct  photographic 
recording)  proved  to  be  the  most  efficient  and  convenient  means  of 
characterizing  the  imaging  results  because  of  its  real  time  processing 
capabilities  in  the  determination  of  focal  plane  location,  intensity 
levels,  and  intensity  xmlformity.  Figure  6-2b  is  a  photograph  of  the  same 
video  monitor  displaying  the  highest  resolution  conjugate  image  that  was 
recorded.  The  resolution  of  the  two  photographs  is  seen  to  be  comparable 
indicating  that  pixel  separation  in  the  camera  detector  plane  probably  was 
a  limiting  factor  in  the  recording  of  images  with  higher  spatial 
frequencies . 

A  second  factor  which  affected  spatial  frequency  resolution  was  the 
effective  aperture  of  the  phase -conjugate  mirror  determined  by  the  pump 
beam  cross  section  at  the  Na  cell.  This  diameter  in  turn  was  determined 
by  the  available  laser  power  from  the  CW  dye  laser,  with  a  compromise 
between  pump  intensity  for  reflectivity  and  pxuap  diameter  for  image 
fidelity  being  necessary.  Employing  the  shortest  focal  length  achromat 
practical  for  L2  (focal  length  <-20.0  cm),  LI  was  translated  along  the 
common  pump  beam  path  to  produce  the  highest  resolution  images  possible 
while  maintaining  acceptable  contrast  levels.  The  corresponding  pump  beam 
l/e  diameter  was  1100  pm.  The  effect  of  this  aperture  size  was 
assessed  by  focussing  the  image  bearing  probe  beam  through  a  1100  pm 
hard  aperture,  and  comparing  the  transmitted  image  with  the  conjugate 
image  reflected  off  of  the  sodium  phase  conjugate  mirror.  The 
phase -conjugate  image  reflected  off  of  the  1100  pm  "mirror"  formed  by 
the  pump  beams  was  noticeably  better  in  fidelity  than  the  image 
transmitted  through  the  hard  aperture.  This  suggests  that  there  is  enough 
intensity  in  the  wings  of  the  pump  beams  to  conjugate  the  image  beam 
efficiently,  i.e.,  the  phase -conjugate  aperture  is  a  soft  one  and  larger 
than  the  l/e  pump  beam  diameter  at  the  pump  powers  used  in  these 
experiments . 

Typical  reflectivities  obtained  with  1100  pm  pump  beams  in  the 


imaging  experiments  were  12-18%,  with  excellent  conjugate  fidelity 
evidenced  by  the  necessity  of  placing  the  resolution  target  (RT)  and 
camera  (C)  equidistant  from  BS2  (Fi^.  6-1).  If  these  distances  were  not 
the  case,  the  ceunera  was  not  at  the  conjugate  image  plane,  resulting  in 
the  conjugate  return  image  being  detected  out  of  focus.  The  excellent 
fidelity  was  also  established  by  testing  image  quality  with  an  optical 
aberrator  in  the  probe  path. 

The  aberration  correction  capability  of  the  four  wave  mixing  cell  was 
demonstrated  with  the  apparatus  in  Fig.  6-3.  Initially,  a  standard  USAF 
target  image  (in  this  case,  the  numeral  4)  was  imprinted  on  the  probe  beam 
by  the  image  plate  (IP)  and  focused  into  the  Na  cell  by  L2.  The  conjugate 
image-bearing  beam  (dotted  line)  retraced  the  input  probe  path  and  was 
reflected  off  the  beamsplitter  to  the  camera  (C) .  A  second  camera  was 
employed  to  monitor  the  quality  of  the  transmitted  image  which  was 
reflected  from  two  A/20  turning  mirrors  and  focused  onto  the  detector 
plane  of  the  transmission  camera  by  achromat  L3. 

Photographs  of  the  unaberrated  phase  conjugate  image  and  transmitted 
image  are  shown  in  Figures  6 -4a  and  6 -4b,  respectively.  Although  the 
transmissive  aperture  of  the  four  wave  mixing  cell  is  significantly  larger 
than  the  phase  conjugating  aperture  (note  the  edges  of  the  bars  are  absent 
in  the  conjugate  image  photograph),  the  central  portion  of  the  transmitted 
image  served  as  a  good  input  image  standard.  After  recording  the 
unaberrated  Images ,  a  65  A  aberrator  plate  (A)  was  inserted  into  the 
input  probe  image  path  so  that  the  distorted  phase  conjugated  image 
returned  through  the  aberrator  plate  before  reaching  the  camera.  The 
resulting  photographs  of  the  aberrated  conjugate  image  (Fig.  6-4c)  and 
corresponding  transmitted  image  (Fig.  6-4d)  are  displayed  below  their 
unaberrated  counterparts.  The  ability  of  the  four  wave  mixing  cell  to 
correct  for  significant  image  distortion  is  then  made  apparent  by 
comparison  of  Fig.  6-4a,  6-4c,  and  6-4d.  The  open  gap  at  the  top  of  the 
target  "4"  is  conservatively  estimated  to  represent  a  10  line  pair  per 
millimeter  spacing;  therefore,  the  present  system  is,  at  a  minimum, 
capable  of  this  degree  of  resolution  through  a  65  wave  aberration. 

6.3  Spatial  Crosstalk  Results 

After  good  image  resolution  was  achieved  with  fairly  high  spatial 
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Figure  6-4.  Aberration  correction  experimental  results,  a)  phase  conjugate 
image  with  no  aberrator  in  propagation  path,  b)  transmitted  probe 
image  with  no  aberrator,  c)  phase  conjugate  image  with  65  X  aberrator 
in  place,  d)  transmitted  probe  image  with  aberrator  in  place. 


frequency  targets,  spatial  crosstalk  experiments  were  undertaken  to 
quantify  reductions  in  image  contrast  and  edge  definition  in  the  presence 
of  noise  beams  that  have  spatial  intensity  distributions  that  .are  Gaussian 
or  image-bearing  .  For  a  Gaussian  noise  source,  modifications  to  the 
basic  imaging  arrangement  of  Fig.  6-1  Included  the  addition  of  a  variable 
beamsplitter  (VBS2,  Fig.  6-5)  and  a  50%  beamsplitter  (BS3)  to  produce  a 
variable  intensity  noise  beam  which  was  separately  focused  into  the  Na 
cell  at  a  pump-noise  angle  of  8.1  mrad  by  achromat  L3.  The  noise  beam 
conjugate  (dotted  line)  retraced  the  noise  beam  path  and  was  detected  by  a 
silicon  photodiode  (PD)  through  BS3.  Maximization  of  the  noise  conjugate 
by  adjustment  of  elements  L3  and  BS3  insured  optimum  probe-noise  overlap 
at  the  Na  cell  and  was  essential  for  consistent  spatial  crosstalk 
results.  The  total  power  allocated  to  the  probe  and  noise  beams  by  VBSl 
could  not  exceed  6-10  mW  without  sacrificing  pump  beam  power. 

Quantification  of  the  changes  in  probe  conjugate  image  quality  due  to 
the  presence  of  Gaussian  or  image-bearing  noise  fields  was  accomplished 
through  the  use  of  video  digitization  of  the  recorded  conjugate 
images. Fig.  6 -6a  and  6 -6b  are  representative  of  the  digital  output  from 
two  cross  section  samples  of  a  2.83  line  pair  per  millimeter  three  bar 
image  recorded  in  the  absence  of  the  noise  beam.  Due  to  the  comparable 
dimensions  of  the  test  target  and  probe  beam  diameter,  the  center  bar  of 
the  three  bar  conjugate  image  was  somewhat  more  intense  than  the  outer  two 
bars.  However,  good  contrast  and  distinct  boundries  are  apparent  from  the 
light- to-dark  values  and  well-defined  light  and  dark  regions  are  obtained 
with  a  4  mW  input  probe  image.  In  the  presence  of  a  2  mW  Gaussian  noise 
source  at  approximately  one -third  the  probe  beam  input  angle  (8.1  mrad  for 
the  Gaussian  noise  beam  versus  25.3  mrad  for  the  image-bearing  probe 
beam),  the  same  cross  section  samples  (Fig.  6-7a  and  6-7b)  show  only  a  15% 
average  decrease  in  peak  intensity  with  no  detectable  degradation  in  edge 
sharpness . 

Since  even  at  small  input  angles,  a  Gaussian  cross  section  noise  beam 
had  little  effect  on  probe  conjugate  image  quality,  the  effects  of  an 
image-bearing  noise  beam  were  studied  with  the  experimental  layout  in  Fig. 
6-8.  A  three  bar  test  target  of  identical  spatial  frequency  but 
orthogonal  orientation  with  respect  to  the  probe  target  (see  Fig.  6-8)  was 
then  Inserted  into  the  noise  beam  path  and  focused  into  the  Na  cell  with 


achromaC  L3.  All  optical  elements  in  the  probe  and  noise  beam  paths  were 
identical  so  that  the  variable  beamsplitter  (VBS2)  directly  controlled  the 
probeinolse  intensity  ratio  at  the  cell.  In  addition,  a  second  camera  was 
positioned  to  record  the  quality  and  intensity  of  the  noise  beam  conjugate 
image  in  order  to  insure  that  both  probe  and  noise  beams  were  conjugated 
from  the  same  region  in  the  conjugate  reflecting  aperture.  Each  of  the 
input  beams  was  focused  separately  into  the  cell  to  produce  the  most 
homogeneously  illuminated  conjugate  image  at  their  respective  cameras  and 
the  images  were  recorded.  Both  beams  were  then  unblocked  and  the  probe 
conjugate  image  recorded.  The  result  of  this  sequence  where  the 
probe: noise  Intensity  ratio  was  approximately  1:6  is  displayed  in  the 
photographs  of  Fig.  6-9a-c  and  in  the  digitized  cross  sections  of  Fig. 

6 -10a  and  6 -10b.  The  average  peak  intensity  is  reduced  approximately  46% 
as  a  result  of  the  presence  of  the  noise  image;  however,  as  in  the  case  of 
a  Gaussian  noise  field,  edge  definition  is  not  compromised  even  by  a 
six- fold  intense  noise  image. 

The  spatial  crosstalk  data  can  be  understood  in  terms  of  the 
intensity  crosstalk  mechanism  described  in  Section  5  if  one  considers  the 
spatial  frequency  content  of  the  Fourier  transform  of  the  probe  and  noise 
images.  Since  the  interaction  occurs  at  or  near  the  focal  planes  of  the 
lenses  that  focus  the  probe  and  noise  beams  into  the  sodium  cell,  the 
interfering  spatial  Images  are  the  Fourier  transforms  of  the  images 
originally  imprinted  on  the  beams.  Fig.  6-lla  and  6-llb  are  photographs 
of  the  Fourier  transforms  of  the  probe  and  noise  images  with  an  1100  /xm 
diameter  aperture  for  comparison  which  simulates  the  1/e^  pump  beam 
diameter  at  the  Na  cell.  Intensity  crosstalk  can  occur  in  regions  of  the 
Fourier  plane  where  there  is  significant  noise  beam  intensity,  which 
because  of  its  localized  distribution,  results  in  spatially  varying 
crosstalk  effect. 

This  localized  spatial  crosstalk  is  strongest  precisely  where  the  low 
spatial  frequency  (e.g.  contrast)  components  of  the  probe  image  are 
located  in  the  sodium  FWM  cell.  Most  of  the  power  in  both  the  probe  and 
noise  beams  is  concentrated  inside  the  first  Airey  disk,  which  in 
transform  space  contains  the  low  spatial  frequency  components  of  the 
oi^isinal  patterns.  As  a  result,  the  intensity  of  these  components  is 
reduced,  leading  to  a  loss  of  image  contrast.  Higher  frequency  Fourier 
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Figure  6-10(b).  Digitized  cross  sections  of  probe  conjugate 
Image  with  noise  Image. 


components,  which  define  the  edge  boundaries,  occur  in  the  outer  wings  of 
the  patterns  in  Fig.  6-11  and  carry  only  a  small  fraction  of  the  intensity 
of  the  original  image.  Therefore,  the  higher  order  components  for  both 
the  probe  and  noise  images  are  much  less  likely  to  be  involved  with  pump 
beam  depletion  and  coupling  saturation,  and  are  conjugated  independently 
of  one  another.  Thus  edge  definition  in  the  conjugate  image  is  not 
expected  to  be  impacted,  a  result  consistent  with  the  observed  "dimming" 
without  loss  of  edge  definition  of  the  conjugate  image  in  the  far  field. 
The  preservation  or  effective  enhancement  of  edges  is  important  for  target 
image  recognition. 

6 . 4  Summary 

Experimental  results  have  shown  that  excellent  image  quality  can  be 
achieved  using  a  Na  vapor  FWM,  even  when  significant  optical  aberrations 
are  present  in  the  propagation  path.  The  resolution  obtained  was 
essentially  instrument-limited;  in  practice,  the  resolution  will  be 
determined  by  the  soft  aperture  defined  by  the  pump  beam  diameter.  In  our 
experiments ,  this  diameter  was  chosen  as  a  compromise  between  pump 
Intensity  and  effective  aperture  size,  with  resulting  reflectivities  in 
the  range  of  10-25  %.  Under  these  conditions,  the  impact  of  self-action 
effects  mentioned  in  Section  4  were  negligible;  such  effects  may  adversely 
affect  image  quality  under  high  reflectivity  conditions. 

Spatial  crosstalk  between  a  noise  and  probe  beam  was  measured  when 
images  were  imprinted  on  the  probe  beam  only,  and  when  both  beams  were 
image -bearing.  In  both  cases,  the  result  was  a  relatively  uniform 
reduction  in  intensity  of  the  probe  conjugate  image,  i.e.,  a  decrease  in 
image  contrast,  but  there  was  no  measurable  loss  of  resolution.  The 
spatial  crosstalk  effect  may  be  explained  as  localized  intensity  crosstalk 
due  to  saturation  of  the  FWM  reflectivity  in  regions  where  noise  beam 
intensity  in  the  FWM  cell  is  largest.  This  tends  to  be  in  regions  where 
low  spatial  frequency  components  of  the  probe  image  Fourier  transform  are 
located,  so  the  effect  is  primarily  a  loss  of  image  contrast  with  little 
loss  of  edge  definition. 


7.0  FWM  CODE  DEVELOPMENT 


The  strategy  for  FWM  code  development  was  to  start  with  a  simple 
one -dimensional  model  of  FWM  in  a  Kerr  medium,  which  could  be  checked 
against  simple  analytical  expressions  for  reflectivity,  and  then  to  extend 
this  model  to  resonantly  enhanced  media  and  to  a  two-dimensional  wave 
optics  model.  The  choice  of  models  developed  was  driven  by  the  issues 
important  for  the  missions  being  considered.  At  the  beginning  of  the 
program,  we  were  considering  use  of  pulsed  lasers  for  FWM;  this  indicated 
the  desirability  of  time -dependent  models.  The  use  of  pulsed  lasers  was 
abandoned  when  the  difficulty  of  meeting  timing  requirements  over  unknown 
mission  distances  became  obvious.  However,  the  use  of  "time-dependent" 
numerical  models  was  retained  for  their  computational  advantages  in 
two-dimensional  problems.  Although  the  equations  for  the  one -dimensional 
model  can  be  solved  by  a  variety  of  numerical  techniques ,  such  as  Runge 
Kutta  or  Adams  Predictor  corrector  methods,  these  techniques  are  not 
applicable  to  a  two-dimensional  model  that  includes  transverse 
derivatives.  The  two  major  approaches  for  solving  the  2D  steady-state 
problem  are  iterating  from  an  initial  guess  or  approximation,  or  solving 
the  time  dependent  problem  until  steady- state  is  reached.  Previous 
experience  in  solving  2D  coupled,  nonlinear  equations  by  iteration  has 
shown  this  technique  to  be  exceedingly  troublesome  and  susceptible  to 
chaotic  behavior  (induced  by  the  numerical  procedure,  not  by  the 
physics).  For  this  reason,  the  2D  FWM  code  uses  a  "time -dependent" 
algorithm.  A  ID  precursor  to  the  2D  code  was  also  written,  using  the  same 
"time  dependent"  algorithm.  This  ID  code  allows  us  to  define  the  region 
of  validity  of  the  2D  code  and  to  determine  the  input  parameters  needed  to 
produce  accurate  predictions  of  FWM  reflectivity.  It  was  also  found 
useful  to  have  a  non- time -dependent  version  of  the  ID  code.  This  code 
used  standard  numerical  techniques  to  solve  for  the  steady- state 
reflectivity,  and  in  general,  ran  faster  than  the  time -dependent  code. 

7.1  ONE -DIMENSIONAL  DFWM  MODELS 

Early  theoretical  treatments  of  DFWM  in  resonant  media  had  a  very 
limited  range  of  validity  because  the  effects  of  pump  beam  absorption  and 
depletion  were  neglected  (Ref.  7-1).  However,  in  cases  where  the  largest 
DFWM  signals  are  obtained  there  is  a  significant  amount  of  pump  absorption 
(-50%  in  our  experiments) .  The  inclusion  of  pump  absorption  into  the 
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differential  equations  which  describe  the  evolution  of  the  four  complex 
field  amplitudes  results  In  a  set  of  equations  which  must  be  solved 
numerically.  This  problem  has  been  solved  previously  and  the  numerical 
results  of  certain  sample  cases  displayed  (Ref.  7-2,  7-3).  As  TRW  is 
Interested  In  examining  DFUM  over  a  large  parameter  space,  a  similar 
computer  code  has  been  developed  to  solve  these  coupled  field  equations 
numerically. 

7.1.1.  CW  Degenerate  FWM  In  Resonant  Media 

Two  different  ID  models  of  DFWM  have  been  developed  at  TRW.  The 
version  to  be  discussed  first,  known  as  IDCW,  is  applicable  to  DFWM  in 
homogeneously -broadened  resonant  media  using  a  CW  narrowband  laser.  The 
coupled  field  equations  are  solved  in  steady  state.  The  relevant  coupled 
field  equations  are  shown  and  discussed  in  Appendix  A.  Briefly,  in  the 
small  signal  limit  where  the  forward  pump  beam  (If)  and  the  backward 
pump  beam  (I^)  are  much  larger  than  both  the  amplified  signal  beam 
(Ig)  and  the  phase  conjugate  intensity  (I^)  throughout  the  entire 
sample,  the  four  pump  equations  (two  amplitude  equations,  two  phase 
equations)  become  decoupled  from  the  four  signal  equations.  An  analytical 
transcendental  expression  for  the  intensity  distribution  of  the  two 
counterpropagatlng  pump  beams  can  then  be  obtained  and  used  to  solve  the 
pump  phase  and  amplitude  equations  (Ref.  7-3).  Finally,  these  solutions 
for  the  spatial  dependence  of  the  pump  waves  are  used  to  determine  the 
spatial  dependence  of  the  nonconstant  coefficients  appearing  in  the  four 
coupled  signal  and  conjugate  equations.  These  four  equations  have  been 
solved  simultaneously  using  three  different  numerical  methods.  All  three 
methods  give  the  same  results  to  many  significant  figures.  Two  of  the 
methods  are  noniterative  in  nature,  and  rely  upon  the  fact  that  in  the  low 
signal  limit  the  DFWM  reflectivity  is  independent  of  the  incident  signal 
beam  intensity.  In  these  two  methods,  an  arbitrary  value  is  assigned  to 
the  amplified  transmitted  signal  beam  intensity  1^  at  z-L  (L  is  the 
length  of  the  nonlinear  medium,  see  Figure  A-1  in  Appendix  A).  The 
problem  simplifies  to  an  initial  value  problem  since  all  eight  variables 
are  now  known  at  the  same  location  (z— L).  The  equations  are  then 
integrated  from  z— L  to  z-0.  In  this  way  the  final  conjugate  intensity  and 
the  value  of  the  initial  signal  necessary  to  produce  such  a  conjugate 
intensity  are  obtained.  The  ratio  of  these  two  quantities  is  the 


reflectivity.  The  two  noniterative  methods  used  to  integrate  the 
equations  were  a  Runge  Kutta-Verner  fifth  and  sixth  order  method  and  a 
variable  order  Adams  predictor  corrector  method  (Ref.  7-4).  Identical 
results  were  obtained  with  both  methods  to  many  significant  figures. 

As  a  further  check  on  the  validity  of  the  numerical  calculations ,  a 
subroutine  capable  of  solving  a  system  of  differential  equations  with 
boundary  conditions  specified  at  two  different  points  was  empljyed.  The 
routine  uses  a  variable  order,  variable  step  size  finite  difference  method 
with  deferred  corrections  (Ref.  7-5).  As  mentioned  above,  all  three 
methods  gave  identical  results . 

The  plots  of  the  DFWM  reflectivity  versus  the  line  center 
absorption  coefficient  (Oq)  shown  in  Figure  7-1  are  numerical 
solutions  obtained  in  Ref.  7-3.  The  symbols  plotted  on  the  figure  are  the 
numerical  results  obtained  by  each  of  the  three  methods  described  above. 

In  most  regions  of  the  figure  the  two  results  agree  very  well.  However, 
discrepancies  exist  in  some  regions,  particularly  when  the  reflectivity 
becomes  large.  We  believe  these  differences  are  due  to  the  fact  that  the 
numerical  integration  used  in  Ref.  7-3  was  not  sufficiently  accurate  in 
this  region  (Ref.  7-6). 


7.1.2 


Lng  In  Resonant  or 


The  other  ID  DFWM  code  which  we  have  developed  was  designed  to  be  a 
precursor  of  the  two  dimensional  (2D)  code  so  that  all  of  the  formalism 
developed  for  the  ID  code  could  be  easily  transferred  to  the  2D  code. 

This  ID  code,  known  as  IDP,  has  a  time  dependent  formalism  which  enables 
it  to  do  both  pulsed  and  CW  calculations.  Eight  coupled  equations  for  the 
amplitude  and  phase  of  each  of  the  four  fields  are  solved  and  include  the 
effects  of  pump  depletion  and  absorption.  The  IDP  code  has  models  for  two 
types  of  FWM  media:  Kerr-like  and  resonantly  enhanced.  The  equations  for 
the  resonantly  nhanced  version  are  given  in  Appendix  A:  the  Kerr- like 
equations  are  similar  and  are  given  in  Ref.  7-7.  The  model  of  Kerr-like 
media  includes  nonsaturable  background  losses  (e.g.  scattering)  and  is 
valid  for  strong  signal  intensities.  The  model  of  resonantly  enhanced  FWM 
is  based  on  the  two-level  model  of  Brown  (Ref.  7-2)  and  includes  saturable 
absorption  and  dispersion.  For  the  resonant  FWM  media,  the  eight 
equations  solved  here  are  similar  to  those  discussed  above  in  Section 
7.1.1.,  but  a  different  method  of  solution  is  used.  Only  five  of  the 


eight  equations  need  to  be  solved  simultaneously  in  order  to  obtain  the 
magnitude  of  the  phase  conjugate  reflectivity.  However,  in  the  2D  DFWM 
model  (Section  7.2)  which  was  developed  at  TRW  in  order  to  examine  the 
influence  of  various  FWM  parameters  on  the  fidelity  of  the  phase 
conjugation,  it  is  necessary  to  solve  all  eight  Equations.  In  order  to 
enable  the  results  of  the  IDP  model  to  be  directly  transferable  to  the  2D 
model  all  eight  equations  are  solved  in  this  ID  case.  Similar  statements 
apply  to  the  Optical  Kerr  version  of  this  IDP  code  except  that  the 
equations,  taken  from  Ref.  7-7,  rigorously  include  the  effects  of  pump 
depletion.  In  contrast,  in  the  resonant  IDP  model  pump  and  signal 
absorption  are  qixantltatively  accounted  for,  whereas  the  treatment  of  pump 
depletion  is  only  qualitative  since  analytical  expressions  for  the 
coupling  coefficients  at  high  signal  Intensity  have  not  been  formulated 
and  numerical  evaluation  requires  prohibitive  computer  capacity.  The 
coupling  coefficients  for  high  signal  intensity  have  been  approximately 
formulated  for  the  case  of  cross  polarized  pumps ,  enabling  the  resonant 
IDP  code  to  give  qualitative  results  for  pump  depletion. 

Figure  7-2  shows  some  sample  results  for  an  optical  Kerr-like 
medium.  In  Part  A  of  the  figure  the  DFWM  is  performed  with  a  CW  laser 
with  a  pump  Intensity  of  100  MW/cm^  and  a  medium  length  of  10  cm.  The 
curves  in  the  figure  plot  the  intensities  of  the  four  light  beams  as  a 
function  of  position  in  the  FWM  medixun.  Under  these  intense  pumping 
conditions  the  gain  experienced  by  the  signal  and  the  conjugate  signal  is 
so  large  that  these  signals  become  comparable  in  magnitude  to  the  pump 
beams.  Due  to  conservation  of  energy  the  pump  beams  are  largely 
depleted.  In  part  B  of  the  figure  a  pulsed  laser  of  1  nsec  duration  (30 
cm  in  length)  is  used  to  perform  the  FWM.  A  square  pulse  shape  has  been 
used.  Once  gain,  large  signals  are  obtained  and  pump  depletion  occurs. 

The  length  of  the  FWM  medium  is  shorter  Chan  the  laser  pulse  lengths .  As 
a  result,  the  front  part  of  the  counterpropagating  pump  pulses  are 
depleted  while  the  rear  parts  are  not  because  they  experience  less  gain 
since  they  do  not  overlap  throughout  the  entire  medium.  This  effect  gives 
rise  to  the  odd  intensity  distribution  shown  in  the  figure. 

7.2  2D  CODE  FOR  CW  OR  PULSED  DFWM  IN  RESONANT  OR  KERR- LIKE  MEDIA 

Many  analytical  and  numerical  techniques  have  been  developed  in 
recent  years  to  study  optical  phase  conjugation  via  backward  four-wave 
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mixing.  So  far,  the  numerical  studies  are  limited  only  to  one -dimensional 
analysis  of  the  effect  in  steady  state  (Ref.  7-2,  7-3,  7-7).  Reflectivity 
is  the  only  resulting  quantity.  Other  important  features  such  as  beam 
misalignment,  self-focusing  and  phase  conjugation  fidelity  have  remained 
unexplored.  These  features  can  only  be  studied  by  including  at  least  one 
additional  transverse  dimension  to  account  for  the  intensity  and  phase 
variation  of  the  fields.  For  this  purpose,  a  two-dimensional  code 
simulating  four-wave  mixing  has  been  developed  at  TRW.  The  original 
version  of  the  code,  which  modeled  FWM  in  a  Kerr  medium,  was  developed  on 
the  Free  Electron  Laser  Resonator  Development  Program  (Ref.  7-8).  The 
extension  to  resonantly  enhanced  media  was  performed  on  NLOT.  The  code 
follows  the  actual  time  history  of  wave  Interaction  in  the  nonlinear 
medium  and  is  capable  of  simulating  the  mixing  process  with  both 
finite-length  pulses  and  CW  conditions. 

The  notation  used  here  for  the  backward  four-wave  mixing 
configuration  is  as  follows.  Beam  1  and  2  are  the  colinear  pump  beams 
propagating  in  opposite  directions.  Beam  3  is  the  signal  beam  incident  at 
a  very  small  angle  with  respect  to  beam  1.  The  plane  defined  by  these 
three  beams  is  the  two-dimensional  working  space  for  the  system,  and  it  is 
assumed  that  there  is  no  variation  of  physical  parameters  in  the  direction 
normal  to  this  plane.  Usually,  direction  z  is  defined  as  the  line 
bisecting  the  angle  between  beam  1  and  3  (Figure  7-3),  although  the 


Figure  7-3.  A  two-dimensional  schematic  system  of  four  wave  mixing. 


Incidence  angle  of  each  beam  can  be  input  Independently.  The  conjugate 
signal  (beam  4)  Is  expected  to  propagate  on  the  same  plane  and  In  nearly 
the  opposite  direction  of  beam  3. 

The  region  of  the  nonlinear  medium  extends  from  0  to  L  In  the  z 
direction  and  has  an  aperture  large  enough  to  cover  the  four  beams .  At 
the  beginning  (t  •  0) ,  the  pulse  front  advance  of  each  beam  Into  the 
medium  can  be  adjusted  to  account  for  pulse  delay  In  any  of  the  beams. 

For  longitudinal  sampling  the  pulses  are  sliced  Into  equally  spaced  sheets 
perpendicular  to  the  z-axls.  The  slice  distance  Is  twice  that  of  the 
propagation  distance  In  each  time  step,  so  that  the  field  on  each  sheet 
will  Interact  with  every  field  sheet  In  the  oppositely  directed  pulse. 

Transversely,  the  fields  are  sampled  In  equally  spaced  grid  points 
across  the  whole  aperture.  The  beam  area  Is  defined  by  either  a  Gaussian 
profile  or  an  apodlzed  (hypergausslan)  square  function.  The  apodlzed 
edges  reduce  the  Fresnel  ripple  In  propagating  through  the  cell.  Because 
the  field  sheets  are  not  perpendicular  to  the  actual  direction  of  beam 
propagation  for  non-zero  Incidence  angle,  a  linear  phase  variation  is 
applied  to  the  field  across  the  beam.  For  pulse  applications,  such  as 
those  calculations  performed  for  TRW's  Free  Electron  Laser  Resonator 
Development  Program  (FELRDP,  Ref.  7-8),  several  different  longitudinal 
pulse  profiles  can  be  used  (Figure  7-4).  These  pulsed  calculations  are 
discussed  here  in  order  to  Illustrate  the  capabilities  of  the  2D  code. 

Before  the  pulses  enter  the  interaction  region,  a  phase  screen 
representing  the  phase  aberration  can  be  applied  to  the  complex  field. 

The  phase  screen  consists  of  up  to  five  sinusoidally  varying  functions. 
With  different  variation  amplitudes  and  periods,  various  phase  screens  can 
be  created  Independently  for  each  beam.  The  phase  screen  of  beam  3  should 
be  applied  to  the  exiting  beam  4  to  check  how  well  the  aberration  is 
corrected.  In  each  case,  the  beam  quality  is  calculated  in  terms  of  the 
Strehl  ratio  calculated  from  the  rms  phase  variation. 

The  pulses  advance  In  the  medium  using  the  near- field  propagation 
scheme.  In  each  propagation  step,  the  one -dimensional  field  array  is 
Fourier- transformed,  multiplied  by  a  transfer  phase  function  and  then 
inverse  Fourier  transformed.  For  CW  cases,  the  propagation  continues 
until  a  steady  state  is  reached.  For  pulse  interaction,  the  propagation 
stops  when  the  rear  end  of  every  pulse  leaves  the  medium  completely.  At 
the  end  of  the  run,  the  reflectivity  and  Strehl  ratio  are  calculated.  The 
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field  arrays  are  stored  for  three-dimensional  plots  of  pulse  intensity  and 
phase  profiles. 

In  the  nonlinear  medium,  every  field  interacts  with  itself  and  with 
each  other  through  the  nonlinear  optical  susceptibility.  The  code  models 
such  Interaction  in  Kerr-like  media  as  well  as  in  saturable  absorbing  (or 
gain)  materials.  In  a  Kerr-like  medium,  the  nonlinearity  can  be  described 
by  an  effective  index  of  refraction  n: 


n  -  n^  +  n2  <E  >. 


(7-1) 


where  <E^>  is  the  time  average  of  the  square  of  the  total  field.  Before 
each  propagation,  the  complex  field  arrays,  Aj  (n^,  n^),  are  updated 
according  to  the  formula 


^3  ■  •*  [(2'o  -  *  “JVn] 
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where  the  Indices  are 
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(7-5) 


Iq  -  Ii  +  l2  +  I3  ^4 


is  the  vacuum  wave  number,  Is  the  Incident  angle  of  beam  j  and 
Az  is  the  propagation  distance. 

Equation  7-2  can  also  be  tranformed  from  complex  field  amplitudes 
Into  coupled  Intensity  and  phase  equations.  Both  methods  are  available  in 
the  code.  Althou^  using  coupled  Intensity  equations  guarantees  energy 
consejTvation,  one  has  to  be  very  careful  in  the  region  where  the  intensity 
Is  zero  because  the  Intensity  variable  appears  In  the  denominator  of  the 
accompanying  phase  equations.  In  general,  the  first  term  In  Equation  7-2 
represents  the  self  phase  modulation  which  would  result  In  self- focusing 
for  a  nonuniform  beam.  The  second  term  represents  the  four  wave 
Interaction  which  would  generate  the  conjugated  return  signal. 

For  the  saturable  absorbers,  the  model  follows  the  equations  given 
in  Appendix  A.  Although  the  coupled  equations  were  derived  In  the  weak 
signal  limit,  we  have  added  depletion  terms  In  the  equations  for  the  pump 
beams.  In  this  way,  the  pump  depletion  effect  for  the  four  wave  mixing  In 
saturable  absorbers  Is  qualitatively  observed. 

In  practically  all  FWM  experiments  Gaussian  laser  beams  are  used, 
whereas  most  theoretical  treatments  of  FWM  assume  that  plane  waves  are 
employed.  In  these  FWM  experiments  with  Gaussian  beams  different  results 
are  often  obtained  than  those  predicted  by  the  plane -wave  FWM  theories. 
These  effects  will  be  particularly  Important  in  saturable  FWM  media,  where 
the  magnitude  of  saturation  will  vary  transversely  across  the  Gaussian 
pump  beams .  The  2D  FWM  code  is  very  useful  in  examining  such  effects . 
Using  the  2D  computer  code,  the  following  example  calculations  were 
performed.  In  the  first  calculation,  the  pump  Intensities  are  1/1"^^^ 

—  2,  the  detuning  Is  zero,  the  absorbance  is  three,  and  the  spot  size  of 
all  three  input  beams  is  200  /im  (radius  of  the  1/e  point  of  the 
E-field).  For  these  conditions  the  DFWM  reflectivity  is  equal  to  1.8%. 

In  a  second  calculation  where  all  the  parameters  are  Identical  except  the 
signal  beam  Is  made  smaller  with  a  spot  size  equal  to  50  pm,  the 
reflectivity  decreases  by  -25%  to  1.4%.  For  plane  waves,  the  DFWM 
reflectivity  Is  Independent  of  the  signal  beam  spot  size.  The  decrease  in 
reflectivity  observed  when  the  signal  beam  spot  size  Is  made  substantially 
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smaller  than  those  of  the  pump  beams  is  due  to  a  reduced  DFWM  reflectivity 
in  the  highly  saturated  central  position  of  the  pump  beams.  The  average 
I/Isat  is  4.3  times  larger  in  the  central  25%  of  the  pump  beams  where 
the  50  fim  signal  beam  overlaps,  than  it  is  in  the  outer  25%  of  the 
Gaussian  pump  beam.  Thus,  the  probe  beam  with  a  spot  size  matched  to 
those  of  the  pump  beams  can  sample  regions  of  the  pumped  FWM  medium  which 
are  not  as  severly  saturated  and  can  extract  additional  conjugate 
intensity  from  these  regions. 

Another  consequence  of  the  relative  signal  beam  spot  size  compared 
to  those  of  the  pump  beams  can  be  seen  in  Figure  7-5  where  the  transverse 
intensity  profile  of  the  phase  conjugate  is  plotted  for  the  two  cases  just 
discussed.  When  the  signal  beam  is  much  smaller  than  the  pump  beams 
(Figure  7-5A),  the  pump  intensity  varies  little  across  the  signal  beam  and 
the  situation  approximates  that  of  plane-wave  pvimp  beams.  In  this  case 
since  the  reflectivity  is  independent  of  the  input  signal  intensity  (no 
pump  depletion) ,  the  signal  conjugate  is  seen  to  be  a  Gaussian  with  the 
same  spot  size  as  the  incident  signal  beam.  In  contrast,  when  the  signal 
beam  is  Identical  in  size  to  the  pump  beams,  the  phase  conjugate 
transverse  intensity  profile  is  highly  non-Gaussian  (Figure  7-5B).  The 
latter  effect  has  been  predicted  theoretically  (Ref.  7-9). 

The  particular  Gaussian  beam  results  discussed  here  apply  only  to 
the  specific  conditions  listed  above.  When  I  «  Iggi-  or  in  regions 
where  the  change  in  reflectivity  with  pump  intensity  is  relatively  flat 
the  variation  in  the  signal  spot  size  will  not  have  an  effect  on  the  phase 
conjugate  reflectivity.  However,  as  TRW  is  primarily  interested  in  phase 
conjugation  in  the  high  reflectivity  regime,  the  Gaussian  beam  effects 
predicted  here  by  the  2D  model  can  have  important  consequences  on  the 
results  of  the  FWM  experiments. 

The  most  advantageous  property  of  a  two-dimensional  code  over  a 
one -dimensional  code  is  its  ability  to  demonstrate  the  fidelity  of  the 
phase  conjugation  process  in  the  four  wave  mixing  process.  In  the 
following  example,  calculated  for  FELRDP  (Ref.  8) ,  2  mm  long  uniform 
pulses  are  used  for  all  three  input  beams  (Figure  7-6A),  A  phase  screen 
prepared  with  five-term  sinusoidal  functions  is  applied  to  the  signal  beam 
before  entering  the  medium.  The  aberrated  beam  is  calculated  to  have  a 
Strehl  ratio  of  0.77  (Figure  7-6B). 

The  intensity  profile  of  the  return  pulse  has  a  higher  trailing  end 
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Figure  7-6.  Dt.,.onstration  of  phase  conjugation  in  four  wave  mixing. 

(A)  Input  signal  pulse  profile,  (B)  Phase  of  aberrated 
signal  beam,  (C)  Return  signal  profile,  (D)  Cot^rected 
phase  distribution  in  the  return  signal. 
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(Figure  7-6C) .  This  can  be  understood  because  the  energy  transfer 
®^^i-clency  Is  proportional  to  the  square  root  of  the  Intensity  product 

AI  - 

Since  (Ij^  +  Ij)  and  (I2  +  14^)  are  constants  of  the  Interaction, 
the  Intensity  product  will  be  higher  when  I3  and  Increase, 

Therefore,  the  energy  transfer  Is  more  efficient  near  the  end  of  the 
pulses.  The  ripples  on  the  top  of  the  Intensity  profile  are  due  to  the 
applied  phase  screen. 

The  phase  screen  Is  applied  to  the  conjugated  return  beam  4  after 
exiting  from  the  medium.  The  phase  of  the  field  becomes  uniform  In  most 
of  the  pulse  area  (Figure  7-6D).  The  phase  Increase  near  both  sides  of 
the  pulse  Is  the  result  of  wave  propagation  and  Is  In  a  very  low 
Intensity  area.  The  correction  of  the  beam  aberration  Is  nearly  complete, 
the  Strehl  ratio  Is  calculated  to  be  0.997. 
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APPENDIX  A.  THEORY  FOR  DFWM  CODES 


The  four-wave  mixing  system  Is  shown  In  Figure  A-1.  The  two  pump 
beams  are  the  forward  ("f")  and  backward  ("b")  beams,  which  propagate  along 
the  Zq  direction.  The  signal  beam  ("s")  and  Its  conjugate  ("c")  are 
counterpropagatlng  along  the  z'  direction.  The  field  equations  use  the  z 
coordinate  through  the  transformations  z^  -  z  cos(tf/2)  and  z'  -  z 
cos(tf/2) . 

Our  starting  point  for  the  Intensity  and  phase  equations  for  each  beam 
used  In  the  code  are  the  four  coupled  non- linear  differential  equations  for 
each  field  derived  In  Boyd  et  al.  (Ref.  A-1).  They  start  from  the  Maxwell 
equations  satisfied  by  the  total  field  as  the  source.  (Polarization  Is 
solved  as  a  steady  state  solution  In  the  Bloch  equation  of  the  two- level 
system  medium.)  With  the  assumptions  of  a  strong  pvimp,  weak  signal,  and  the 
slowly-varylng-envelope  approximation  for  the  field  amplitude,  they  arrived 
at  the  following  field  equations: 


3Z  ®f  *  *®f®f 

3z  '  ■^“b^b 


(A-1) 


(A-2) 


(A-3) 


(A-4) 


where 

“f.b  ■  “o  (1  * 

f  ,b 


a  »  a^(l  +  160)12)  (B/C^) 
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our-wave  mixing  system  in 


“o  =  %/il  +  6<oT2)  e  a^/Cl  +  6^) 


a  being  the  line-center  field  absorption  coefficient,  which  becomes  the 


gain  coefficient  if  it  is  less  than  0,  and 


'f.b= 


B  =  1  +  Ifa  + 


C  =  {B^  - 


1/2 


=  {1  ^  (If  +  ib)  +  (If  -  ibrl 


2,1/2 


where  Eg  is  the  medium  saturation  field  amplitude. 

From  the  equations  A.l  through  A. 4,  we  have  eight  coupled  equations  of 
intensity  and  phase  for  all  four  beams.  The  intensities  and  phases  are 
defined  by: 


:  _  =  p  p  ^ I ^f  K ) 

f,b  '^f.b  ® 


e  =  e  I *^5  r ^ 

S,c  ®s,c  ® 
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•  c  S  ~  Is.c'^^sat 


where,  for  a  -  f,b  or  s,c, 
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The  eight  coupled  equations  are: 
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where 


“l  ■  2-0 

V  =  4a^j/C^ 

♦  =  *5  +  ^c  ■  *f  '  *b 
U  =  2a^  B/C^  »  U4 

with 


o  ,  6,  B,  and  C  defined  above. 

Q  »  »  » 


Equations  A- 5  through  A- 12  are  truely  equivalent  to  A-1  through  A- 4  if  one 
omits  the  underlined  terms  in  equations  A-5  through  A-8.  These  terms  are 
added  rather  than  derived  from  the  field  equations.  This  is  an  approximate 
way  to  account  for  the  pump  depletion  following  the  derivation  of  similar 
equation  for  Kerr-like  media  by  Hardy  and  Sllberberg.  (Ref.  A-2). 

For  application  to  two-level  systems,  inputs  needed  are  resonant  wave 
length,  pump  laser  wave  length,  line  center  field  absorption  (gain) 
coefficient  •  gain  medium  saturation  intensity  Igat’ 

normalized  detuning  parameter  S.  At  resonance  £  -  0,  hence  our  results 
are  not  affected  by  a  change  of  sign  in  (  .  This  code  also  applies  to 
Kerr- like  media  if  one  inputs  the  susceptibility. 

Since  the  derivation  is  based  on  the  assumption  of  a  weak  signal,  a 
more  rigorous  treatment  is  needed  for  exploring  the  strong  signal  regime. 
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